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Foreword 

Dr.  Edward  S.  Epstein 

Director,  National  Climate  Program  Office 

National  Oceanic  and  Atmospheric  Administration 

The  occasion  of  the  Third  Climate  Diagnostics  Workshop  signaled  as 
well  the  beginning  of  the  National  Climate  Program  under  Public  Law  95-367. 
The  National  Climate  Program  Act,  signed  into  law  on  September  17,  1978, 
encourages  development  of  national  applications  and  services  of  climate 
information  --  information  derived  from  active  research.  These  annual 
Workshops  serve  both  research  and  services  in  this  program.  They  are, 
first  and  foremost,  intended  to  review  the  climate  of  the  preceding  year 
and  examine  particularly  unusual  recent  climatic  events.  This  review 
affords,  too,  an  opportunity  to  examine  diagnostic  tools  and  techniques, 
to  exchange  current  impressions  on  the  mechanisms  behind  short-term  cli- 
mate variations,  and  to  introduce  more  scientists  and  institutions  to 
the  study  of  climate. 

Of  perhaps  particular  interest  to  those  who  assembled  in  Miami  and 
those  who  will  find  these  Proceedings  useful,  1978  also  marked  the 
establishment  of  the  NOAA  Climate  Analysis  Center  in  the  National  Weather 
Service;  some  of  the  concepts  of  tasks  for  this  Center  came  from  previous 
Workshops,  and  conversely,  some  Workshop  concepts  from  early  planning  for 
the  Center.  The  Climate  Analysis  Center  should  become  a  major  bridge  be- 
tween research  and  services  under  the  National  Program.  The  CAC  has  been 
established  to:  maintain  awareness  of  current  climate  and  weather 
anomalies  and  fluctuations;  produce  climate  projections  and  extended  range 
outlooks;  disseminate  climate  information  to  support  applications,  includ- 
ing impact  assessments  and  research;  and  develop  and  test  improved  diag- 
nostic and  climate  projection  techniques  and  products. 

Thus,  these  Workshops  have  provided  those  of  us  planning  the  National 
program  with  an  idea  of  the  limits  of  our  abilities  and  what  is  needed  to 
improve  them.  Views  expressed  by  participants  of  previous  Workshops  have 
improved  our  plans  and  helped  develop  National  Program  concepts.  As  the 
National  Climate  Program  continues  to  evolve,  planners  will  draw  upon 
these  Workshops  again  for  guidance,  for  information  and  for  furthering  co- 
operative looks  at  the  global  climate.   I  hope  other  participants  will  find 
these  meetings  even  more  advantageous,  and  that  individual  and  collective 
research  aims  are  advanced  accordingly. 

I  wish  to  take  this  opportunity  to  thank  Dr.  Eric  Kraus,  Director  of 
the  Cooperative  Institute  for  Marine  and  Atmospheric  Sciences,  and  offi- 
cials of  the  University  of  Miami,  for  hosting  the  meeting  in  such  excellent 


facilities  and  for  helping  to  organize  the  Workshop.  My  appreciation 
is  extended  to  Ms.  Josie  Sanchez  of  CIMAS  who  served  ably  as  the  key  for 
local  arrangements  and  logistical  support.  Of  course,  the  success  of  a 
meeting  depends  on  those  who  come  to  participate.  I  have  had  an  oppor- 
tunity to  review  the  contents  of  these  Proceedings  and  find  this 
collection  a  stimulating  reminder  of  the  excellent  work  that  is  being 
done  and  of  the  science  we  can  look  forward  to.  My  thanks  to  all 
participants  for  their  contributions  which  have  made  this  a  valuable 
and  successful  meeting. 
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Preface 


The  "notes"  which  comprise  these  Proceedings  are  to  guide  the 
reader  to  current,  ongoing  research,  or  to  notions,  questions  or  problems 
that  beg  elaboration  or  solution.  Some  editorial  license  has  been  taken 
where  a  few  papers  have  been  moved  out  of  the  sequence  in  which  they  were 
given  at  the  Workshop;  this  to  improve  continuity.  Notes  on  a  study  by 
Quiroz  have  been  added  to  augment  information  on  the  most  recent  climate 
and  to  follow  up  on  many  participants'  interests  in  the  unusual  1976-77 
winter,  as  expressed  during  both  1977  and  1978  Workshops.  Interest  in 
the  winter  season  continues  (Wagner,  p.  1-1;  Diaz  and  Quayle,  p.  3-1; 
Dickson,  P.  4-1;  Quiroz,  p.  5-1;  Namias,  p.  6-1;  Seckel  ,  Bakun,  and 
Parrish,  p.  7-1;  Winston  and  Krueger,  p.  8-1,  Berg  and  Matson,  p.  10-1; 
Douglas,  p.  28-1,  and  Barnett,  p.  33-1)  because  of  the  potentially  great 
impact  on  human  activities  and  because  it  is  an  extreme  within  the  annual 
cycle  which  serves  as  an  important  point  of  calibration  in  testing  hypoth- 
eses of  climate  variability. 

Data  for  1978  have  been  examined  by  several  of  the  group  assembled, 
providing  views  of  the  most  recent  seasonal  climate:  Wagner  (p.  1-1); 
Angell  (p.  2-1);  Diaz  and  Quayle  (p.  3-1);  Dickson  (p.  4-1);  Quiroz 
(p.  5-1);  Namias  (p.  6-1);  Seckel,  Bakun  and  Parrish  (p.  7-1);  Winston 
and  Krueger  (p.  8-1);  Kukla  and  Gavin  (p.  9-1);  Berg  and  Matson  (p.  10-1); 
Douglas  (p.  28-1)  and  Ernst  and  Sherman  (p.  36-1).  Other  analyses  of 
atmospheric  data  reveal  a  biennial  variation  in  surface  temperatures  over 
the  U.S.  (Rasmusson,  et.  al .  p.  22-1)  and  suggestions  of  latitudinal,  inter- 
annual  continuity  in  700  mb  anomalies  (Hanson,  p.  23-1).  Large-scale  at- 
mospheric moisture  fluxes  are  examined  by  Rosen  (p.  24-1). 

Considerable  attention  was  given  during  the  Workshop  to  ocean-related 
climate  characteristics.  Miyakoda  and  Rosati ,  prompted  by  last  year's  Work- 
shop discussions  of  incomparability  between  sea  surface  temperature  data 
sets,  analyzed  ship  and  satellite  derived  data  over  Pacific  regions  in  an 
attempt  to  clarify  differences  (p.  11-1).  Variability  of  sea  surface 
temperature  was  examined  by  Merle  (p.  12-1);  Thang  and  Lamb  (p.  13-1); 
Brown  and  Evans,  who  also  compared  ship  and  satellite  measurements,  (p.  14-1); 
Steiner,  who  has  looked  for  possible  interrelationships  between  SST,  wind 
and  clouds  (p.  25-1);  Fletcher  (p.  26-1);  Starr,  Bryson  and  Kutzbach  (p.  27-1); 
Douglas  (p.  28-1),  and  Barnett  (p.  33-1).  Heat  transport  in  the  oceans  was 
reviewed  by  Oort  (p.  18-1),  estimated  by  Hastenrath  in  relation  to  large- 
scale  atmosphere  and  ocean  heat  budgets  (p.  21-1),  and  examined  through 
modeling  by  Bryan  (p.  19-1)  and  Schopf  (p.  20-1).  Heat  exchanges  at  the 
ocean-atmosphere  interface  were  calculated  by  Clark  (p.  16-1)  and  Bunker 
(p.  17-1). 
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Connections  between  terrestrial  hydrology  and  ocean  thermodynamics 
and  the  potential  for  affecting  regional  climate  were  discussed  by  Rooth 
(p.  15-1).  A  hydrologic  model  used  to  examine  moisture  exchanges  within 
the  climate  system  was  discussed  by  Eagleson  (p.  29-1),  and  Schaake  (p.  31-1) 
outlined  the  state-of-the-art  in  hydrologic  modeling  for  streamflow  pre- 
dictions using  climate  outlooks  prepared  by  the  National  Weather  Service. 
Although  not  formally  presented,  a  paper  has  been  added  to  the  proceedings 
(submitted  for  the  meeting  by  Kuril  ova,  et.  al . ,  P.  30-1)  in  response  to 
a  view  of  many  participants  that  insufficient  attention  has  been  placed 
on  understanding  the  hydrologic  cycle  as  part  of  the  climate  variability 
problem.  Robert  Hirsch,  in  material  presented  during  the  poster  session 
(not  formatted  for  inclusion  in  the  Proceedings)  discussed  activities  of 
the  U.S.  Geological  Survey  and  the  monthly  publication  of  "Water  Resources 
Review,"  which  reports  (among  other  items)  streamflow  and  ground-water 
conditions  in  the  U.S.  and  parts  of  Canada;  the  "Review"  is  free  on  ap- 
plication to  the  Water  Resources  Review,  U.S.  Geological  Survey,  Reston, 
Virginia  22092. 

Development  of  credible  climate  prediction  techniques  is  an  aim  of 
the  National  Climate  Program,  and  an  aim  of  several  of  the  Workshop  partici- 
pants. Among  those  who  addressed  the  prediction  aspects  of  climate  research 
were  Schaake  (p.  31-1),  Dittberner  (32-1),  Barnett  (33-1),  and  Schofield 
and  Barnett  (34-1),  while  Pittock  (p.  35-1)  ably  demonstrated  why  caution 
is  required. 

As  a  final  note,  the  assistance  of  Drs.  D.  Gilman,  R.  Harnack, 
S.  Hastenrath,  E.  Kraus,  C.  Rooth  and  J.  Winston  as  discussion  leaders  is 
gratefully  acknowledged.  Thanks  are  extended  to  all  Workshop  participants 
for  their  excellent  presentations,  lively  discussions  and  timely  insights 
into  climate  processes  and  prediction  research. 

William  A.  Sprigg 
November  1 978 
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A  REVIEW  OF  THE  SEASONAL  CIRCULATION  AND 
WEATHER  FROM  FALL  1977  THROUGH  SUMMER  1978 

A.  James  Wagner 
National  Weather  Service,  NOAA 

FALL  1977  -  Except  for  a  small   area  of  weak  blocking  north  of  Siberia, 
700  mb  heights  were  generally  below  normal   at  high  latitudes,  and 
fast  zonal  westerlies  prevailed  at  middle  latitudes  around  most  of 
the  hemisphere.     An  exception  was  eastern  North  America,  where  there 
was  a  tendency  for  recurrent  blocking. 

The  negative  height  anomaly  near  Vancouver  Island  was  associated  with 
frequent  storminess  in  the  area  and  the  advection  of  cool   maritime  air 
from  the  Gulf  of  Alaska  into  the  Pacific  Northwest  and  adjacent  areas 
of  Canada.     Most  of  the  southern  two-thirds  of  the  United  States  had 
warmer  than  normal   temperatures,  a  continuation  of  the  previous  summer 
pattern.     The  temperature  anomaly  of  +4°F  in  west  Texas  was  due  in  part 
to  a  heat  wave  which  produced  record  temperatures  well   in  excess  of 
100°F  during  the  last  third  of  September. 

Precipitation  was  above  normal  over  much  of  the  country,  the  only 
exceptions  being  part  of  the  Southwest,  where  a  severe  drought  con- 
tinued, and  the  western  portions  of  the  southern  and  central  Great 
Plains.     The  severe  summer  drought  over  much  of  the  Southeast  was 
relieved.     Heavy  rains  from  Hurricane  Babe  during  the  second  week  of 
September  fell   from  the  Louisiana  Coast  northward  and  eastward  to  the 
Tennessee  Valley  and  lower  Appalachians.     The  drought  in  the  mid- 
Atlantic  states  was  not  relieved  until  October,  however. 

WINTER  1977-78  -  Many  aspects  of  the  mean  circulation  were  similar  to 
that  of  the  previous  record  cold  winter,  although  subtle  differences 
between  the  two  patterns  led  to  most  of  the  records  in  the  1977-78 
season  being  due  to  heavy  snow  rather  than  cold.     A  large  area  of  nega- 
tive height  anomaly  again  covered  much  of  the  north-central   Pacific, 
and  an  extensive  area  of  blocking  was  centered  over  Canada.     A  deep 
trough  over  the  eastern  United  States  was  associated  with  frequent  and 
often  intense  storminess  contributing  to  record  and  near-record  snow- 
fall  totals  over  the  eastern  third  of  the  United  States. 

One  of  the  important  differences  between  the  past  two  winters  is  that 
the  strongest  westerlies  were  farther  south  in  1977-78,  with  the  low 
latitude  portion  of  the  western  North  American  ridge  weaker  than 
normal.     Instead  of  being  deflected  northward  toward  Alaska  as  in  the 
previous  winter,  Pacific  storms  entered  California,  giving  more  than 
twice  the  normal  winter  precipitation  to  much  of  the  Southwest.     The 
drought  was  decisively  broken,  but  at  the  cost  of  flooding  and  mud- 
slides in  southern  California.     The  main  positive  700  mb  height  anomaly 
was  centered  over  northern  Canada  instead  of  British  Columbia  as  in 
the  previous  winter.     Easterly  anomalous  flow  over  the  northern  Great 
Plains  contributed  to  upslope  flow  of  Arctic  air,  with  the  result  that 
cold  and  snowy  conditions  prevailed,  in  sharp  contrast  to  generally  mild, 
dry  weather  the  previous  winter. 
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Mild  Pacific  air  dominated  all  areas  west  of  the  Continental   Divide, 
where  storm- related  cloudiness  also  prevented  the  normal   nocturnal 
cooling  over  the  Plateau.     All  areas  east  of  the  Divide  were  colder 
than  normal,  but  east  of  the  Mississippi   River  it  was  not  quite  as 
cold  as  the  previous  winter.     Much  of  the  Great  Plains  was  considerably 
colder  in  1977-78.     This  was  due  in  part  to  the  eastern  North  American 
trough  being  farther  west  in  the  more  recent  winter.     This  position 
also  brought  more  moisture  into  the  East,  where  the  previous  winter 
had  been  quite  dry.     Extensive  areas  had  two  to  more  than  three  times 
their  normal  winter  snowfall   in  1977-78.     Record-breaking  blizzards 
paralyzed  the  Midwest  in  January  and  southern  New  England  in  February. 

SPRING  1978  -  A  wave  number  two  pattern  at  high  latitudes  was  associated 
with  blocking  over  the  Beaufurt  Sea  and  Northern  Scandinavia  while  an 
intense  vortex  was  located  over  Baffin  Island.     A  weak  area  of  blocking 
was  centered  just  north  of  Minnesota.     This  feature,  together  with  the 
absence  of  appreciable  troughing  in  the  West,  produced  a  generally 
milder  than  normal   spring  over  much  of  the  West.     Fast  mid-latitude 
westerlies  prevailed  over  both  oceans,  with  deeper  than  normal  Aleutian 
and  Icelandic  lows.     A  stronger  than  normal   trough  over  Europe  contri- 
buted to  an  unusually  cold  and  wet  spring  there. 

In  contrast  to  the  fast  mid-latitude  westerlies  over  the  oceans,  below 
normal   heights  across  the  southern  United  States  and  the  Canadian  block 
meant  weak  westerlies  in  that  area.     Storms  crossing  the  U.S.  moved 
slowly  and  produced  large  quantities  of  precipitation.     Pacific  storms 
continued  to  produce  more  than  twice  the  normal   seasonal   rainfall   along 
the  California  coast.     The  only  noticeably  dry  area  was  in  Texas  and 
Louisiana,  where  some  areas  had  less  than  half  the  normal   precipitation. 

SUMMER  1978  -  The  strong  polar  low  persisted  over  northern  Canada  into 
the  summer,  moving  slightly  westward  from  its  spring  position.     It  was 
surrounded  by  a  ring  of  positive  height  anomalies,  south  of  which  an 
active  storm  track  prevailed  across  northern  Eurasia.     The  area  of  posi- 
tive height  anomalies  broadened  to  include  middle  latitudes  from  the 
central   Pacific  to  the  central  Atlantic,  with  several  moderately  strong 
anomaly  centers  across  the  United  States  from  lower  California  to  the 
mid-Atlantic  States.     Each  was  associated  with  a  center  of  positive 
temperature  anomaly  and  at  times,  drought. 

The  overall   temperature  anomaly  pattern  was  similar  to  that  of  the  pre- 
ceding summer,  although  areas  of  heat  and  drought  were  for  the  most  part 
much  less  severe.     The  seasonal   mean  maps  do  not,  however,  reflect  the 
fact  that  parts  of  central  Texas  had  a  record-breaking  heat  wave  during 
July,  with  readings  as  high  as  114°F  at  Wichita  Falls.     During  the  first 
week  of  August,  torrential   rains  from  moisture  brought  inland  from 
Tropical   Storm  Amelia  at  the  end  of  July  caused  widespread  flooding  along 
a  north-south  line  through  central   Texas.     Some  totals  were  in  excess  of 
20  inches.     By  the  end  of  the  summer,  rather  severe  drought  conditions 
existed  over  Oklahoma,  northeast  Texas,  and  portions  of  adjacent  states. 
Sections  of  the  Southeast  and  Southwest  were  also  dry,  but  drought  con- 
ditions in  these  areas  were  less  severe  than  during  the  previous  summer 
because  fairly  good  rains  had  fallen  in  the  preceding  months  of  1978. 
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An  Update  of  Temperature  and  Vortex-Area  Data 
Through  the  Winter  of  1977-78 

J.K.  Angell 
Air  Resources  Labs,  ERL,  NOAA 
Silver  Spring,  MD  20910 


The  following  are  the  main  results  obtained  from  an  update,  through 
the  northern  winter  of  1977-78,  of  global  temperature  variations  at 
the  surface  and  in  surface-100  mb  and  100-30  mb  layers,  as  well  as  from 
an  update  of  the  size  and  location  of  the  300  mb  north  circumpolar  vor- 
tex.  Some  other  points  of  interest  are  also  recorded.   Figure  1  shows 
the  global  distribution  of  radiosonde  stations  used  initially  for  the 
temperature  analysis.   In  recent  years  some  substitutions  have  had  to 
be  made  in  the  equatorial  zone. 

1.  The  large  3-yearly  temperature  oscillations  in  the  tropics, 
apparently  related  to  El  Nino  (EN)  occurrences,  have  recently  decreased 
in  amplitude  and  become  less  uniform  (Figs.  2  and  3). 

2.  There  is  evidence  of  a  similar  3-yearly  temperature  oscillation 
in  the  surface-100  mb  layer  of  north  polar  latitudes,  almost  in  phase 
with  the  oscillation  in  the  tropics  (Fig.  3). 

3.  In  recent  years  there  has  been  a  greater  temperature  decrease 
in  north  temperate  latitudes  than  in  north  subtropics,  resulting  in  an 
increase  in  meridional  temperature  gradient  between  the  two  zones  (Figs . 
2  and  3). 

4.  In  north  temperate  latitudes  the  overall  decrease  in  temperature 
has  been  associated  with  an  overall  increase  in  total  ozone.   On  the 
shorter  time  scale  there  is  evidence  of  a  lag,  i.e.,  the  relatively 

high  total  ozone  values  in  1977  appear  to  be  related  to  the  cold  year 
1976  (Fig.  3). 

5.  In  the  Northern  Hemisphere  since  about  1965  the  surface  temp- 
erature has  tended  to  increase  slightly  while  the  surface-100  mb  temp- 
erature has  tended  to  decrease  slightly,  implying  an  increase  in  lapse 
rate  (Fig.  4). 

6.  In  the  100-30  mb  layer  of  the  Northern  Hemisphere,  the  tempera- 
tures were  relatively  high  following  the  eruptions  of  Agung  and  Fuego 
(Fig.  4),  but  for  the  world  as  a  whole  the  temperatures  were  much  higher 
following  Agung  (Fig.  5). 

7.  The  surface  temperature  decrease  following  Agung  was  larger 
in  the  Northern  Hemisphere  than  for  the  world,  but  note  that  in  both 
cases  the  decrease  was  not  unique  (Figs.  4  and  5). 
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8.  The  large  interannual  variability  of  temperature  in  recent 
years  is  shown  by  the  observation  that  1976  (X)  was  a  record  cold  year 
globally,  whereas  1977  (XX)  was  a  relatively  warm  year  globally  (Fig.  6). 
This  alteration  was  worldwide,  being  found  in  the  tropics  and  both 
extratropical  regions. 

9.  In  north  temperate  latitudes  the  winter  of  1977-78  (XX)  was 
indicated  to  be  not  as  cold  as  the  winter  of  1976-77  (Fig.  7),  even 
though  the  area  of  the  300  mb  north  circumpolar  vortex  was  indicated  to 
be  greater  in  the  winter  of  1977-78  (Fig.  8). 

10.  As  perhaps  a  contributing  factor  to  the  discrepancy  in  9  above, 
it  has  been  found  from  the  70-year  Historical  Weather  Map  Series  that 
the  central  latitude  of  the  subpolar  low-pressure  belt  and  subtropical 
high-pressure  belt  have  averaged  about  1  degree  of  latitude  less  at  the 
time  of  11-year  sunspot  maximum  than  sunspot  minimum  (Fig.  9).   There  is 
evidence  also  of  a  22-year  cycle  in  the  central  pressure  of  subpolar  low 
and  subtropical  high  (Fig.  9). 

11.  The  eccentricity  of  the  north  circumpolar  vortex  was  the  smal- 
lest of  record  (vortex-center  nearest  the  North  Pole)  during  1977  Fig.  8). 

12.  In  the  25-55  km  layer  of  the  western  hemisphere  the  temperature 
is  indicated  to  have  decreased  by  about  3°C  between  sunspot  maximum  in 
1969  and  sunspot  minimum  in  1976  in  polar  (P) ,  subtropical  (S)  and  equa- 
torial (E)  latitudes  (Fig.  10) .   It  will  be  of  interest  to  see  if  this 
temperature  increases  again  with  the  recent  upsurge  in  solar  activity. 

13.  There  is  evidence  of  about  a  2-season  lag  between  the  3-yearly 
temperature  oscillation  in  the  equatorial  zone  and  in  both  temperate 
latitudes  (Fig.  11).   Hence,  if  the  tropical  oscillation  persists,  there 
would  appear  to  be  some  possibility  of  foreshadowing  seasonal  temperatures 
in  the  extratropics.   On  the  basis  of  the  somewhat  limited  data  record, 
summer  temperatures  are  better  foreshadowed  than  winter  temperatures  in 
this  way. 

14.  The  warm  periods  in  the  3-yearly  tropical  temperature  oscilla- 
tion have  been  associated  with  enhancements  in  the  rate  of  increase  of 
carbon  dioxide  (Fig.  12). 
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Fig.    8.      Comparison  with  other  years   of    the   area,    eccentricity,    and 
displacement  along   axes   90°W-90°E   and  0°-180°    of   the  300  mb   north  cir- 
cumpolar  vortex  during  1976    (X)    and   1977    (XX),    and   at   right   the   compari- 
sons  for   the  northern  winter  of   1976-77    (X)    and   1977-78    (XX) .      A  positive 
eccentricity  signifies   that   the  vortex  center   is   further   from  the  North 
Pole   than  average  and   similarly,    for   the  bottom  diagrams,    that   the  vortex- 
center   is  displaced   further   in  the  direction  of   90°E  and  180°    than  average. 
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WESTERN  HEMISPHERI 
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1975 


Fig.  10.   Temperature  trends  as  a  function  of  height  for  north  polar  (P) , 
north  subtropical  (S)  and  equatorial  (E)  zones  of  the  western  hemisphere. 
A  1-2-1  smoothing  (1-1  at  beginning  and  end  of  record)  has  been  applied 
twice  to  successive  year-average  values. 
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Fig.  11.   Recent  temperature  variations,  for  the  surface-100  mb  layer,  in 
tropical  and  north  and  south  temperate  latitudes,  where  again  a  1-2-1 
smoothing  has  been  applied  twice  to  successive  seasonal  values.   The 
dashed  lines  connect  corresponding  temperature  maxima  and  minima. 
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Fig.  12.   Comparison  of  the  rate  of  increase  of  carbon  dioxide  at  Mauna 
Loa  and  the  South  Pole  with  the  temperature  variations  in  tropics  and 
extratropics.   Smoothing  as  in  Fig.  2. 
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Contribution  to  the  1978  Climate  Diagnostics  Workshop 

A  Comparative  Analysis  of  the  1977-78  Winter  Season  in  the  United  States 

Henry  F.  Diaz  and  Robert  G.  Quayle 

National  Climatic  Center 

Asheville,  N.  C.  28801 

The  record-setting  winter  of  1977-78  was  unusual  in  a  number  of  ways.   On 
the  temperature  side,  record  cold  weather  prevailed  in  the  Nation's 
midsection  from  the  Ohio  valley  westward  to  the  Central  Plains.   It  followed 
on  the  heels  of  the  severe  winter  of  1976-77  which  had  set  many  records  for 
prolonged  cold,  and  had  produced  one  of  the  coldest  months  in  the  Nation's 
history  (January).   The  occurrence  of  such  extreme  back-to-back  winter 
seasons  has  had  only  two  or  three  precedents  in  the  instrumental 
meteorological  record:   1903-04  &  1904-05;  1855-56  &  1856-57  and  to  a  lesser 
degree  of  severity  a  series  of  cold  winters  from  1884-85  to  1887-88. 

Several  low  monthly  and  seasonal  temperature  records  were  set  in  winter  of 
1977-78.   It  was  the  coldest  February  on  record  in  six  states  and  the  second 
or  third  coldest  in  an  additional  sixteen.   It  was  the  coldest 
December-February  period  in  twelve  states  and  the  second  or  third  coldest  in 
another  eight.   The  January-February  mean  was  the  lowest  on  record  in  19 
states  and  second  or  third  coldest  in  four  more  states. 

On  the  precipitation  side,  the  extreme  2-year  drought  in  the  far  West  was 
decisively  broken  with  accumulations  in  some  areas  exceeding  twice  the 
normal  seasonal  precipitation  amounts.   The  Northeast  also  experienced  heavy 
precipitation,  as  record  snowfall  exceeding  40  inches  in  parts  of  southern 
New  England  occurred  in  February.   January  1978  was  one  of  the  wettest  on 
record  in  the  Northeast. 

Population-weighted  heating  degree  days  are  commonly  used  as  an  index  of 
fuel  demand  for  heating  purposes.   Because  of  a  westward  shift  of  the  region 
of  largest  negative  anomalies,  the  1977-78  heating  degree  day  season  did  not 
exceed  the  1976-77  season,  finishing  a  close  second  in  an  81-year  period  of 
record. 
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The  Low  Index  Winter  1977-78  Vs.  The  High   Index  Winter  1975-76 
A  Comparison  of  Energetics 

Robert  R.  Dickson 
National  Weather  Service,  NOAA 


For  many  years  it  has  been  common  to  characterize  the  large- 
scale  atmospheric  circulation  in  terms  of  the  strength  of  the  wes- 
terlies at  middle  latitudes--!' .e.  the  zonal   index.     A  high  index 
situation  involves  fast  westerlies  at  middle  latitudes  and  the  en- 
hanced transport  of  maritime  air  across  the  continents.     A  low  index 
regime,  on  the  other  hand,  displays  weak  mid-latitude  westerlies, 
blocking  ridges  at  high  latitudes  and  increased  meridional   transport 
of  air  masses. 

In  the  brief  interval   from  winter  1975-76  to  that  of  1977-78 
the  atmosphere  has  fluctuated  from  a  well   developed  high  index  situa- 
tion to  a  pronounced  low  index  pattern  in  the  Western  Hemisphere 
(Figures  1-3).     Over  the  United  States  this  has  produced  a  change 
from  widespread  warmth  in  winter  1975-76  to  extreme  temperature  con- 
trasts—warm in  the  west,  cold  in  the  east—during  the  1977-78 
winter  (Fig.  4). 

The  Energy  Analysis  Program  of  the  National  Meteorological 
Center  (Miller  et  al ,  1975)  provides  data  by  which  the  average 
energetics  behavior  of  these  two  "opposite"  types  of  winter  can  be 
viewed.     The  time  variation  of  available  energetic  parameters  pre- 
ceding, during  and  following  the  two  winters  is  shown  in  Figure  5. 
These  quantities  are  integrated  over  the  Northern  Hemisphere  from 
1000  to  100  mb,  and  averaged  over  monthly  and  mid-month  to  mid- 
month  periods.     It  should  be  noted  at  the  outset  that  the  zonal 
kenetic  energy  (Kz)  of  the  "low  index"  winter,  1977-78,  exceeded 
that  of  the  "high  index"  winter  1975-76--by  as  much  as  25%  in 
February.     Beginning  in  fall   and  continuing  through  most  of  the 
winter  to  spring,  1977-78  was  a  greater  converter  of  zonal   to  eddy 
available  potential  energy  and  eddy  to  zonal   kinetic  energy 
(Cn  and  C|<)  and  Kz  was  higher.     During  winter  1977-78,  zonal  avail- 
able potential  energy  (Az)  remained  greater  despite  generally 
stronger  conversions  through  the  energy  cycle  to  sustain  the 
stronger  zonal   flow.     Factors  involved  in  this  replenishment  of  AE 
probably  included  radiative  effect  of  the  late  winter  snow  cover 
(Fig.  6),  as  well   as  the  insolation-  excluding  effects  of  the 
extensive  cloud  cover  which  must  have  occurred  at  mid-latitudes. 

The  latitudinal   distribution  of  vertically  integrated  values 
of  the  various  energy  parameters  and  the  northward  eddy  momentum 
transport  for  February  (Krueger  1978),  month  of  peak  Kz  in  the  two 
winters  under  study,  is  given  in  Figures  7  and  8.     The  excess  in 
Northern  Hemisphere  Kz  during  winter  1977-78  (Fig.  5)  is  now  seen 
to  be  highly  concentrated  in  a  latitudinal   sense  with  a  peak  value 
in  February  1978  more  than  twice  that  of  February  1976.     This  con- 
centration was  supported  by  increased  energy  conversions  south  of 
45°N  latitude   (Figs.  7C  and  D).     The  low  latitude  maximum  in  the 
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conversion  from  eddy  to  zonal   kinetic  energy  involves  an  augmented 
eddy  momentum  transport  out  of  the  tropics   (Fig.  8).     Eddy  kinetic 
energy  in  February  1978  (Fig.   7B)  was  concentrated  just  north  of 
the  Kz  max  where  kinetic  energy  was  fed  from  the  zonal   flow  to 
eddies   (negative  CK,  Fig.   7D).     In  contrast,  in  February  1976,  uti- 
lization of  available  potential  energy  was  displaced  northward  and 
momentum  transport  out  of  the  tropics  and  the  associated  contri- 
bution to  zonal   kinetic  energy  were  reduced. 

Thus  the  low  index  winter  is  seen  as  one  with  a  concentration 
of  zonal   kinetic  energy  at  relatively  low  latitudes  due  to  the 
cooperative  effect  of  a  southward  displaced  mid-latitude  baroclinic 
zone  and  accompanying  energy  conversions  together  with  the  enhanced 
export  of  kinetic  energy  from  the  tropics.     The  two  events,  futher- 
more,  are  not  independent.     An  initial  substantial   displacement  of 
the  winter  mid-latitude    westerlies  south  of  normal   at  appropriate 
longitudes  brings  westerly  troughs  into  the  tropics.     As  pointed 
out  by  Riehl    (1969)  such  troughs  in  the  upper  troposphere  are  a 
significant  means  for  momentum  transport  out  of  the  tropics  with 
maximum  transport  possible  at  those  longitudes  where, jointly,  the 
subtropical  jet  stream  reaches  fartherest  north  and  the  polar 
troughs  intrude  fartherest  south--i.e.  near  the  east  coasts  of 
Asia  and  North  America  and,  to  a  lesser  extent,  over  the  Near 
East.     Reference  to  the  mean  700  mb  height  maps  for  the  two  winters 
(Figs.  2  and  3)  reveals  an  increased  opportunity  for  such  inter- 
actions in  the  low  index  winter  of  1977-78  and  a  diminished 
opportunity,  in  the  high  index  winter  of  1975-76. 
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LEGEND  TO  FIGURES 


Fig.  1       700  mb  mean  geostrophic  wind  speed  profiles  for  the 
Western  Hemisphere  during  the  winters  of  1975-76  and 
1977-78. 

Fig.  2      Mean  700  mb  height  (dam)  and  DN  (m)  for  winter  1975-76. 

Fig.   3      Mean  700  mb  height  (dam)  and  DN   (m)  for  winter  1977-78. 

Fig.  4      Mean  surface  air  temperature  DN  categories  for  winters 
of  1975-76  and  1977-78. 

Fig.  5      Time  variation  of  monthly   (and  mid-month  to  mid-month) 
mean  values  of  energetics  parameters  for  the  Northern 
Hemisphere,  1000  to  100  mb,  for  portions  of  1975-76 
and  1977-78. 

Fig.  6      Snow  and  ice  boundaries  for  February  9-15,  1976  and 
February  13-19,  1978. 

Fig.   7      Latitudinal  profiles  of  February  mean  values  of 

vertically  integrated  energetics  parameters,  1000  to 
50  mb,  for  1976  and  1978. 


Fig.  8      Latitudinal   profiles  of  vertically  integrated  north- 
ward transport  of  eddy  angular  momentum,  1000  to  50  mb, 
for  February  1976  and  1978. 
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Fi^.  6      Snow  and  ice  boundaries  for  February  9-15,  1976  and 
February  13-19,   1978. 
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Fig.  8      Latitudinal  profiles  of  vertically  integrated  north- 
ward transport  of  eddy  angular  momentum,  1000  to  50  mb, 
for  February  1976  and  1978. 
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Fig.  8      Latitudinal  profiles  of  vertically  integrated  north- 
ward transport  of  eddy  angular  momentum,  1000  to  50  mb, 
for  February  1976  and  1978. 
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A  HEMISPHERIC  VIEW  OF  THE  TROPOSPHERE  AND  STRATOSPHERE 
DURING  THE  WINTERS  1975-76,  1976-77,  AND  1977-78 

R.  S.  Quiroz 

National  Meteorological  Center,  NOAA,  Washington,  D.C.  20233 

ABSTRACT 

Zonal  averages  of  temperature,  height,  and  mean  zonal  wind  are 
portrayed  for  selected  altitudes  and  latitudes  of  the  stratosphere  and 
troposphere  for  the  last  3  winters.  From  these  data  the  following 
findings  are  presented:  (1)  Coupling  in  the  stratospheric  thermal 
fields  at  low  and  high  latitudes  was  strongest  in  1976-77,  with  a 
negative  correlation  of  0.92  between  low-  and  high-latitude  oscillations 
(2)  Coupling  in  the  tropospheric  thermal  fields  in  middle  and  high 
latitudes,  as  indicated  by  1 00-mb  height  data,  was  also  strongest  in 
1976-77.   (3)  Tropospheric  jet  stream  and  stratospheric  polar  night 
jet  strengths  were  inversely  related  in  the  3  winters,  with  a  delay  of 
maximum  winds  in  the  troposphere  (minimum  winds  in  the  stratosphere) 
from  January  in  1976-77  to  February  in  1977-78.   (4)  In  1976-77  and 
1977-78,  strong  high-latitude  stratospheric  warmings  were  succeeded  by 
polar  tropospheric  warming  and  mid-latitude  tropospheric  cooling,  with  a 
delay  of  about  one  month  in  these  occurrences  in  winter  1977-78  relative 
to  1976-77.  The  negative  temperature  anomaly  defined  by  40°N  zonal  mean 
temperatures  (compared  to  normal)  and  integrated  over  the  periods  of 
significant  cooling  in  December-February  is  about  two-thirds  as  great  in 
1977-78  as  in  1976-77.  No  physical  interpretation  of  these  results  is 
attempted  here;  it  is  meanwhile  evident  that  meriodional  heat  transports 
were  important  in  both  the  troposphere  and  stratosphere. 
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"19  78  Break  in  the  West  Coast  Drought" 

Jerome  N ami as 

Scripps  Institution  of  Oceanography 
La  Jolla,  California  92093 


This  paper  describes  how  two  years  of  back-to-back  drought 
in  California  and  the  Far  West,  during  the  winters  of  19  76  and 
19  77,  abruptly  changed  to  an  anomalously  wet  winter  in  19  78. 
After  first  describing  some  of  the  conditions  of  the  general 
circulation  and  the  sea-surface  temperatures  preceding  and  during 
the  19  77  drought,  the  differing  patterns  antecedent  and  during  the 
winter  of  1978  are  detailed. 

These  differences  were  largely  associated  with  a  strong  low 
latitude  westerly  flow  pattern  in  1978  in  contrast  to  a  strong 
meridional  ridge  in  the  preceding  year,  and  also  a  different  SST 
anomaly  distribution.   These  striking  differences  began  to  show  up 
as  early  as  September  19  77  over  much  of  the  Northern  Hemisphere 
as  well  as  over  the  North  Pacific.   Specifications  and  predictions 
of  upper  level  heights  from  sea-surface  temperatures  strongly  in- 
dicate the  break  in  the  drought  regime,  and  many  charts  are  shown 
objectively  relating  precipitation  to  pressure  patterns  and  to 
teleconnections . 

Some  physical  explanations  of  the  above  material  are  offered. 
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"1978  BREAK  IN  THE  WEST  COAST  DROUGHT" 

CLIMATE  DIAGNOSTICS  WORKSHOP  -  OCTOBER  31  -  NOVEMBER  2,  1978 

MIAMI,  FLORIDA 
JEROME  NAMIAS 

LEGENDS  TO  FIGURES 




Fig.  1     Percentage  of  normal  precipitation  over  the  West  for  the  past 
three  winters . 

Fig.  2     700  mb  height  and  departures  (Broken  Lines)  from  the  long  term 
mean  for  Winter  1977. 

Fig.  3     Sea  level  pressure  and  anomalies  (Broken)  in  mb  for  Winter  1977. 

Fig.  4     700  mb  mean  height  patterns  for  the  Winters  of  1977  and  1978. 

Fig.  5     700  mb  height  anomalies  and  sea-surface  temperature  (SST) 

anomalies  (Shaded)  for  the  Winters  of  1977  and  1978.   Slant 
shading  indicates  Below  Normal  Temperatures,  Stippling — Above 
Normal . 

Fig.  6     700  mb  height  and  standardized  anomalies  for  the  Winter  1977-78. 

Fig.  7     Mean  sea-level  pressure  for  January,  1978. 

Fig.  8     Correlations  between  precipitation  in  the  North , Central  and 

Southern  California  Water  Sheds  (Boxes)  and  the  mean  sea-level 
pressure  elsewhere.   Correlations  are  for  26  years  of  data. 

Fig.  9     Sea-surface  temperature  anomalies  for  the  Septembers  of  1976 
and  1977. 

Fig.  10    Sea-surface  temperature  anomalies  for  the  Falls  of  1976  and  1977. 

Fig.  11    700  mb  height  and  standardized  departures  from  normal  for  Sept., 

1977. 
Fig.  12    700  mb  height  standardized  anomalies  for  the  Falls  of  1976  &  1977 

Fig.  13    Winter  1978  700  mb  anomalies  Specified  from  the  contemporary 
Sea-surface  temperature  anomalies. 

Fig.  14    Predicted  700  mb  anomalies  for  the  1977-78  Winter  made  at  the 
end  of  the  preceding  Sept.,  Oct.,  and  November. 

Fig.  15    700  mb  height  and  anomalies  for  the  Winter  of  1977-78  predicted 
assuming  that  the  standardized  anomalies  of  Fall  1977  are  con- 
served. 
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LEGENDS  TO  FIGURES 
J.  Namias 

Fig.  16   Winter  teleconnec tions  (cross  correlations)  of  700  mb  anomalies 
from  a  center  at  40°N-140°W. 

Fig.  17   Winter  teleconnections  of  700  mb  anomalies  from  a  center  at 
65°N-115°W. 

Fig.  18   Winter  teleconnections  of  700  mb  anomalies  from  a  center  at 
45°N-125°W. 
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Large-scale  Atmospheric  Circulation  and  the  West  Coast  Upwelling  Regime 

in  1976/77  and  1977/78 

Gunter  R.  Seckel,  Andrew  Bakun  and  Richard  Parrish 

Pacific  Environmental  Group,  SWFC 

National  Marine  Fisheries  Service,  NOAA 

Monterey,  California  93940 


Anomalies  in  the  large-scale  atmospheric  circulation  affect  the 
upwelling  regime  off  the  west  coast  of  North  America.  In  turn,  these 
anomalies  affect  the  stocks  of  valuable  fishery  resources  in  the  U.S. 
200  mile  fisheries  management  zone. 


THE  UPWELLING  REGIME 

During  spring  and  summer  the  North  Pacific  High  and  the  Continental 
Low  pressure  systems  produce  prevailing  northwesterlies  over  the  off- 
shore region  of  North  America.  During  fall  and  winter  the  Aleutian 
Low  and  a  pressure  ridge  extending  from  the  continent  over  central 
California  and  westward  over  the  eastern  Pacific  produce  prevailing 
southwesterlies  off  Oregon  and  Washington  and  northwesterlies  off 
southern  California.  The  wind  stress  climatology  (Fig.  1,  2,  3,  4) 
(Nelson,  1977)  shows  that  maximum  northwesterlies  begin  to  strengthen 
off  southern  California  early  in  the  year.  In  April  maximum  north- 
westerlies occur  off  central  California  and  in  summer  they  are  found 
off  northern  California.  During  summer  the  entire  offshore  region 
from  Canada  to  Baja  California  is  under  the  influence  of  the  pre- 
vailing northwesterlies.  From  October  to  April  southwesterlies 
prevail  off  Oregon  and  Washington. 

Northwesterlies  off  the  west  coast  produce  an  offshore  component 
of  surface  flow  (Ekman  transport) .  This  water  is  replaced  in 
part  with  colder,  nutrient  rich  "upwelled"  water  from  below  the 
surface  layer.  Southwesterlies  produce  an  onshore  component  of 
surface  flow,  "piling  up"  warm  surface  water  against  the  coast. 
This  process  can  be  termed  "downwelling." 

The  upwelling  regime  also  affects  the  along  shore  California  and 
Davidson  currents.  Upwelling  during  spring  and  summer  changes 
the  baroclinicity  in  such  a  manner  as  to  increase  the  southward 
flow  of  the  California  Current.  Onshore  surface  flow  and  down- 
welling  would  tend  to  increase  the  northward  flowing  Davidson 
current . 
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The  offshore  or  onshore  Ekman  transport  has  been  computed  from  sea- 
level  atmospheric  pressure  distributions.  These  transports  are  used 
as  coastal  divergence  or  upwelling  indices  (Bakun,  1973,  1975).  Daily 
upwelling  indices  based  on  6-hourly  pressure  analyses  are  computed 
at  three-degree  latitude  intervals  along  the  west  coast  from  Baja 
California  to  the  Gulf  of  Alaska.  Examples  for  1975  are  shown  in 
Figure  5  and  illustrate  the  upwelling  regime  associated  with  the 
prevailing  wind  systems  described  above. 


FALL  AND  WINTER  1976/77  and  1977/78 

The  winter  atmospheric  circulation  during  1976/77  and  1977/78 
(Fig.  6,  7)  produced  upwelling  regimes  that  are  distinct  from  the 
more  usual  pattern  that  occurred  in  1975.  The  anomolous  patterns 
are  most  distinct  off  central  California,  lat.  36°N  (Fig.  8). 

Anomalously  low  wind  speeds  and,  therefore,  low  index  values 
prevailed  during  the  fall  of  1976  through  February  1977.  In 
October  and  November  1977  stronger  than  normal  upwelling  persisted. 
However,  from  December  1977  through  the  winter  of  1978,  the  pattern 
for  Central  California  resembled  that  normally  found  off  Oregon 
and  Washington  with  recurring  strong  downwelling  events.  The  strong 
downwelling  events  also  occurred  off  southern  California  at  lat.  33°N. 


FISHERIES  IMPLICATIONS 

The  biological  significance  of  upwelling  in  terms  of  sea  surface 
temperature  and  biological  productivity  is  well  known.  Other  aspects 
of  the  upwelling  system  may  be  particularly  important  during  the 
egg  and  larvae  stages  of  fishery  resources  (Parrish,  Bakun  and  Nelson, 
manuscript) . 

An  offshore  component  of  surface  flow  would  tend  to  disperse  eggs 
and  larvae  into  areas  remote  from  the  favorable  coastal  environment. 
Onshore  surface  flow  would  act  to  hold  them  in  the  near-coastal 
environment  in  which  they  were  spawned.  We  observe  striking  con- 
formities between  fish  reproduction  strategies  and  the  seasonality 
and  geography  of  onshore -offshore  transport.  Off  Oregon  and 
Washington,  important  species  of  fish  spawn  during  the  winters 
when  onshore  wind  driven  flow  prevails  rather  than  during  the  more 
productive  upwelling  season.  In  the  southern  California  Bight, 
where  a  closed  gyral  circulation  provides  favorable  drift  conditions, 
these  types  of  species  are  much  reduced  and  the  fauna  is  dominated 
by  spring  spawners.  In  the  region  of  vigorous  upwelling  off 
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Northern  California  which  is  characterized  by  strong  offshore  surface 
transport  through  most  of  the  year  there  is  a  paucity  of  locally 
spawning  species.  Rather,  the  fish  stocks  which  harvest  the  massive 
productivity  of  that  region  are  mainly  migrating  species  which 
spawn  in  the  Southern  California  Bight. 

During  the  slack  wind  conditions  from  November  1976  to  February 
1977  organisms  were  subjected  to  much  lower  than  normal  levels  of 
onshore  and  offshore  excursions  of  surface  waters  and  to  turbulent 
mixing  and  related  dispersion  of  planktonic  aggregations.  The 
dominant  winter  spawners  of  the  central  California  coast  are  the 
rockfishes.  Several  species  apparently  had  excellent  spawning  success 
as  unusually  large  numbers  of  juvenile  rockfishes  were  reported  by 
personnel  of  the  California  Department  of  Fish  and  Game.  In  marked 
contrast,  the  recurring  strong  storm  events  during  1977/78  with  on- 
shore-offshore excursions  of  wind-driven  surface  flow  would  have  made 
it  difficult  for  planktonic  larvae  to  remain  in  favorable  bathy- 
metric  ranges.  Strong  turbulent  mixing  would  have  tended  to  destroy 
patch  structure  in  the  plankton.  Anomalous  coastal  convergence  or 
piling  of  surface  waters  against  the  coast  throughout  the  California 
region  would  tend  to  induce  northward  alongshore  flow  and  anomalously 
warm  advection.  The  onshore  Ekman  transport  events  at  30°N  are  without 
precedent  in  the  1967-1978  data  base.  This  anomalous  onshore  transport 
during  the  northern  anchovy's  spawning  season  may  have  concentrated 
both  anchovy  larvae  and  forage  near  the  coast  with  resulting  good 
spawning  success  in  relation  to  spawning  stock  size.  The  large 
numbers  of  young -of -the -year  fish  in  southern  California  waters  as 
well  as  the  unusually  large  concentrations  of  young  fish  as  far  north 
as  Monterey  Bay  during  the  summer  and  fall  of  1978  appear  to  support 
this  hypothesis. 
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Figure  1.  Resultant  surface  wind  stress  vectors  for  January.  Stresses  off 
Oregon  and  Washington  reflect  the  prevailing  southwesterly  winds 
to  be  expected  from  the  atmospheric  pressure  distribution.  South 
of  Point  Conception  stress  vectors  reflect  relatively  weak  north- 
westerlies,  with  maximum  stress  occurring  off  Baja  California. 
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Figure  2.  Resultant  surface  wind  stress  vectors  for  April.  The  prevailing 
winter  stress  pattern  off  Oregon  and  Washington  is  beginning  to 
weaken  while  northwester lies  off  California  continue  to  strengthen 
with  highest  stress  occurring  off  central  California. 
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Figure  3.  Resultant  surface  wind  stress  vectors  for  July.  The  entire  west 
coast  is  dominated  by  northwesterlies.  The  area  of  high  stresses 
extends  from  Point  Conception  to  Oregon  with  highest  stresses 
occurring  off  northern  California  near  Cape  Mendocino. 
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Figure  4.  Resultant  surface  wind  stress  vectors  for  October.  Off  Oregon 
and  Washington  the  winter  stress  distribution  pattern  is  be- 
ginning to  set  in.  Off  California  stress  magnitudes  continue 
the  decrease  that  began  in  September. 
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Tropical  Circulation  and  Radiation  Fluctuations: 
interactions  with  Mid-Latitude  Circulations  During  1976-1978 

Jay  S.  Winston  and  Arthur  F.  Krueger 
NOAA-NESS,  Washington,  D.  C. 


1 .  Summary 

The  remarkable  displacement  of  westerlies  toward  the  equator 
during  the  past  two  Northern  Hemisphere  winters  was  associated  with 
a  weakening  of  the  large  scale  standing  eddies  over  the  tropics  and 
an  eastward  displacement  of  major  convective  cloudiness  and  rainfall 
in  the  tropical  Pacific.  These  events  in  turn  appear  to  have 
contributed  to  a  stronger  Hadley  circulation.  The  changes  in  the 
tropics  were  first  noted  during  the  1976  southern  winter  when  the 
South  Pacific  anticyclone  and  tradewinds  weakened  and  equatorial 
sea  surface  temperatures  increased.   In  much  of  the  tropical  central 
Pacific  the  cloudiness  regime  was  even  more  intense  in  the  1977-78 
(northern)  winter  than  in  1976-77. 

2.  Relevant  zona  I  I y  averaged  features  of  the  circulation 

Zonal  kinetic  energy  distributions  for  the  last  three  Januaries 
are  shown  in  Fig.  I  (top).  Note  in  particular  the  very  high  kinetic 
energy  between  25N  and  35N  during  these  past  two  winters  in  contrast 
to  1976.  This  was  accompanied  by  a  larger  meridional  temperature 
gradient  in  these  latitudes  (Fig.  2).  This  expansion  of  the  wester- 
lies and  temperature  gradient  equatorward  was  accompanied  by  a  similar 
shift  of  the  subtropical  anticyclones  and  also  the  baroclinic 
instability.  The  baroclinic  activity  can  be  monitored  with  the 
conversion  of  zonal  to  eddy  available  potential  energy  (C^)  which 
is  shown  in  the  upper  part  of  Fig.  3.   Note  the  larger  values  in  1977 
and  1978  extending  from  40N  to  about  1 5N  in  contrast  to  1976. 
Larger  conversions  from  eddy  to  zonal  kinetic  energy  (CK)  also  occurred 
these  past  two  winters  in  the  subtropics  (lower  part  of  Fig.  3). 
While  the  westerlies  were  receiving  larger  amounts  of  kinetic  energy 
on  their  southern  flanks  they  were  losing  energy  on  their  northern 
flanks.   Note  in  particular  the  larger  negative  values  of  CK  near 
40N  in  association  with  large  positive  values  of  C^,  indicating  that 
the  eddies  in  middle  latitudes  grew  in  response  to  both  baroclinic 
and  barotropic  processes.  The  stronger  eddy  processes  occurring  at 
lower  latitudes  these  past  two  winters  may  have  played  a  role  in 
reducing  the  amplitude  of  the  large  standing  eddies  over  the  tropics. 
This  loss  in  tropical  eddy  kinetic  energy  is  shown  in  the  lower  part 
of  Fig.  I.   To  a  large  extent  this  reflects  a  loss  of  kinetic  energy 
of  the  lowest  two  wave  numbers  as  shown  in  Fig.  4  and  represents  a 
weakening  of  the  Walker  circulation. 

These  events  occurring  in  the  wintertime  circulations  appear  to 
signify  a  stronger  Hadley  circulation  occurring  during  these  last  two 
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years.  There  is  indirect  evidence  of  this  in  the  outgoing  longwave 
radiation  shown  on  the  left  side  of  Fig.  5.  Note  in  particular  the 
low  values  between  5S  and  5N  during  the  last  two  Januaries  in 
contrast  to  January  1976,   Since  this  parameter  is  largely  modulated 
by  major  cloudiness  over  the  tropics  this  suggests  increased  cloudi- 
ness, rainfall,  and  upward  motion  during  the  past  two  winters.   In 
contrast,  at  I 5N  increased  values  in  longwave  cooling  occurred  the 
last  two  winters  suggesting  decreased  cloudiness  and  greater 
subsidence.   Also  note  the  generally  higher  temperatures  over  the 
tropics  and  subtropics  in  contrast  to  1976  (right  side  of  Fig.  5). 
This  increased  temperature  (or  total  potential  energy)  very  likely 
results  from  the  increased  release  of  heat  of  condensation  suggested 
by  the  decreased  longwave  radiation. 

The  eddy  momentum  flux  has  also  been  used  to  infer  the  strength 
of  the  mean  meridional  circulation.   In  particular  an  approximate 
balance  between  the  Coriolis  term  and  horizontal  eddy  flux  has  been 
noted  in  the  upper  troposphere  by  Newell  and  his  colleagues.  Assuming 
this  balance  we  have  made  estimates  of  200  mb  meridional  flow  (Fig.  6), 
Note  the  three  cell  structure  and  also  the  stronger  northward  motion 
at  low  latitudes  these  past  two  winters. 

3.   Regional  fluctuations  in  longwave  radiation  (cloudiness)  and 
circulation  over  the  tropics 

Examination  of  patterns  of  radiation  and  circulation  over  the 
past  few  years  has  revealed  that  the  major  fluctuations  in  the  tropics 
have  occurred  over  the  Pacific  region.   By  the  southern  winter  of 
1976  the  South  Pacific  anticyclone  and  tradewinds  had  weakened 
substantially.  This  is  shown  by  several  of  the  time  series  in 
Fig.  7.  Note  the  Easter-Darwin  pressure  anomaly  difference  which  is 
a  simple  measure  of  the  Southern  Oscillation.  This  peaked  during 
1975,  then  declined  slowly,  dropping  below  normal  during  1976  where 
it  has  essentially  remained  except  for  several  short  upward  excur- 
sions. Variations  in  the  strength  of  the  tradewinds  for  two  points 
along  the  equator  are  shown  just  below  with  the  scale  reversed  - 
i.e.  easterlies  above  the  axis.  Note  the  strong  trades  during  1975 
and  1976.  At  I 70E  they  declined  to  near  zero  during  1976  and  even 
became  westerly  during  this  past  northern  winter.  Since  strong 
southern  tradewinds  accompany  a  well-defined  equatorial  dry  zone 
extending  well  west  of  180  we  can  also  use  satellite  observations  to 
monitor  this.  Two  such  time  series  in  longwave  cooling  are  shown  at 
the  top  of  Fig.  7.  At  the  equator  and  180  the  decline  in  cooling 
(increased  major  cloudiness)  set  in  by  June  1976.   South  of  the 
equator,  towards  French  Polynesia  (I0S,  I60W),  it  set  in  late  in 
1976. 

The  spatial  variation  of  this  outgoing  longwave  radiation  over 
the  Pacific  for  the  past  three  winters  is  shown  in  Fig.  8.   Note 
the  eastward  spread  of  lower  values  (e.g.,  below  250  W  m~2)  between 
ION  and  I0S  from  (northern)  winter  of  1976  (Dec.  1975-Feb.  1976)  to 
winter  of  1977  and  the  further  eastward  spread  by  winter  of  1978. 
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These  lower  values  signify  substantial  increases  in  major  cloudiness 
through  the  heart  of  the  central  tropical  Pacific.   Note  also  the 
accompanying  major  eastward  retreat  of  the  equatorial  dry  zone. 
Taking  the  280  W  m   contour  as  an  indicator  of  the  extent  of  a 
nearly  unadulterated  dry  zone  we  see  that  its  western  edge  retreated 
from  near  I75W  in  winter  1976  to  I60W  in  1977,  and  then  to  near 
1 45W  in  1978.   Fig,  9  illustrates  the  patterns  of  year-to-year 
differences  in  these  longwave  radiation  values,  The  eastward  shift 
of  the  large  area  of  decreases  is  striking,  Note  that  much  of  the 
central  Pacific  region  (roughly  from  I 50W  to  I65E)  experienced  falls 
in  longwave  cooling  for  two  consecutive  years.  Accompanying  the 
eastward  march  of  the  falls  is  a  large  area  of  rises  over  the  western 
Pacific  in  1978-77  signifying  substantial  decreases  in  the  major 
cloudiness  over  the  normally  convectively  active  "maritime  continent". 

The  200  mb  circulation  anomalies  for  the  past  three  winters  are 
shown  in  Fig.  10.*   In  many  respects  the  last  two  are  similar  and 
quite  different  from  that  for  1975-76.  Note  the  relative  easterly 
flow  over  the  equatorial  Pacific,  with  anticyclonic  anomalous  flow 
both  north  and  south  in  contrast  to  the  equatorial  westerly  flow 
and  cyclonic  anomalous  flow  north  and  south  during  the  1975-76  winter. 
Interestingly,  westerly  anomalous  flow  appeared  over  the  equatorial 
Atlantic  the  past  two  winters.   Another  feature  worth  noting  is  the 
anomalous  flow  from  the  eastern  tropical  Pacific  toward  the  north  in 
all  three  years.   In  the  past  two  winters  this  flow  penetrated  into 
southern  and  southeastern  United  States,  suggesting  frequent 
incursions  of  deep  tropical  air  into  the  subtropical  edge  of  the 
southward-displaced  westerlies  over  the  United  States  and  the  Atlantic 

The  type  of  shorter-period  activity  that  contributed  to  the 
anomalous  seasonal  mean  patterns  of  radiation  and  circulation  of 
the  past  few  winters  is  illustrated  in  Figs.  II  and  12.  The  GOES 
infrared  imagery  for  February  19,  1978  (Fig.  II)  shows  large  clusters 
of  convective  activity  south  of  the  equator  in  the  central  Pacific. 
(Three  tropical  cyclones  were  located  in  the  southern  intertropical 
convergence  zone  just  north  of  Polynesia.)  High  level  outflow  from 
this  region  could  be  seen  in  the  GOES  imagery  to  extend  across  the 
equator  to  join  up  with  the  Northern  Hemisphere's  subtropical  jet 
stream.  Note  the  long  cirrostratus  -  altostratus  band  extending 
from  the  South  Paci f i c, across  Mexico  out  over  the  Atlantic,  a  distance 
of  12,000  km!   Six  day  mean  patterns  of  circulation  and  longwave 
radiation  for  February  16-21,  1978  (Fig.  12)  illustrate  the  abnormal 
eastward  displacement  of  the  southern  ITCZ  and  the  associated  flow 
patterns.   The  west-east  axis  of  very  low  outgoing  longwave  radiation 
extending  from  about  I 60E  to  I 50W  is  located  generally  near  the 
east-west  1000  mb  trough  (indicated  by  the  H  labels)  extending  from 
Indonesia  to  French  Polynesia.  The  200  mb  flow  over  this  major 


*  The  anomalies  represent  departures  from  6  year  mean  seasonal  values 
based  on  data  from  1968-74.   H  refers  to  centers  of  clockwise  circu- 
lation, L  refers  to  centers  of  counterclockwise  circulation. 
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cloudiness  and  low  level  convergence  is  ant icyclon ic.  At  the  eastern 
end  of  this  active  cloud  region  the  200  mb  flow  is  most  strongly 
anticyclonic  and  there  is  a  concentrated  flow  from  the  ITCZ  north- 
ward across  the  equator.  North  of  the  equator  this  flow  divides, 
one  branch  continuing  mainly  eastward  into  the  subtropical  jet  stream 
over  Mexico,  while  the  other  branch  returns  to  the  south  circulating 
about  the  pronounced  low  latitude  trough  south  of  the  equator. 
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Figure    I 
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Figure  2 
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Figure  4 


200mb  Eddy  Kinetic  Energy  over  the  Tropics 
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Figure   5 
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Figure  6 
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Figure   7 
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OUTGOING  LONGWAVE  RADIATION  (WM~2) 
W I NTER  1 976 


OUTGOING  LONGWAVE  RADIATION  (  W  M  "  2  ) 
WINTER  1977 


OUTGOING  LONGWAVE  RADIATION  (WM'2) 
WINTER  1978 
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Figure  9 
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Figure    12 
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SNOW  AND  SEA  ICE  COVER  FLUCTUATIONS  IN  1977-78 

G.  Kukla  and  J.  Gavin 

Lamont-Doherty  Geological  Observatory 
Palisades,  New  York   10964 


Variation  of  snow  and  pack-ice  in  both  hemispheres  was 
measured  from  U.S.  Navy  and  NOAA  charts  separately  for  vari- 
ous geographic  segments  (Figs.  1  and  2) .  The  methodology 
used  is  described  in  detail  elsewhere  (Kukla  and  Gavin,  1978) . 

NOAA  charts  were  occasionally  cross-checked  and  comple- 
mented by  information  from  Air  Force  charts  and  DMSP  satellite 
images.  The  continued  improvement  in  the  quality  of  the  snow 
charts  through  1978  should  be  noted.  Data  checks  are  made  on 
a  continual  basis  by  exchanging  results  with  that  generated 
by  another  method  (Wiesnet  and  Matson,  1975). 

The  northern  hemisphere  snow  charts  were  incomplete  for 
the  spring  of  1978  due  to  instrument  failure.  Air  Force 
photographs  were  analysed  in  order  to  complete  the  set. The  charts 
have  been  made  available  to  NOAA  for  further  distribution. 

Snow  and  ice  coverage  during  1977  was  on  the  heavy  side 
in  the  northern  hemisphere,  but  light  in  the  southern  one.  In 
February  1978  the  largest  average  monthly  snow  and  ice  extent 
in  the  northern  hemisphere  observed  over  the  past  decade  was 
reached.  The  largest  12  month  running  average  was  reached  in 
the  interval  ending  August  1978.  (Figs.  3  and  4) . 
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Compared  to  the  19  73-77  season  ~   mean  ,  Arctic  ice  was 
above  normal  in  winter  and  spring  but  below  in  summer  and 
autumn.  Snow  was  below  normal  in  spring  and  summer  but  above 
average  in  autumn  and  winter.  (Fig. 10) 

Large  year  to  year  regional  Arctic  ice  variations  can 
be  observed.  (Fig.  7  ) .  Note  the  opposite  trends  between  the 
Bering  Sea  (BER)  and  the  North  Atlantic  sector  (NAT)  which 
includes  the  Barents,   Norwegian  and  Greenland  Seas. 


Kukla,  G.  and  Gavin,  J.  1978.  Proc.  of  the  Intern.  Workshop 

on  World  Glacier  Inventory,  17-22  Sept.,  1978, Switzerland. 

Wiesnet,  D.  R.  and  Matson,  M.  1975.  Monthly  Weather  Review, 104, 
323-835. 
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FIGURE  1:   Geographic  segments  analysed  by  grid  count  for  both 
the  northern  and  southern  hemisphere. 

FIGURE  2 :   Geographic  segments  for  the  northern  hemisphere 
(planimetered) 

FIGURE  3:   12  month  running  means  of  snow  and  ice  cover  in  the 
northern  hemisphere (NH-SIC) ;  of  sea  ice  in  concen- 
tration >4  octas  plus  the  Antarctic  continent  in 
the  southern  hemisphere  (SH-SIC) .  Plotted  on  terminal 
date.  Also  shown  Antarctic  sea  ice  and  land  ice  extent 
at  the  end  of  November  (SH-NOV) .  Scale  in  million  km2. 
Dots  mark  annual  average. 

FIGURE  4 :   Mean  monthly  snow  and  sea  ice  covers  in  the  Northern 

Hemisphere  from  1967-78.  Ice  ^4  octas,  large  polynyas 

2 

excluded,  Greenland  included.  Scale  in  million  km  . 

FIGURE  5:    Mean  monthly  snow  and  ice  cover  in  the  Northern  Hemis- 
phere from  1972-78  for  individual  segments. (Jan-Jun) 

FIGURE  6:    Same  as  Fig.  5   (Jul-Aug)  Scale  in  million  km2. 

FIGURE  7:    Mean  monthly  Arctic  ice  cover  from  1972-78  for  BER, 

LAB,  and  NAT.   Ice   >  4  octas.  Scale  in  0.1  million  km2 

FIGURE  8 :  Mean  monthly  sea  ice  extent  in  Antarctic  segments 
ANR  and  ANW  from  1973-78.  Minimum  concentration  4 
octas.  Units  in  millions  of  square  kilometers. 

FIGURE  9 :    12  month  running  means  of  Antarctic  sea  ice  for  5 
segments.  Plotted  at  terminal  date.  Upper  curves 
show  ice  over  1  octa  and  polynyas,  lower  curves  show 

sea  iceiM  octas  and  exclude  polynyas.  Lines  mark 

2 

beginning  of  the  year.  Scale  in  million  km  . 
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FIGURE  10:   Seasonal  difference  in  percent  from  the  1973-77  mean 
shown  for  Arctic  ice  and  northern  hemisphere  snow  and 
ice  combined.  Winter (Dec-Feb) ; Spring (Mar-May) ; Summer 
(June-Aug) ; Autumn (Sept-Nov) . 

FIGURE  11:  Difference  in  percent  from  the  1973-77  mean  for 
northern  hemisphere  snow  and  ice  for  individual 
segments.  Units  in  5%  increments. 
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The   12 -year  Satellite  Snowcover  Record  and 
the  Record  Winter  Snowcover  of  1977-78 

l  2 

Craig  P.   Bergx  and  Michael  Matson 

The  winter  of  1977-78  was  the  second  consecutive  year  of  record- 
breaking  winter  snowcover  in  North  America  according  to  NOAA/NESS 
satellite-derived  measurements.     February  1978,   in  particular,  had  the 
maximum  snowcover  on  record  for  the  12 -year  North  American  data  set  from 
1967-1978. 

The  data  set,  now  complete  from  November  1966  through  September  1978, 
consists  of  monthly  mean  snowcover  maps  and  their  corresponding  areal 
snowcover  values.     The  snowcover  data  have  been  corrected  to  correspond 
to  true  surface  area   (Matson,   1977) .     December  1967  and  January  1968 
imagery  from  the  ESSA  Series  satellites  was  previously  missing  making  the 
monthly  maps  for  these  two  months  unavailable   (Wiesnet  and  Matson,   1976) . 
Through  the  courtesy  of  Sea  Ice  Consultants,  new  images  were  obtained 
and  the  monthly  mean  snowcover  maps  and  snowcover  area  were  derived. 

The  updated  plot  of  the  12 -month  running  mean   (figure  1)   shows  an 
upward  trend  due  to  the  record  winter  snowcover  of  1977-78.     The  Eurasian 
running  mean  snowcover,   though     declining  in  1977,  now  shows  a  slight 
upward  trend.3  The  Northern  Hemispheric  running  mean  snowcover,   also 
declining  in  1977,    likewise  shows  an  upward  trend  as  a  result  of  the 
extensive  snowcover  during  the  winter  of  1977-78. 


•"■Lieutenant,  NOAA  Corps 
^Physical  Scientist,  NOAA/NESS 

The    Scanning    Radiometer    on    the    NOAA-5    satellite    failed    on    16 

MAR    78.     During    the    months    of    March    and    April    1978,     limited 
Defense    Meteorological    Satellite    Program    (DMSP)    data,     supplemented 
with    the    10-year    mean    monthly    snowline,    were    used    to    derive    the 

Eurasian    weekly    snowline.  ,~  , 


North  American  and  Eurasian  winter  satellite  snowcover  values  for 
the  twelve-year  period  are  graphically  illustrated  in  figure  2. 
February  1978  snowcover  in  North  America  was  the  highest  and  January 
1978  snowcover  the  second  highest  for  any  month  in  the  12 -year  snowcover 
record.  Eurasian  winter  snowcover  in  1977-78  remained  about  the  same 
as  the  1977  Eurasian  winter  snowcover. 

Figure  3  permits  a  comparison  of  the  snow  lines  of  the  record 
February  1978  Northern  Hemispheric  snowcover  with  the  10-year  (1967-76) 
mean  February  snowcover.  The  February  1978  snow  line  in  the  southern 
and  midwestern  regions  of  the  U.S.  is  generally  500  to  1000  kilometers 
south  of  the  10-year  satellite-derived  mean  snow  line;  whereas,  in  the 
West  the  snow  line  differed  insignificantly  from  its  10-year  mean. 
Eurasian  February  1978  snowcover,  when  compared  to  the  10-year  mean, 
shows  large  areas  of  more  extensive  snowcover  in  Eastern  Europe,  Central 
Asia,  Korea,  and  Japan. 

Table  1  contains  the  10-year  mean  monthly  winter  continental 
snowcover  and  the  1976-77  and  1977-78  winter  snowcovers.  January 
February,  and  March  1978  winter  snowcover  values  in  North  America  are 
greater  than  those  for  1976-77.  December  snowcover  in  both  winters 
was  within  the  expected  limits. 

Table  1  also  shows  that  December  and  January  of  both  winters  were 
unusual  snowcover  months  in  Eurasia,  whereas  February  and  March  of  both 
winters  were  within  the  expected  limits  of  snowcover  extent.  The  1977-78 
Eurasian  snowcover  was  more  extensive  in  December,  February,  and  March 
than  in  1976-77,  with  the  peak  February  snowcover  of  1978  being  less 
extensive  than  the  peak  January  snowcover  of  1977. 
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Figures  4-7  show  the  monthly  mean  snowcover  maps  for  the  Northern 
Hemisphere  for  the  winter  of  1977-78  for  comparison  to  the  10-year  mean 
monthly  snowcover  maps  for  the  Northern  Hemisphere.  Note  the  extensive 
snowcover  in  North  America  during  January,  February,  and  March. 

Several  statistical  relationships  between  snowcover  data  of  the 
winter  months  in  Eurasia  and  the  Northern  Hemisphere,  as  first  reported 
by  Wiesnet  and  Matson  (1976) ,  have  been  updated  to  include  eleven  data 
points.  Figures  8-10  are  the  11-point  updates  of  these  linear  regression 
relationships.  The  snowcover  estimates  for  Eurasia  and  Northern 
Hemisphere  shown  in  Table  2  were  derived  using  the  updated  regression 
equations  based  on  the  1977-78  snowcover  data.   In  figure  11,  we  compare 
these  estimates  with  the  measured  snowcover  values.  North  American 
snowcover  estimates  are  the  arithmetic  difference  of  Eurasian  from 
Northern  Hemispheric  snowcover  values,  i.e.,  e^  =  e^  -  eg,  where 
e  is  the  snowcover  estimate. 

References 
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Fig.  1.    12-month  running  mean  of  monthly  mean  snowcover  for  North 

America,  Eurasia,  and  the  Northern  Hemisphere  for  the  period 

1966-1978. 
Fig.  2.    Monthly  mean  winter  snowcover  of  North  America  and  Eurasia 

for  the  period  1966-1978. 
Fig.  3.    Comparison  of  the  satellite-derived  mean  snow  line  for 

February  1978  and  the  10-year  satellite-derived  mean  snow  line. 
Fig.  4.    Comparison  of  December  1977  mean  snow  line  and  10-year  mean  snow 

line,  1966-1975. 
Fig.  5.    Comparison  of  January  1978  mean  snow  line  and  10-year  snow 

line,  1967-1976. 
Fig.  6.    Comparison  of  February  1978  mean  snow  line  and  10-year  mean 

snow  line,  1967-1976. 
Fig.  7.    Comparison  of  March  1978  mean  snow  line  and  10-year  mean  snow 

line,  1967-1976. 
Fig.  8.    Regression  analysis  of  January  snowcover  vs.  December  snowcover 

for  Eurasia  and  the  Northern  Hemisphere. 
Fig.  9.    Regression  analysis  of  February  snowcover  vs.  December  plus 

January  snowcover  for  Eurasia  and  the  Northern  Hemisphere. 
Fig.  10.   Regression  analysis  of  January  through  March  snowcover  vs. 

December  snowcover  for  Eurasia  and  the  Northern  Hemisphere. 
Fig.  11.   Winter  1978  snowcover  estimates  compared  to  winter  1978 

snowcover  measurements. 
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THE  12-YEAR  SATELLITE  SNOWCOVER  RECORD 
AND  THE  RECORD  WINTER  SNOWCOVER  OF  1977-78 

Craig  P.  Berg  (LT,  NOAA  Corps)  and  Michael  Matson 

ABSTRACT 
The  winter  1977-78  satellite-derived  snowcover  data  were  examined 
and  compared  to  the  12-year  records  to  determine  possible  trends  and 
relationships.  Winter  1977-78  was  the  second  consecutive  year  of  record- 
breaking  continental  winter  snowcover  in  North  America  as  reflected  in 
graphical  analysis  of  both  winter  snowcover  and  the  upward  trend  in  the 
12-month  running  mean.  The  regression  analysis  equations  have  been 
updated  and  now  include  eleven  data  points;  correlation  coefficients 
remain  significant  for  possible  continental  snowcover  estimates. 
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Fig.    4 
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Fig.    5 
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Fig.   9 
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Table  2.   Winter  Snow  Cover  Estimate  Data  in  the  Northern  Hemisphere 
and  Eurasia  from  Regression  Analysis 


Snowcover 
Period 


Estimate 
(xl06km2) 


JAN  78 

45.8 

FEB  78 

47.0 

JAN-MAR  78 

132.1 

Measured 
(xl06km2) 


Percent 
Difference 


NORTHERN  HEMISPHERE 

46.4  -1.3 

47.3  -0.6 

134.8  -2.0 


Standard  Error 
of  Estimate 
(xl06km2) 


1.58 
1.60 
2.80 


JAN  78 

30.6 

FEB  78 

31.2 

JAN -MAR 

78 

87.7 

EURASIA 

28.9 

+5.9 

29.7 

+5.1 

84.1 

+4.3 

1.16 
1.51 
3.09 
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COMPARISON  OF  WINTER  1978  SNOW  COVER    ESTIMATES 
VERSUS  WINTER  1978  SNOW  COVER    MEASUREMENTS 
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The  Variation  of  Sea  Surface  Temperature  in  1976  and  1977 

The  Data  Analysis 

K.  Miyakoda  and  A.  Rosati 

The  objective  of  this  paper  is  to  investigate  the  variation  of  SST  over  the 
world's  oceans  for  the  years  of  1976  and  1977.  In  order  to  study  the  spacial 
distribution  of  the  SST,  ship  and  satellite  data  at  1°  quadrangles  were 
collected.  Two  points  were  investigated;  (i)  the  difference  of  monthly  mean 
SST  data  between  the  two  sources;  and  (ii)  map  analysis  over  the  world  ocean. 
The  SST  charts  were  then  compared  with  independent  analyses. 

Using  monthly  mean  SST  data  derived  from  NOAA-5  soundings,  the  quality 
of  the  data  was  examined  relative  to  the  observations  of  ships  of  opportunity. 
Table  1  shows  the  month-to-month  variations  for  two  regions  in  the  North 
Pacific.  Concerning  the  difference  of  the  spatial  mean  SST,  negative  values 
are  dominating  at  the  high  latitudes,  40°-60°N,  implying  that  the  satellite 
data  are  cooler  than  ship  data,  probably  due  to  cloud  coverage.  On  the  other 
hand,  in  the  lower  latitudes,  20°-40°N,  the  differences  are  more  positive. 
The  standard  deviation  ranges  from  +  1.18°  to  +  2.3°C,  +  1.49°C  being  the 
arithmetic  average  for  the  entire  period  and  domain.  Any  systematic  seasonal 
variation  in  the  satellite  data  deviation  was  hardly  detectable. 

Figure  1  shows  the  statistics  over  the  domains  of  the  North  Pacific 
(20°-60°N,  120°E-120°W)  and  the  North  Atlantic  (20°-60°N,  70°W-10°W)  for  the 
annual  averages  of  1976  and  1977  separately.  The  numbers  are  the  mean 
differences  and  the  standard  deviation,  and  sample  numbers  are  shown  in 
parenthesis.  In  April  1977,  a  new  method  to  detect  clouds  and  improve 
atmospheric  attenuation  correction  was  introduced  into  the  satellite 
retrievals.  This  resulted  in  a  reduction  of  the  sample  size  but,  supposedly 
improved  the  quality  of  the  data.  The  mean  difference  ranges  from  -1.3°C 
to  +.8°C,  and  the  standard  deviation  ranges  from  +  1.43°C,  which  is  not 
small  in  view  of  +  2.3°C  as  the  root-mean-square  of  the  monthly  mean  SST 
anomaly. 

We  understand  that  the  quality  of  the  satellite  SST  and  the  data 
processing  techniques  are  still  under  continuing  development.  In  years  to 
come,  we  have  to  rely  on  satellite  data  over  a  wide  area  of  the  ocean.  Yet 
there  seems  to  exist  a  fundamental  difference  between  the  satellite  and  the 
ship  observations;  one  is  a  skin  temperature  while  the  other  is  the  temperature 
at  a  depth  3^6  meters.  It  is  essential  to  develop  an  improved  technique  by 
which  one  can  only  extract  the  benefits  from  satellite  data,  and  assimilate 
or  amalgamate  it  properly  with  ship  data. 

Adopting  an  analysis  program  of  Levitus  and  Oort  and  using  ship  and 
satellite  data,  we  produced  monthly  mean  SST  charts  of  l°xl°  grid  over  the 
world  ocean  for  1976  and  1977.  Figure  2  shows  an  amalgamation  of  satellite 
and  ship  analysis  for  Jan.  1977. 

The  charts  agreed  well  with  the  independent  analysis  at  Fleet  Numerical 
Weather  Central,  Monterey,  so  far  as  the  large  scale  features  are  concerned. 
However,  GFDL's  analysis  contained  a  great  deal  of  small-scale  irregularities, 
whereas,  the  other  was  considerably  smoother.   It  would  appear  that  screening 
the  data  is  very  important  in  the  analysis,  and  a  proper  judgement  has  to  be 
made  as  to  whether  a  datum  is  the  reflection  of  an  eddy  or  an  error.  It  is 
hoped  that  satellite  data  will  give  some  insight  in  this  regard. 
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FIGURE  CAPTIONS 


Fig.  1.   The  mean  difference  and  the  standard  deviations  of  the  difference 
between  the  satellite  and  ship  SST  for  the  Pacific  Ocean  and  the 
Atlantic  Ocean  in  1976  (upper)  and  1977  (lower).  The  geography 
is  shown  to  indicate  the  locations. 

Fig.  2.   The  Amalgamated  SST  anomaly  (°C)  for  January  1977.  The  negative 
anomalies  are  shaded. 
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Table  1 


(T 


SAT 


T^-up)  ±  Standard  Deviation 


T*sUJJjL 

!<JtiU<U 

20°-40°     ,  160°E-160°W 

40 

>-60°N   , 

160 

°E-160°W 

1976 
Jan 

-1.85 

+ 

1.35 

.84  ±  1.47 

Feb 

-1.75 

+ 

1.55 

.04  ±  1.17 

Mar 

-1.30 

+ 

1.25 

.19  ±  1.17 

Apr 

-1.20 

+ 

1.2 

.79  ±  1.4 

Hay 

-   .97 

+ 

1.5 

.94  ±  1.2 

June 

-1.19 

+ 

1.85 

1.22  ±1.3 

July 

-   .8 

+ 

1.5 

1.01   ±  1.43 

Aug 

-   .35 

+ 

1.8 

.7     ±  1.18 

Sept 

-   .71 

+ 

1.45 

-   .58  ±  1.25 

Oct 

-1. 

+ 

1.35 

.46  ±  1,31 

Nov 

-1.4 

+ 

1.6 

,91    ±  1.49 

Dec 

-1,08 

+ 

1.5 

.99  ±  1.23 

1977 
Jan 

-2.1 

+ 

1.45 

+ 

.12 

+ 

1.23 

Feb 

-2.5 

+ 

1.3 

- 

.56 

+ 

1.27 

Mar 

-1.9 

+ 

1.4 

- 

.34 

+ 

1.33 

Apr 

-1.95 

± 

1.95 

- 

.38 

+ 

1.53 

May 

-   .585 

± 

1.75 

- 

.07 

+ 

1.43 

June 

+  .385 

± 

2.05 

+ 

.54 

+ 

1.83 

July 

+2.7 

± 

2.3 

+ 

.74 

+ 

1.77 

Aug 

+  .38 

± 

1.7 

- 

.6 

+ 

1.33 

Sept 

-   .32 

± 

1.5 

.55 

+ 

1.27 

Oct 

-    .95 

± 

1.9 

- 

.27 

+ 

1.43 

Nov 

-   .67 

± 

1.8 

+ 

.18 

+ 

1.43 

Dec 

-1. 

± 

1.6 

+ 

.35 

+ 

1.5 

1978 
Jan 

-1.65 

± 

1.55 

.33 

+ 

1.3 
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ANNUAL  AND  INTERANNUAL  VARIABILITY  OF  SST 
IN  THE  EASTERN  TROPICAL  ATLANTIC  OCEAN 

J.  Merle 

CIMAS,  University  of  Miami 
4600  Rickenbacker  Causeway,  Miami,  FL  33149 


In  the  eastern  equatorial  Atlantic  Ocean  the  annual  signal  is  larger  than 
the  interannual  anomalies  but  its  characteristics  have  important  spatial 
variations  (Fig.  1).   Two  regions  are  specially  affected  by  large  amplitude 
both  of  the  annual  signal  (Fig.  2)  and  of  the  interannual  anomalies:   south 
of  the  Equator  (from  0  to  4°S)  and  E.W.  coast.   In  these  areas  the  thermo- 
cline  shallows  and  even  surfaces  in  summer  lowering  suddenly  the  SST. 
During  the  other  months  a  thin  low  saline  and  warm  mixed  layer  covers  the 
whole  Gulf  of  Guinea  giving  a  long  warm  season  only  cut  by  a  slight  secon- 
dary cold  season  in  December- January .   This  warm  superficial  layer  call 
to  our  mind  a  kind  of  El  Nino  phenomenon.   This  warm  season  is  possibly 
due  to  an  eastward  warm  advection  and  its  local  consequences:   rains  and 
run  off  of  the  great  rivers.   This  advection  could  be  complex  involving 
the  three  branches  of  the  counter-current  system  in  an  alternative  northern 
and  southern  influence  explaining  the  strong  superficial  semi-annual  compo- 
nent responsible  of  a  bimodal  temperature  cycle.   This  warm  and  low  saline 
surface  layer  appears  to  be  the  leading  factor  of  the  SST  variability  in 
the  annual  and  interannual  time  scale  whatever  the  dynamic  process  conducting 
its  motion  and  the  motion  of  the  sub-surface  layers. 

Interannual  anomalies,  as  in  the  equatorial  region  of  the  Pacific  Ocean, 
have  a  large  spatial  extension  at  the  scale  of  the  basin  and  affects  the 
region  for  several  months,  sometimes  more  than  a  year  (Fig.  4).  The 
magnitude  of  these  anomalies  is  intensified  both  at  the  Equator  and  near 
the  E-W  coast.  A  long-term  trend  shows  an  increase  of  SST  from  the  begin- 
ning of  the  century  to  the  sixties  and  a  decrease  until  now.  This  trend 
is  similar  to  the  one  observed  by  Fieux  and  Stommel  (1975,  1976)  both  in 
the  North  Atlantic  Ocean  and  in  the  Arabian  Sea  (Fig.  3). 
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Fig.  1.   Annual  variability  and  interannual  variability  (SD)  of  SST  along 
a  ship  line  crossing  the  Eastern  Tropical  Atlantic  from  20°S-10°E  to  10°N 
-20°W.   Interannual  variability  is  found  to  be  5  times  smaller  than  the 
annual  variability  except  in  the  North-South  transition  area  (8°N)  where 
the  annual  signal  is  minimum  (see  Fig.  2). 


Fig.  2.   Amplitude  in  degree  C  of  the  annual  component  of  the  SST  annual 
signal  in  the  Tropical  Atlantic. 


Fig.  3.   Long-term  trend  of  SST  anomalies  in  MSQ300  (0=10W-0-10S)  compared 
with  the  one  found  in  the  North  Atlantic  and  Arabian  Sea  (Fieux  and  Stommel) 


Fig.  4.   Interannual  anomalies  of  SST  in  several  areas  of  the  Eastern  Tropical 
Atlantic  including  a  coastal  station.   The  coherence  of  the  anomalies  seem 
remarkable. 
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RECENT  VARIABILITY  OF  THE  GULF  STREAM  EXTENSION  FRONT 


C.-P.  Thang  and  Peter  J.  Lamb 

Cooperative  Institute  for  Marine  and  Atmospheric  Studies 

University  of  Miami/NOAA 

Miami,  Florida  33149 


Sea  surface  temperature  (SST)  variations  in  the  vicinity  of  the  boundary 
between  the  Gulf  Stream  Extension  and  Labrador  current  waters  during  1962-76 
have  been  identified  from  a  preliminary  analysis  of  individual  monthly  data 
for  2.5°  latitude-longitude  grid  intersections.   This  area  to  the  south  and 
east  of  Newfoundland  experienced  the  maximum  SST  variance  and  strongest  ave- 
rage meridional  SST  gradients  in  the  Atlantic  north  of  the  Equator  during 
1962-76.   The  present  note  focuses  on  the  variability  of  the  winter  SST  field; 
variations  characteristic  of  other  seasons  appear  less  pronounced  and  systematic 

Time  series  of  the  yearly  minimum  SST  were  established  at  2.5°  latitude 
intervals  along  five  separate  meridional  transects  in  the  region  between  35°- 
65°N  and  25°-60°W.   To  the  south  of  Newfoundland,  and  immediately  to  its  east, 
two  principal  trends  are  evident  (Figs.  1-3).   Along  the  equatorward  side  of 
the  transects  (37. 5°-42. 5°N) ,  warming  of  up  to  3.5°C  occurred  between  1969 
and  1975  or  1976.   Its  magnitude  decreased  eastward.   From  1970-1974,  this 
warming  coincided  with  cooling  of  similar  magnitudes  along  the  poleward  side 
of  the  transects  (45°-52 ,5°N) .   As  a  result,  the  maximum  meridional  SST  gra- 
dient was  strengthened  by  up  to  50%,  but  little  if  any  change  occurred  in  its 
latitudinal  position.   At  the  more  poleward  latitudes  50°N  and  52.5°N,  the 
aforementioned  cooling  commenced  as  early  as  1966. 

Further  east  (Figs.  4-5),  the  two  mid-Atlantic  transects  (longitudes 
28.75°W  and  38.75°W)  reveal  smaller  changes.   There  is  only  very  limited 
evidence  of  post-1969  warming  along  the  equatorward  side  of  these  transects, 
with  the  values  barely  reaching  1°C.   Cooling  of  up  to  1.5°C  occurred  after 
1968  at  most  latitudes  from  42.5°N  poleward,  with  more  pronounced  SST  de- 
creases being  confined  to  extreme  northern  latitudes.   The  latter  produced  a 
marked  poleward  displacement  of  the  maximum  meridional  SST  gradient  after  1970. 

In  contrast  to  the  post-1966  trends  described  above,  the  earlier  years 
of  the  study  period  were  characterized  by  more  irregular  variation  (Figs.  1-5). 
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Arabian  Sea  Summer  Cooling 

This  study  is  focused  on  the  large  scale  interannual  variability 
of  the  western  Indian  Ocean's  response  to  the  Summer  Monsoon.  Saha 
(1974)  and  Shukla  (1976)  have  investigated  the  role  of  SST  in  modu- 
lating the  monsoonal  rainfall  over  western  India.  Other  work  (Duing 
and  Leetmaa,  1978)  has  studied  the  overall  upper  mixed  layer  heat  budget 
in  this  region.  We  wish  to  determine  the  extent  and  variability  of  the 
monsoon  induced  ocean  thermal  contrast. 

We  utilize  satellite  remotely  sensed  surface  temperatures  in  this 
study,  hence  validation  experiments  between  various  satellite  data 
products  and  available  ship  observations  are  an  integral  part.  The 
work  involves  most  aspects  of  remotely  sensed  data  acquisition  and 
processing.  In  the  following  discussion  we  shall  concentrate  on  the 
analysis  aspects. 

The  discussion  is  broken  into  two  sections:  data  preparation  and 
analyses  in  Appendix  A  and  a  discussion  of  pertinent  results  in  the 
present  section.  We  shall  freely  abstract  results  from  the  Appendix; 
the  reader  is  urged  to  study  the  technique  outlined  in  Appendix  A. 

Four  sources  of  SST  observations  are  used  in  this  study:  100  year 
ship  means,  ship  of  opportunity  results,  NESS  First  Guess  Field  (FGF), 
Scanning  Radiometer  (SR) ,  and  NESS  G0SSTC0MP  SR  retrievals.  Our 
intention  is  to  use  ship  measurements  as  a  verification  of  satellite 
retrievals  while  also  (in  the  case  of  100  year  means)  providing  long- 
term  mean  monsoon  SST  response  patterns.  The  study  spans  the  months  of 
April  through  September  for  the  years  1975-1977.  Problems  exist  due  to 
variation  in  data  density  from  year  to  year  within  a  source  of  obser- 
vations and  also  between  different  sources.  (Examples  of  data  density 
variation  are  examined  in  Appendix  A  for  ship  and  FGF  retrieval  percent- 
ages). To  minimize  difficulties  due  to  the  diversity  in  data  density, 
monthly  mean  fields  were  calculated  and  uniform  optimal  interpolation 
techniques  were  applied  by  month  to  each  data  set. 

Validation  of  the  satellite  observations  with  ship  measurements 
consisted  of  comparing  the  interpolated  monthly  maps  of  ship  and 
satellite  observations.  Two  types  of  satellite  observations  are  used: 
the  FGF  are  connoted  as  "monthly"  in  Table  1,  G0SSTC0MP  as  noted. 

Table  1  presents  a  comparison  between  ship  and  satellite  SST 
observations  for  1977.  The  preparation  of  the  data  sets  is  discussed 
in  Appendix  A.  The  standard  deviation  of  the  satellite  and  ship 
observations  averages  approximately  1.5  C.  Mean  differences  are  approxi- 
mately 1  C.  This  should  be  compared  with  Brower  et  al_.  (1975)  who  show 
XMS  variability  between  ship  of  opportunity  temperatures  to  be  ^1.5  C. 
Hence  in  the  Indian  Ocean  satellite  observations  appear  to  be  a  reasonable 
method  of  monitoring  SST  (or  at  least  they  are  no  worse  than  ships  of 
opportunity) . 
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A  common  grid  has  been  defined  for  the  area  93  x  67  elements  (46.5° 
longitude  x  33.5°  latitude,  0.5  degree  square).  The  upper  left-hand 
corner  is  located  at  38°E,  23°N.  Ship  and  satellite-acquired  data  of 
various  types  must  be  mapped  onto  this  basis  set.  Experience  shows  that 
due  to  sampling  vagaries  it  will  be  difficult  to  densely  fill  such  a 
grid  with  the  ship  observations.  Interpolation  of  unequally-spaced 
observations  to  the  regular  grid  is  accomplished  by  an  optimal  inter- 
polation scheme  (Gandin,  1963). 

A  100  year  compilation  of  ship  temperature  (Navy  Atlas  Project: 
Indian  Ocean,  1976)  has  been  used  as  a  climatic  mean  temperature  set. 
These  temperatures  (1°  x  1°  grid)  are  interpolated  two-dimensional ly  with 
a  Bessel  interpolation  scheme  to  0.5°  x  0.5°  resolution.  The  results  are 
shown  in  Figure  1  and  the  interpolation  technique  is  discussed  in  detail 
in  Appendix  A. 

Figure  1  presents  the  climatic,  i.e.,  100  year  ship  mean  temperatures 
for  the  months  April  through  September.  Upwelling  along  the  African  coast 
at  8°N  to  12°N  is  well  resolved,  starting  in  May  and  reaching  a  maximum  in 
August  along  Somali  and  Arabia.  The  cooling  in  the  central  part  of  the 
western  Indian  Ocean  is  well  resolved  with  the  relative  high  temperatures 
being  constrained  to  a  region  southwest  of  India.  Upwelling  along  the 
Indian  West  Coast  is  found  also.  An  abrupt  transition  between  seasonal 
warming  and  monsoon-induced  cooling  is  observed.  These  results  compare 
\/ery  well  with  the  satellite  observations  we  have  previously  shown  for 
1976.  The  maximum  thermal  contrast  occurs  between  the  months  of  April 
and  August  in  the  climatic  data.  Therefore  we  shall  limit  the  discussions 
of  the  multi-year  satellite  observations  to  monthly  differences  between 
April  and  August. 

Figure  2a  shows  the  climatic  mean  result  for  the  April -August 
differences.  Figures  2b,  c,  and  d  illustrate  the  April  minus  August 
differences  for  the  raw  composite  satellite  retrievals  for  the  years  1975, 
1976  and  1977.  Duing,  in  the  previous  section,  presents  results  for  1963. 
Although  the  general  pattern  does  not  change  from  year  to  year,  the 
absolute  magnitude  and  area!  extent  of  given  isotherms  show  large  variation. 
A  relative  low  exists  southwest  of  India  in  all  years  with  a  similar 
isotherm  spreading  to  the  south  and  west.  This  can  be  inferred  from  the 
climatic  patterns  in  Figure  1  by  noticing  the  development  during  the  summer 
months  of  the  large  scale  cooling.  Marked  changes  occur  in  the  northwestern 
quadrant  yearly  and  from  climatology.  Examination  of  the  +2.0  C  isotherm 
shows  greater  than  a  factor  of  2  areal  variation  over  the  time  period 
studied. 

The  interannual  differences  of  Arabian  Sea  summer  cooling  are  large, 
both  in  amplitude  and  area.  In  fact,  they  are  large  enough  to  substantially 
alter  the  heat  budget  calculations  presented  by  Duing  and  Leetmaa  (1978). 
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As  a  consequence,  we  may  speculate  that  the  two  main  contributors  to 
summer  cooling,  i.e.,  cross-equatorial  advection  by  the  Somali  Current 
and  coastal  upwelling,  also  undergo  substantial  interannual  variation. 
We  are  planning  to  intercompare  the  interannual  cooling  results  with 
precipitation  data  from  the  Indian  Subcontinent  between  1975  and  1977 
in  order  to  test  the  model  predictions  of  Shukla,  1975. 


14-4 


<_> 

o 

I-* 

ON 

m 

m 

o 

ON 

CM 

o 

oo 

B 

o 

* 

o 

ON 

ON 

o 

ON 

I— ( 

i—i 

2 

o 

O 

s^ 

SB 

i 


o 

o 

H 

CO 

c 

CO 

CO 

o 

<U 

O 

s 

o 


ON 


co 


vO 
I 


vO 


CN 


O 


CN 
I 


M 
3S 
CO 


ON 


On 
vO 


CO    ~* 


NO 


CM 
CO 


en 


m 


5 


<0 
H 


55 


CO 

(U   * 


•-H     00 

r-^   co 


NO 


co 


co 


vO 


CO 


Os 


33 
CO 


ON 

in 

-* 

00 

CO 

iT\ 

CO 

-* 

m 

CO 

1 

H 

co 
co 
o 

O 


o 

u~> 

CN 

<t 

c 

■— i 

i—i 

O 

o 

oo 

r* 

co 

• 

• 

• 

• 

• 

a> 

* 

VO 

ON 

CO 

o 

00 

s 

o 

1 

' 

i— 1 

1 

1 

■u 

| 

i-i 

a) 

>> 

a 

3 

55 

U 

>N 

c 

r-4 

00 

O 

a 

CO 

3 

s 

3 

£ 

<< 

X 

►-> 

•-> 

< 

14-5 


B 


20°S 


20°  N 


10°N 


40"E 


60°E 


20°S  ** 


Figure  1.  Long-term  variability  for  the  months  April-September  based  on 
100  means  of  ship  observations:  a,b,c,d,e,f  refer  to  monthly  means  for 
April,    May,    June,    July,    August    and   September,    respectively. 
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APPENDIX  A 
(Arabian  Sea  Cooling) 

Data  sources  and  subsequent  processing  are  presented.  Emphasis  is  placed 
on  deriving  an  appropriate  interpolation  algorithm  which  can  be  used  on 
ship  and  satellite  SST  observations.  Observations  are  discussed  and 
followed  by  a  delineation  of  the  interpolation  procedure. 

Ship  Observations 

World  Meterological  Organization  Temperature  Daily  File  (TDF11)  ship 
observations  for  the  grid  area  were  obtained  from  the  National  Climatic 
Center  (NCC)  for  the  months  April  through  September  in  the  years  1975,  1976 
and  1977.  These  observations  were  sorted,  checked  for  gross  discrepancies, 
and  mapped  to  the  common  grid.  Retrievals  by  month  and  year  for  the  grid 
are  shown  in  Table  A.l. 

TDF11  observations  are  mostly  intake  or,  rarely,  bucket  measurements.  The 
statistics  are  highly  variable;  Brower  et^  al".  (1975)  show  that  the 
standard  deviation  between  ships  approaches  1.5  C.  Other  work  (Barnett, 
private  communication)  in  the  central  Pacific  has  shown  these  measurements 
are  generally  warmer  than  the  surroundings,  highly  ship  dependent,  and  are 
biased  by  s   0.4  C. 

Satellite  Observations 

NOAA  Scanning  Radiometer  (SR)  infrared  observations  were  chosen  as  the  data 
source  for  this  aspect  of  the  analysis.  The  retrievals  are  mapped  to  0.5 
x  0.5  grids  and  have  been  archived  daily  for  the  entire  study  period.  Two 
types  of  satellite  observations  are  analyzed:  raw  and  "analyzed".  Both 
products  are  produced  by  NESS  from  the  SR  data.  The  first  comprise  the 
"First  Guess  Field"  from  which  NESS  calculates  the  latter,  GOSSTCOMP 
product. 

The  First  Guess  Field  is  generated  on  a  daily  basis  from  the  SR  and 
Vertical  Temperature  Profiling  Radiometer  (VTPR)  results.  Clear  areas  of 
image  are  identified  and  a  10  x  10  (100  km.  x  100  km.)  matrix  extracted 
about  grid  points.  VTPR  water  vapor  correction  and  cloud  contamination 
corrections  are  applied  and  a  boxcar  average  is  calculated  to  produce  a 
corrected  SST  estimate.  FGF  algorithms  have  changed  several  times  since 
the  field' 8  inception.  The  algorithm  is  documented  by  Leese  et^  al. 
(1971a).  These  changes  represent  another  bias  which  must  be  removed. 
Retrievals  by  month  are  given  in  Table  A. 2. 

GOSSTCOMP  (Global  Operational  SST  Composite)  derivation  is  discussed  by 
Leese  et^  £l_.  (1971a).  Essentially,  the  FGF  is  combined  with  previous  day's 
observations,  constraints  on  maximum  spatial  gradient  and  allowable  daily 
changes  are  applied  and  the  resultant  field  is  interpolated  to  produce  a 
dense  mapped  SST  grid.  The  constraints  have  also  been  subject  to  change 
over  the  period  of  interest  and  represent  unknown  biases.  Leese  et  al. 
(1971b)  give  reliability  estimates  for  GOSSTCOMP.  Since  all  grid  points 
are  filled  there  always  exists  100%  retrieval.  We  have  discarded  data  more 
than  one  day  old  in  the  monthly  mean  generation  in  order  to  reduce  biases 
due  to  the  interpolation  of  old  data  by  NESS. 
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The  goal  of  this  study  is  the  understanding  of  the  larger  scale,  seasonally 
varying  signatures  in  the  SST  field.  Monthly  mean  data  sets  are 
constructed  for  each  of  the  three  data  types  by  year.  A  total  of  54  data 
matrices,  many  of  them  sparsely  populated,  now  must  be  interpolated  for  map 
making  and  analysis  purposes . 

Interpolation 

Removal  of  the  long-term  mean  does  not  produce  a  non— zero  mean  data  set. 
For  the  three  years  under  study  the  mean  Indian  Ocean  temperature  differs 
from  the  100  year  mean  by  0.2  C  to  1.0  C  depending  on  the  month  and 
geographical  location.  Subsequently  a  mean  square  difference  plane  is 
removed  from  the  remaining  anomaly  and  the  remaining  field  has  zero  mean. 

Gandin  (1963),  Thiebaux  (1973),  have  developed  optimal  interpolation  in 
which  the  mean  square  error  of  the  interpolated  value  is  a  minimum.  Such 
an  interpolation  utilizes  an  auto— co variance  matrix  of  the  data  points  to 
compute  a  set  of  weights  with  which  a  new  value  is  computed  from  it's 
nearest  neighbors.   Functionally  we  define  the  interpolation  as: 

Y  -  BX  A.l.A 

where  , 

B  ■  C  .  C  7.  A.l.B 

oi    1J 

with  C  the  covariance  matrix  of  the  item  with  respect  to  its  neighbors,  C  a 
covariance  matrix  for  the  X  items,  X  a  vector  of  the  surrounding 
observations  and  Y  the  interpolated  value.  We  make  an  assumption  of  field 
homogeneity,  but  not  isotropy.  That  is,  we  assert  that  the  covariance 
matrix  is  independent  of  position  in  the  data  field,  but  that  the  elements 
have  covariances  dependent  on  the  vector  orientation  of  the  point  set  (the 
covariance  matrix  is  independent  of  field  position,  however  the  value  for  a 
given  element  may  depend  on  the  relative  orientation  between  two  points  as 
well  as  their  relative  separation:  interested  readers  are  referred  tc 
Gandin,  (1963)  for  a  complete  discussion). 

Ideally  the  covariance  matrix  should  always  be  computed  from  the  data  set 
to  be  interpolated.  The  long-term  and  satellite  data  set  densities  art 
high  enough  for  this  approach  to  succeed.  However,  the  ship  observatiot 
sets  are  quite  sparse  with  monthly  composite  minima  and  maxima  of  11%  ant 
28%  occupation,  respectively.  The  lack  of  observations  precludes  i 
reasonable  (i.e.  low  variance)  estimate  of  the  monthly  auto— covariance 
matrices  in  the  Indian  Ocean.  Months  taken  by  year  could  be  averaged 
together  or  yearly  averages  computed.  This  approach  was  rejected  on  the 
basis  of  expected  and  observed  high  interannual  variability  in  the  SST  data 
fields  and  their  respective  auto— covariance  matrices. 

Obvious  alternative  candidates  for  the  ship  auto— covariance  fields  are 
matrices  computed  from  the  long-term  or  ship  data  sets.  Sample  covariance 
matrices,  normalized  to  1  at  zero  lag,  so  they  are  auto—correlation 
functions,  for  the  climatic  and  satellite  observations  are  presented  in 
Figure  A.l  for  June,  1977.  They  are  quite  different  with  the  rapid  fall 
off  of  the  long-term  auto—correlation  function  ascribed  to  the  monthly 
averaging  over  many  (>100)  realizations  destroying  the  inherent  fine 
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structure.  Due  to  this  we  "believe  the  satellite  observations  are  more 
relavent  to  the  ship  measurements  and  we  shall  use  the  satellite  covariance 
matrices  to  interpolate  the  ship  observations. 

The  data  set  (54  matrices  +  6  100  year  mean  matrices)  has  been  uniformly 
mapped  and  interpolated.  A  two-dimensional  Gaussian  smoothing  function  is 
uniformly  applied  to  the  data  sets  to  damp  out  high  frequency  spatial 
variability.  For  these  cases  the  function  has  an  "E"  folding  length  of  1.5 
grid  elements  (165  km.). 

A  number  of  non— physical  biases  with  annual  time  scale  exists  in  the  data 
sets;  ship  population  statistics  vary,  and  satellite  sensing  systems  and 
retrieval  and  analysis  algorithms  change.  To  compensate  for  these  changes 
(only  at  the  zero  order)  all  data  sets  are  referenced  to  the  pre— summer 
monsoon  condition  by  removing  the  April  observation  grid.  The  resultant 
data  sets  cover  the  period  May— September,  1975,  1976  and  1977.  We  will 
concentrate  on  the  August  differences  for  the  data  types. 


Table  A.l 
Ship  SST  Retrievals  (%) 

Year/Month   April      May       June        July  August  September 

1975  11.5      11.6      14.4        14.5        15.0        17.7 

1976  17.2      17.3      18.9        19.1        17.8        17.5 

1977  18.4      19.3      25.5        27.9        25.7        21.7 

Table  A. 2 
Satellite  SST  Retrievals  (%) 

Year /Month   April      May       June        July  August  September 

1975  51.8      49.2      28.7        34.4        34.7        37.1 

1976  81.7      91.7      84.5        73.2        75.5        90.0 

1977  95.4      95.1      91.5        80.5        90.5        94.8 
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Figure  A.l.    June  covariance  matrices  for  a)  climatic  mean  and  b)  raw 
satellite  data  for  June,  1977. 


14-10 


References 


Brower,  R.L.,  W.G.  Pichel,  C.C.  Walter  and  T.L.  Signore.  1975.  Current 
status  and  quality  of  global  operational  sea  surface  temperatures 
from  satellite  infrared  data.  NESS  International  report,  10  pp. 

Duing,  W.  and  A.  Leetmaa.  1978.  Arabian  Sea  Cooling;  a  preliminary  heat 
budget. 

Gandin,  L.S.  1963.  Objective  Analysis  of  Meteorological  Fields.  Israel 
Program  for  Scientific  Translations,  Jerusalem,  1965.  242  pp. 

Leese,  J.W.,  W.  Pichel,  B.  Goddard  and  R.  Brower.  1971a.  An  experi- 
mental model  for  automated  detection,  measurement,  and  quality 
control  of  sea  surface  temperature  from  ITOS-SR  data.  Proceedings 
of  the  Seventh  International  Symposium  on  Remote  Sensing  of  the 
Environment,  March  1971,  Ann  Arbor,  Michigan,  Willow  Run  Laboratories, 
University  of  Michigan,  pp.  625-646. 

Leese,  J. W.,  W.  Pichel,  B.  Goddard  and  R.  Brower.  1971b.  Factors 

affecting  the  accuracy  of  sea  surface  temperature  measurements  from 
ITOS-SR  data.  Conference  Preprint  No.  90  on  Propagation  Limitation 
in  Remote  Sensing,  June  1971,  NATO  Advisory  Group  for  Aerospace 
Research  and  Development,  Paris,  France,  NATO,  Paris,  France.  13  pp. 

Saha,  K.  1974.  Some  aspects  of  the  Arabian  Sea  summer  monsoon.  Tell  us, 
26:  464-476. 

Shukla,  J.  1975.  Effect  of  Arabian  Sea  surface  temperature  anomaly  on 
Indian  summer  monsoon:  a  numerical  experiment  with  the  GFDL  model. 
J.  Atmos.  Sci.,  32:  503-511 

Thiebaux,  H.J.  1973.  Maximally  stable  estimation  of  meteorological 
parameters  at  gridpoints.  J.  Atmos.  Sci.,  30:  1710-1714. 


14-11 


SALINITY  FLUCTUATIONS  IN  THE  LABRADOR  SEA,  AND  THEIR  IMPLICATIONS 
FOR  HIGH  LATITUDE  OCEANIC  CLIMATE  FEEDBACK 

Claes  Rooth 
CIMAS,  University  of  Miami 


SUMMARY:   Hydrographic  data  time  series  from  Ocean  Station  Bravo  suggest  that 
over  a  five  year  period  from  1967  to  1971,  the  normal  winter  time  convection  in 
the  Labrador  Sea  was  shut  off  due  to  low  surface  salinities.   Two  models  of  sea- 
sonal heat  cycles  interacting  with  a  steady  fresh  water  supply  suggest  that  bi- 
modal  state  distributions  (intransitivity)  is  possible  in  some  parameter  range. 
Large  differences  in  the  severity  of  ice  conditions  occur  between  regimes. 

BACKGROUND:   John  Lazier1  of  the  BIO  in  Canada  has  processed  the  hydrographic 
data  time  series  from  Ocean  Station  Bravo.   Figure  1,  produced  by  him,  shows 
a  pronounced  transient  in  surface  conditions  in  the  Labrador  Sea,  where  for 
several  years  in  sequence  surface  salinities  were  too  low  to  allow  deep  convective 
mixing  at  freezing  point  temperatures.   This  period  coincided  with  a  period  of 
extreme  ice  conditions  around  Iceland  and  in  the  Davis  Straits,  and  terminated  in 
the  winter  1971-72,  just  ahead  of  the  Atlantic  leg  of  the  GEOSECS,  a  benchmark 
survey  of  ocean  chemistry  and  tracer  distributions.   As  was  shown  by  Dorsy  and 
Peterson2   concentrations  of  tritium  in  surface  waters  of  the  East  Greenland 
Current  (1972)  were  anomalously  high  compared  to  expectations  based  on  stationary 
Arctic  Circulation.   Derived  from  thermonuclear  tests,  but  also  from  various 
aspects  of  nuclear  technology  applications,  tritium  in  the  Arctic  is  completely 
dominated  by  USSR  activities.   It  is  of  interest  here,  that  testing  of  devices 
in  the  atmosphere  ceased  in  1964.   It  has  been  accepted  since  the  drift  of  the 
FRAM,  that  sea  ice  drift  from  the  Siberian  coast  line  to  the  Greenland  Sea  takes 
a  couple  of  years,  suggesting  that  near  zone  fallout  cannot  be  invoked  as  an 
explanation  of  the  observed  tritium  anomaly.   A  third  suggestive  coincidence  was 
brought  up  at  this  meeting  by  Kukla  ,  who  finds  in  satellite  observations  of 
snow  cover  a  substantial  negative  anomaly  in  the  late  sixties,  abruptly  termina- 
ted in  1971.   According  to  Kukla,  the  anomaly  is  due  to  late  snow  cover  onset  in 
fall,  suggestive  of  an  elongated  runoff  season,  and  perhaps  implying  enhanced 
supply  of  fresh  water  and  terrestrial  fallout  tritium  previously  tied  up  in 
"permanent"  snow  on  the  North  Siberian  slopes. 

The  significance  of  fresh  water  supply  in  the  Arctic  for  the  maintenance  of  the 
ice  regime  has  been  long  appreciated,  c.f.  review  by  Weyl4  .   The  seasonal  vari- 
ations in  the  near  surface  halocline  caused  by  freezing  out  of  salt  in  the  forma- 
tion of  sea  ice  are  thought  in  many  instances  to  limit  the  extent  of  the  latter. 
What  is  explored  here  is  the  effect  of  the  mean  annual  transport  processes  of 
heat  and  salt  in  a  basin  where  convection  may  be  switched  off  by  the  fresh  water 
supply. 

MODEL  PHYSICS  AND  RESULTS.   My  basic  model  is  that  of  a  two-layer  isolated  basin, 
which  is  in  diffusive  lateral  exchange  with  a  large  reservoir  (ocean)  of  heat  flow 
(not  necessarily  of  zero  mean)  and  to  a  fresh  water  inflow,  and  may  be  forced 
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either  by  fluctuations  in  the  latter,  or  in  the  surface  salinity  of  the  outer 
basin.   The  analysis  concerns  itself  with  possible  forms  of  the  stationary  annual 
cycle,  and  intransitivity  of  the  mean  state,  and  is  thus  applicable  to  both  types 
of  climatic  forcing.   The  use  of  a  non-linear  equation  of  state  is  essential,  since 
it  allows  the  generation  of  available  potential  energy  by  the  onset  of  convection. 
It  is  indeed  found  that  for  a  range  of  parameter  combinations,  either  a  seasonally 
convective  state  with  small  salinity  anomaly  in  the  basin  may  exist,  or  alternati- 
vely a  permanent  halocline. 

Provided  that  these  box  model  results  can  be  shown  to  realistically  represent  the 
first  order  effects  in  real  world  marginal  seas,  an  interaction  pathway  between 
high  latitude  terrestrial  hydrology  and  ocean  thermodynamics  has  been  demonstrated 
which  has  a  potential  for  controlling  regional  climatic  anomalies  on  time  scales 
of  several  years  and  larger  . 
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The  interannual  variability  of  heat  exchange  across  the 
air-sea  interface  has  been  investigated  using  9  years  of  data 
collected  by  ships-of-opportunity  operating  in  the  North 
Pacific . 

Figures  1  and  2  show  the  seasonal  values  of  standard 
deviations  and  ranges,  respectively,  for  the  total  heat  flux 
across  the  interface,  the  evaporative  flux,  and  the  solar 
insolation.   These  graphs  reveal  several  interesting  aspects  of 
heat  flux  interannual  variability: 

1)  The  largest  standard  deviation  and  range  values  occur 
in  spring  and  fall,  while  the  lowest  values  occur  in 
summer. 

2)  The  increase  in  total  heat  flux  variability  from  winter 
to  spring  is  due  to  an  increase  in  insolation 
variability,  while  that  of  evaporative  heat  flux  has 
decreas  ed . 

3)  The  decrease  in  total  heat  flux  variability  from  spring 
to  summer  is  due  to  a  decrease  in  the  variability  of 
both  insolation  and  evaporation. 

4)  The  increase  in  interannual  variability  of  total  heat 
flux  from  summer  to  fall  is  due  to  an  increase  in  the 
variability  of  both  insolation  and  evaporation. 

5)  The  decrease  in  variability  of  total  heat  flux  from 
fall  to  winter  is  due  to  a  decrease  in  the  variability 
of  insolation. 

6)  The  maximum  interannual  variability  occurs  in  the 
western  North  Pacific  in  all  seasons  except  summer  when 
it  is  equivalent  in  both  the  eastern  and  western 
Pacific. 
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1)  The  western  area  received  or  lost  heat  at  near- normal 
rates,  while  the  eastern  area  exchanged  heat  with  the 
atmosphere  at  near-normal  rates  until  October  1976;  it 
then  began  losing  heat  at  below-normal  rates  (anomalies 
of  40-80  cal/cm  /day)  through  February  1977.   The 
differential  anomalous  cooling  between  the  two  areas 
during  this  period  is  reflected  in  both  surface  and 
sub-surface  ocean  thermal  data. 

2)  The  area  averages  of  total  heat  flux  (Q7-)  are 
predominantly  related  to  those  of  latent  heat  flux  (0^) 
and  somewhat  less  to  those  of  sensible  heat  flux  (Qj). 
(Correlation  coefficients  between  time  sequences  of  Q  j- 
and  Qfi  and  Q^  are  0.99  and  0.97,  respectively.) 
Fluctuations  in  the  solar  insolation  (Q^O  had  little 
effect  on  the  total  heat  flux  anomalies  for  these  two 
open  ocean  areas. 

3)  The  fluctuating  nature  of  heat  exchange  in  the  western 
area  reflects  the  variations  in  position  of  the  upper 
level  trough  and  storm  tracks  throughout  the  period 
(the  most  intense  anomalous  ocean  cooling  occurred  to 
the  west  of  180°  on  the  western  flank  of  the  intense 
Aleutian  low) . 

4)  In  the  eastern  area,  the  almost  total  absence  of  storm 
tracks  after  September  1976  caused  below-normal  latent 
and  sensible  heat  loss  and  led  to  above-normal 
temperatures  in  the  ocean  surface  layer  down  to  60 
meters . 
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Figure  1 


Figure  2 


Figure  3 

Figure  4 

Figure  5 

Figure  6 
Figure  7 


Seasonal  values  of  the  standard  deviations  for  total 
heat  flux  (Qr),  evaporative  heat  flux  ( Q^ )  and 
insolation  ( Q^ )  for  the  western  and  eastern  mid- 
latitude  North  Pacific.   Seasonal  values  were 
computed  over  the  9-year  period  1969-1977. 

Seasonal  values  of  the  ranges  for  total  heat  flux 
(Qj.),  evaporative  heat  flux  ( Q^> )  and  insolation  (Q  •) 
for  the  western  and  eastern  mid-latitude  North 
Pacific.   Seasonal  values  were  computed  over  the 
9-year  period  1969-1977. 

Area  averages  of  anomalous  total  heat  flux  (Q7-) 
between  ocean  and  atmosphere.   Areas  1  and  2  are 
shown  in  Figure  7« 

Area  averages  of  anomalous  evaporative  heat  flux  (Qg  ) 
between  ocean  and  atmosphere.   Areas  1  and  2  are 
shown  in  Figure  7« 

Area  averages  of  anomalous  sensible  heat  flux  (Qy) 
between  ocean  and  atmosphere.   Areas  1  and  2  are 
shown  in  Figure  7* 


Area  averages  of  anomalous  insolation  (Q^*) 
1  and  2  are  shown  in  Figure  7. 


Areas 


Areas  1  and  2  over  which  the  heat  flux  anomalies  in 
Figures  3  through  6  were  averaged. 
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A  25-year  Oceanic  Flux  Data  Set  of  the  Atlantic  Ocean  and  Analyses 

Andrew  F.  Bunker 

Woods  Hole  Oceanographic  Institution 

Woods  Hole,  Massachusetts  02543 

Energy  fluxes  at  the  surface  of  the  Atlantic  Ocean  from  80°N  to 
60° S  and  averages  of  meteorological  variables  have  been  computed  on  a 
monthly  basis  for  the  25-year  period  from  1948  through  1972.   Fluxes  of 
sensible,  latent,  and  radiational  heat  and  momentum  were  calculated 
using  bulk  aerodynamic  equations  with  variable  exchange  coefficients 
and  empirical  radiation  equations.   The  data  used  were  the  surface 
marine  observations  collected  and  recorded  on  magnetic  tape  by  the 
National  Climatic  Center.   Averages  of  fluxes  and  variables  were  formed 
over  complete  Marsden  Squares  (10X10  degree  squares)  for  each  month  dur- 
ing the  period.   These  averages  together  with  identification  and  number 
of  observations  have  been  recorded  on  magnetic  tape.   The  list  of  var- 
iables and  fluxes  is  shown  on  Fig.  1.   The  tapes  are  available  to  all 
at  cost  of  duplication.   In  addition  to  monthly  averages  other  tapes 
contain  statistics  describing  the  data  set.   These  include  means,  stan- 
dard deviations,  third  order  orthogonal  polynomials,  and  other  quanti- 
ties by  month,  season,  year,  and  entire  set. 

Exploratory  analyses  using  various  statistical  techniques  have 
been  carried  out  in  preparation  for  a  systematic  study  of  the  data  set. 
Figure  2  presents  the  25-year  average  of  the  annual  net  heat  gain  by 
the  ocean  which  depicts  the  energy  exchange  regimes  of  the  ocean.  The 
purpose  of  the  study  is  to  determine  the  physical  processes  producing 
the  observed  temporal  and  spatial  variations  of  the  fluxes  and 
variables . 

A  survey  of  charts  of  standard  deviations  of  fluxes  and  variables 
quickly  pin-points  regions  in  which  the  variability  is  large  and  sug- 
gests relations  between  components  of  the  system.   Figure  3  shows  the 
variation  of  the  oceanic  heat  gain  during  winter.   A  region  of  high 
variability  extends  from  the  Norwegian  Sea  southwestward  to  the  Carib- 
bean Sea.  Much  of  this  variability  is  caused  by  the  variability  of 
the  meridional  wind  component  as  shown  in  Fig.  4.  Variability  of  the 
meridional  and  zonal  winds  in  the  neighborhood  of  continents  result  in 
large  air-sea  temperature  difference  variability  as  shown  in  Fig.  5. 
Other  charts  available  but  not  presented  here  show  contributions  by 
variability  of  cloud  cover,  pressure  and  sea  temperature.   In  the 
tropics  where  the  trade  winds  blow  with  little  variation  the  heat  ex- 
change variation  is  produced  primarily  by  variations  in  cloudiness  and 
sea  temperature  variations  presumably  caused  by  upwelling  and  current 
variations.   In  the  westerlies  of  South  America  during  the  southern 
hemisphere  summer  the  largest  variation  appears  in  a  region  with  few 
observations  and  hence  is  suspect.   If  the  variation  is  real,  other 
charts  show  that  it  is  caused  by  greater  variability  in  cloudiness, 
meridional  winds  and  air-sea  temperature  differences. 

Anomalies  of  the  heat  gain  by  the  ocean  during  January  1950  are 
presented  in  Fig.  6  to  show  the  large  scale  and  continuity  of  patterns. 
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Charts  for  other  months  show  similar  large  scale  characteristics  al- 
though the  patterns  are  dissimilar.   A  mean  sea  level  pressure  chart 
has  been  inserted  in  the  lower  left  to  demonstrate  the  relationship 
between  pressure  pattern  and  heat  gain  through  the  advection  of  warm 
or  cold  air. 

Time  variations  of  fluxes  have  been  explored  by  a)  plotting  the 
spatial  distribution  of  linear  trends,  b)  spectral  and  cospectral  anal- 
yses of  time  series  of  heat  gain  and  variables,  and  c)  computing  eigen- 
vectors and  plotting  time  series  of  their  combination  vectors.  Pig.  7 
exhibits  the  linear  trend  over  the  25-year  period  of  the  oceanic  heat 
gain  by  the  ocean  in  W  m  yr~  .   It  is  seen  that  most  of  the  North  and 
South  Atlantic  has  a  small  negative  trend.   Only  in  the  northwestern 
regions  of  the  North  Atlantic,  regions  off  the  west  coast  of  Africa, 
and  a  small  region  inside  the  horn  of  South  America  is  there  a  positive 
trend . 

Many  fluxes  and  variables  have  been  analyzed  to  determine  spectral 
and  cospectral  variations.   Figure  8  presents  the  spectra  of^ oceanic 
heat  gain,  meridional  winds,  and  latent  heat  flux  in  Marsden  Square  115. 
The  Gulf  Stream  flows  eastward  through  the  northern  section  of  the  square 
and  the  Sargasso  Sea  occupies  the  southern  part.   The  energy  density  of 
heat  gain  and  latent  heat  flux  show  little  or  no  slope  for  periods  from 
2  months  to  5  years  but  it  increases  toward  longer  periods.   The  merid- 
ional winds  are  quite  variable  throughout  the  spectral  range  and  have  a 
trend  toward  lower  energy  density  in  the  long  periods .   In  the  2  to  5 
year  period  range  the  peaks  and  troughs  of  the  heat  gain  and  latent  heat 
flux  vary  in  unison  while  the  meridional  winds  are  out  of  phase.   This 
negative  covariance  is  shown  more  clearly  in  Fig.  9,  the  cospectrum  of 
the  meridional  winds  and  the  latent  heat  flux. 

Several  time  series  for  a  few  areas  which  have  the  greatest  varia- 
bility have  been  analyzed  by  the  empirical  orthogonal  function  tech- 
nique.  Presented  here  are  some  results  for  January  time  series.   In 
Fig.  10  the  upper  left-hand  chart  presents  the  covariance  of  the  oceanic 
heat  gam.xn  Marsden  Square  116  with  the  gain  in  other  squares.   It  is 
seen  that  squares  in  the  Gulf  of  Mexico  and  containing  the  Gulf  Stream 
have  high  covariances.   Cold  shelf  and  slope  water  off  Nova  Scotia  have 
a  smaller  covariance  while  squares  in  the  eastern  Atlantic  have  negative 
covariance  values.   The  other  charts  show  the  first  5  eigenvectors.  The 
percentages  written  to  the  upper  right  of  each  chart  give  the  percentage 
of  the  variance  explained  by  each  eigenvector.   A  physical  interpreta- 
tion can  be  given  to  some  of  the  eigenvectors.   The  first  appears  to  be 
associated  with  the  warm  Gulf  Stream  and  possibly  the  Brazil  Current 
off  Argentina.   The  second  appears  to  be  more  regional  than  physical 
since  it  has  high  values  in  the  northeastern  Atlantic  and  the  equator- 
ial regions.   The  third  eigenvector  seems  to  be  associated  with  the  cold 
Labrador  Current  and  cold  water  upwelling  off  West  Africa,  along  the 
equator,  and  off  Southwest  Africa.   The  fourth  eigenvector  is  puzzling 
in  that  it  has  large  negative  values  over  most  of  the  South  Atlantic. 

Combination  vectors  associated  with  eigenvectors  1  and  3  are  the 
only  vectors  which  show  any  readily  discernable  trend  patterns.   In 
Fig.  11  it  is  seen  that  the  first  vector  series  is  positive  from  the 
beginning  until  1955,  reaches  a  minimum  in  the  mid-60' s,  then  slowly 
rises.   The  third  series  starts  positive,  drops  rapidly  to  a  low  in 
1952,  becomes  positive  in  1963,  reaches  a  peak  in  1967  and  then  de- 
creases. 

As  noted  earlier,  a  more  systematic  study  of  the  data  set  will  be 

carried  out  in  the  future. 
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Net  Annual  Heot  Gain 
by  the  Ocean 
W  m-2 
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Figure  2.  Net  annual  heat  gain  by  the  North  and  South  Atlantic  Ocean 
in  Watta  m"2. 
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Figure  3.   Standard  deviation  of  net  wintertime  heat  gain  by  the  ocean 
in  W  m"2. 
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Figure  4.  Standard  deviation  of  meridional  wind  component  in  winter 
in  m  s~  . 
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Figure  5.   Standard  deviation  of  wintertime  air-sea  temperature  differ- 
ence in  °C. 
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Figure  6.  Anomaly  of  net  heat  gain  by  the  Ocean»January  1950  in  W  m" 
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Figure  7.  Linear  trends  in  oceanic  heat  gain  in  W  m'2  yr"1. 
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Short  Review  of  the  Evidence  for  a  Large  Poleward 
Transport  of  Heat  by  the  Oceans 

Abraham  H.  Oort 

Geophysical  Fluid  Dynamics  Laboratory/NOAA 

Princeton  University 

Princeton,  New  Jersey  08540 


The  various  methods  to  determine  the  north-south  oceanic  heat  flux 
can  be  grouped  into  four  categories.  They  are  (1)  the  direct  method,  (2) 
model  methods,  (3)  the  surface  heat  balance  method,  and  (4)  the  planetary 
energy  balance  method.  The  characteristic  features  of  each  method  are 
summarized  in  Table  1,  and  typical  results  can  be  found  in  Fig.  1. 

Using  each  of  the  first  three  approaches  transport  values  are  found 
to  range  generally  between  1  and  2  x  10^5  watt  in  tropical  latitudes. 
However,  the  fourth  method  (see  Fig.  2)  leads  to  much  higher  values  of 
more  than  3  x  1015  watt  at  low  latitudes  for  the  annual  mean,  and  a  sea- 
sonal amplitude  variation  of  nearly  5  x  1015  watt.  In  the  author's 
opinion,  this  serious  discrepancy  cannot  be  explained  by  errors  in  the 
planetary  energy  balance  computations  alone,  and  the  first  3  methods 
probably  underestimate  the  real  transport.  Some  of  the  possible  reasons 
are  alluded  to  in  Table  1  (see  column  labelled  "difficulties"). 

Results  from  (still  primitive)  oceanic  general  circulation  models 
(Bryan,  1978a)  suggest  the  possible  importance  of  zonal  mean  overturnings, 
often  associated  with  Ekman  drift.  This  gives  some  guidelines  for  a 
possible  direct  measurement  program. 

The  time  seems  ripe  to  embark  on  pilot  studies  to  directly  measure 
in  situ  the  oceanic  heat  flux.  A  good  candidate  project  might  be  to 
sample  the  Atlantic  Ocean  intensively  at  a  particular  tropical  latitude 
for  an  extended  period  of  several  years.  Besides  large  seasonal,  also 
sizeable  interannual  variations  in  the  oceanic  heat  transport  may  be 
expected. 
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OCEANIC  MERIDIONAL 
HEAT  TRANSPORT 
VALUES 
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Fig.  la.  Various  heat  transport  values  for  the  North  Atlantic  Ocean 
and  the  Northern  Hemisphere  after  Baker  (1978).  The  lower 
four  curves  are  based  on  model  method  2b,  the  Sellers/ 
Budyko  curve  on  the  surface  heat  balance  method  3,  and  the 
two  upper  curves  on  the  planetary  energy  .balance  method  4. 
(The  differences  between  the  last  2  curves  are  mainly  due 
to  improvement  of  the  satellite  radiation  estimates  since 
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F1g.  lb.    Northward  heat  transport  in  the  North  Atlantic  and  Pacific 
Oceans  as  calculated  by  Schopf  (1977)  using  model  method  2  c 
based  on  the  observed  3-dimensional  density  structure. 
Geostrophic  and  Ekman  fromulations  are  used  to  determine  the 
shear  flow  component,  and  a  barotropic  model  calculation  to 
define  the  vertical  mean  flow  component. 
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Fig.  lc.    Northward  heat  transports  as  calculated  by  Takano  ejt  al_. 

(1973)  and  by  Bryan  and  Lewis  (1978)  using  model  method  2  c 
with  only  ocean  surface  conditions  specified  from  observa- 
tions. The  curves  represent  average  conditions  for  two 
steady-state  (viscous  -  no  transient  eddies)  January  and 
July  model  runs. 
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ENERGY  BALANCE  FOR  POLAR  CAP 


(a)  ANNUAL  MEAN  CONDITIONS 


To=-Fta-Ta 


(b)  SEASONAL  MEAN  CONDITIONS 


V=-FTA-TA+(SA+S0+SL+S,) 


Fig.  2.    Schematic  diagram  of  the  energy  flow  in  the  Northern 

Hemisphere  climate  system  for  annual  mean  (a)  and  shorter 
time  (b)  averaged  conditions.  Fx/\  denotes  the  net  radiative 
flux  at  top  of  the  atmosphere,  Ta  and  To  the  horizontal 
transports  of  energy  in  the  atmosphere  and  oceans,  resp., 
and  S/\,  Sj,  S|_,  Sq  the  rates  of  storage  of  energy  in  the 
atmosphere,  cryosphere,  land  surface  and  oceans,  resp. 
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THl   HEAT  gALANCE  OF  A  01/WSTIC  ••'•::•    q:    tl.:e:   :•;:.:.--.■;    —l ANTIC 
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Geophysical  Fluid  DynamicE  Laboratory/NCAA 

Princeton  University 

Princeton,  New  Jersey  08540 


Due  to  a  lack  of  suitable  dots  on  ocean  currents,  the  role  of  the 
ocean  in  the  poleward  transfer  of  heat  is  almost  always  es time- ted  in- 
directly. Most  estimates  are  carried  out  using  the  heat  balance  of 
the  ocean  surface.   Recently,  Vonder  Haar  and  Oort  (1973)  and  Oort  and 
Vender  Haar  (1976)  have  used  the  heat  balance  of  the  entire  atmosphere 
and  data  on  ocean  heat  storage  to  determine  the  transport  of  heat  by 
ocean  currents.   Before  planning  new  field  studies  to  obtain  the  data 
^--bdeti   for  direct  r  eas  'j'  ■-:  encs  of  heat  transport  it  is  very  important  tc 
make  sure  that  all  the  information  possible  has  been  extracted  from 
existing  hydrographic  data. 

With  this  in  rind  a  new  diagnostic  method  has  been  developed  for 
computing  large  scale  ocean  currents.   Traditional  methods  for  computing 
ocean  circulation  are  based  on  geostrophic  calculations  from  temperature 
and  salinity  measurements,  and  determining  the  total  transport  by  choosing 
an  arbitrary  reference  level  or  solving  a  vertically  integrated  vorticity 
equation.  Diagnostic  calculations  of  the  latter  type  have  been  mace  by 
Sarkisyan  (1971)  and  Holland  and  Hirschman  (1972).   The  difficulty  with 
such  methods  is  that  the  velocity  structure  of  the  solution  is  not 
compatible  with  the  conservation  of  heat  and  salinity  as  shewn  by  Schopf 
(1977).   Better  success  has  been  achieved  with  a  new  method  which  includes 
the  conservation  equations  for  potential  temperature  and  salinity,  as 


19-1 


well  as  the  equations  of  motion.  The  solutions  are  constrained  to  be  in 
the  close  neighborhood  of  observed  values  of  temperature  and  salinity  by 
adding  a  Newtonian  damping  term  proportional  to  the  local  departure  of 

the  solution  from  data  given  by  Levitus  and  Oort  (1977). 

o      o 

Calculations  are  carried  out  on  a  l.OxlJD  grid  for  the  North  Atlantic 

with  12  levels  in  the  vertical.  Numerical  integration  with  respect  to 
time  of  the  time-dependent  model  leads  to  a  steady  state  in  a  relatively 
short  period  of  time  (one  or  two  years).  The  following  conclusions  may 
be  drawn  about  the  heat  and  fresh  water  transport  in  the  model. 

a)  Mesoscale  eddies  must  have  an  effective  diffusion  of  more 

3  2 
then  10  m  /S  to  be  an  important  mechanism  in  poleward  heat 

transport.  At  low  latitudes  the  dominant  effect  is  over- 
turning in  the  meridional  plane  due  to  Ekman  transport  poleward 
below  the  trade  winds. 

b)  Horizontal  gyres  due  to  the  Gulf  Stream  and  Norwegian  Current 
are  less  effective  in  transporting  heat.  At  high  and  middle 


latitudes  mesoscale  eddies  are  probably  important, 
c)   The  transport  of  fresh  water  amounts  to  0.5x10  metric  tons/S, 
which  is  in  rough  agreement  with  atmospheric  water  balance 
data. 
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Low-Latitude  Seasonal  Ocean  Heat  Transport 


Paul  S.  Schopf 
Laboratory  for  Atmospheric  Sciences  (GLAS) 
Goddard  Space  Flight  Center 
National  Aeronautics  and  Space  Administration 


In  Hellerman's  new  compilation  of  global  wind  stress  on 
a  monthly  basis,  it  is  noticed  that  the  annual  cycle  of  mean 
zonal  wind  stress  has  a  strong  signal  near  the  equator,  which 
may  be  roughly  approximated  by: 


A  sin  (10)  cosf-cOt). 


,o 


where  1,=  Tf  /22   and  o>  =  21T/1  year.   This  reflects  the 
strengthening  of  the  low  latitude  easterlies  in  the  winter 
hemisphere,  and  their  weakening  in  the  summer. 

The  oceanic  response  of  this  forcing  consists  of  Ekman 
flow  plus  freely  propagating  planetary  waves.   The  Ekman  flow 
is  symmetric  about  the  equator,  with: 

B  sin  (ly0)  cos  f-o>t) 

V  ^    : n  ^x 

Y  sin  (0) 

showing  a  "sloshing"  motion  in  the  upper  water  column.   Numerical 
model  studies  of  a  linear  ocean  circulation  have  included 
planetary  waves  with  the  Ekman  flow.   They  show  that  the  planetary 
waves  do  not  significantly  alter  the  net  amount  of  heat  transport, 
although  they  change  the  response  locally. 

A  basin  extending  over  100   of  longitude  and  from  50  N 
to  50°S  was  studied.   The  vertical  temperature  profile  was 
chosen  so  as  to  agree  with  the  mean  values  obtained  from  Levitus 
and  Oort's  (1977)  data  for  the  Pacific  (over  the  same  latitude 
span) .   The  Ekman  response  of  this  basin  indicated  an  annual 
cycle  in  heat  transport  of  1.5  x  lO^S  w  across  the  equator. 

Currently  a  three  layer  model  is  under  construction  at 
GLAS  to  enable  us  to  include  nonlinear  effects  and  a  thermo- 
dynamically  active  mixed  layer.   The  above  work  carried  heat 
adiabatically ,  and  allowed  no  net  heat  exchange  across  the 
surface  of  the  ocean  (i.e.,  it  ignores  the  dynamic  effects 
of  the  surface  heat  flux  on  the  flow) .   The  new  model  will  allow 
for  rectification  of  the  heat  flux  due  to  planetary  waves,  as 
well  as  the  generation  of  more  realistic  boundary  currents. 
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On  the  Heat  Budget  of  the  Tropical 
Atmosphere  and  Oceans 

by 

Stefan  Hastenrath 
University  of  Wisconsin 

Abstract: 

Heat  budget  estimates  for  the  global  tropics  are  derived  from  recent  cal- 
culations of  the  oceanic  heat  budget  and  satellite  measurements  of  net  ra- 
diation at  the  top  of  the  atmosphere.  Annual  mean  heat  export  from  the  zone 
30  N  -  30  S,  amounts  to  about  101  x  1014  W  (=  100  units).  Of  this  total  39 
and  61  units  are  performed  within  the  oceanic  water  body,  and  the  atmospheric 
column  over  sea  and  land,  respectively.  In  the  zone  0  -  10  N,  to  which  the 
planetary  cloud  band  ("ITCZ")  is  essentially  limited  throughout  the  year, 
atmospheric  heat  export  reaches  only  13  units,  as  compared  to  an  oceanic  export 
of  18  units  from  the  zone  0  -  10S.  In  particular,  oceanic  export  in  the  belt 
0  -  5  S  alone  contributes  11  units,  being  90  percent  of  the  net  radiative  heat 
gain  at  the  top  of  the  atmosphere  in  this  latitude  zone.  Accordingly,  atmo- 
spheric heat  export  from  the  realm  of  the  "ITCZ"  related  to  hot  tower  mechanisms 
seems  to  play  a  more  modest  role  in  the  global  heat  budget  than  heretofore 
believed.  By  comparison,  oceanic  export  from  the  cold  water  zones  immediately 
to  the  South  of  the  Atlantic  and  Pacific  Equator  emerges  as  an  important  factor 
in  global  energetics. 

Oceanic  meridional  heat  transport  in  the  Pacific  is  directed  from  the 
tropics  into  either  hemisphere;  in  the  Atlantic  it  is  northward  from  high 
Southern  latitudes  all  the  way  to  the  Arctic;  and  it  is  directed  southward  in 
the  Indian  Ocean.  Oceanic  heat  gain  in  the  Pacific  offsets  dificits  in  the 
higher  Southern  latitudes  of  the  Atlantic  and  Indian  Ocean  sectors,  as  well 
as  in  the  Atlantic  as  a  whole.  Meridional  heat  transport  for  all  oceans  com- 
bined is  largest  around  30  N  and  20  S,  where  it  accounts  for  53  and  35  percent 
of  the  total  poleward  transport.  Atmospheric  transport  is  largest  and  ef- 
fects the  bulk  of  the  total  transport  in  mid-latitudes. 

Greatly  different  estimates  of  net  radiation  at  the  top  of  the  atmosphere, 
and  of  oceanic  and  atmospheric  heat  export  must  be  regarded  as  compatible  within 
the  broad  error  limits  indicated  at  present  for  all  three  terms. 
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Biennial  Variations  in  Surface  Temperatures  over  the 
United  States  as  Revealed  by  Singular  Decomposition 


Eugene  M.  Rasmusson,  John  B.  Jalickee 
Phillip  A.  Arkin  and  Wen-Yuan  Chen 

Center  for  Environmental  Assessment  Services 

EDIS/NOAA 


1.  Introduction 

The  evidence  for  quasi-biennial  variations  in  the  surface  air 
temperature  derives  mostly  from  single  station  analyses.   In  this  study 
we  examined  these  variations  over  the  United  States  during  the  period 
1931-75  using  the  dense  array  of  mean  monthly  data  from  the  344  climate 
divisions  which  cover  the  conterminous  48  states.   The  method  of  analysis 
is  known  as  singular  decomposition  (Lawson  and  Hansen,  1974).   In  the 
meteorological  context  it  is  often  referred  to  as  Empirical  Orthogonal 
Function  (EOF)  Analysis  (Lorenz,  1956).   Our  results  will  therefore  be 
described  in  EOF  terminology. 

2.  Basic  EOF  Analyses 

Table  1  shows  the  variance  explained  by  the  EOF  analyses. 

Table  1 
VARIANCE  EXPLAINED  BY  THE  5  LEADING  TERMS  OF  THE  SINGULAR  DECOMPOSITION 


TERM 


DATA  SET  1 

1.  Departures  from         47.7 
Monthly  Means 

2.  Normalized  45.2 
Departures 

3.  Bandpass  Filtered        40.6 
Normalized  Departures 

4.  Hilbert  Singular         47.5 
Decomposition  of  3. 

Above 


23.8 

8.9 

4.5 

2.6 

(71.5) 

(80.4) 

(84.9) 

(87.4) 

20.8 

10.0 

5.3 

3.6 

(66.0) 

(76.0) 

(81.3) 

(84.9) 

26.3 

15.1 

5.9 

3.1 
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Row  1  of  Table  1  lists  the  explained  variance  for  the  EOF  analysis 
of  the  non-seasonal  departures.   Similarly,  the  results  for  the  EOF 
analysis  of  normalized  non-seasonal  departures  are  given  in  row  2.   Row 
3  shows  the  results  for  the  data  set  which  results  from  filtering  the 
normalized  non-seasonal  departures  using  a  band-pass  filter  with  maximum 
response  at  24  months,  and  half  power  points  at  18  and  36  months  (Landsberg 
et  al. ,  1963). 

We  shall  first  summarize  the  results  for  the  normalized  departures 
(row  2).   The  3  leading  EOF's  reflect  the  larger  spatial  scale  variability. 
The  smaller  spatial  scale  variance  (24  percent)  is  concentrated  over  the 
western  U.S.  and  Florida  (Fig.  1),  and  is  mostly  resolved  by  EOF*s  4  and 
5. 

The  spectra  for  the  amplitude  time  series  of  EOF^  1,  2,  and  3  are 
shown  in  Fig.  2.   The  dotted  lines  are  the  5  and  95  percent  confidence 
limits  assuming  the  time  series  to  be  a  first  order  Markov  process  (Gilman 
et  al. ,  1962).   These  limits  should  be  used  only  as  a  general  guide. 

Clearly  evident  in  the  EOF  3  time  series  is  a  quasi-biennial  signal 
near  a  period  of  26  months.   The  only  other  EOF  3  spectral  estimate 
reaching  the  95  percent  confidence  limit  occurs  at  the  third  harmonic  of 
a  26  month  period.   Similarly,  the  single  EOF  2  spectral  estimate  which 
exceeds  the  95%  limit  occurs  at  the  fifth  harmonic  of  a  26  month  period. 
The  spectrum  for  EOF  1  shows  5  peaks  exceeding  or  nearly  reaching  the 
95%  confidence  limit.   Since  the  time  series  were  not  detrended,  the 
very  low  frequency  peak  is  probably  not  meaningful.   The  peaks  at  periods 
of  6.  5  and  4.  3  months  again  correspond  to  higher  harmonics  of  a  26  month 
period.   The  remaining  peaks  at  slightly  more  than  4  and  9  months  are 
unexplained,  and  may  be  merely  due  to  chance. 

The  quasi-biennial  spike  in  the  EOF  3  spectrum  contains  approximately 
10  percent  of  the  total  variance  of  this  time  series.   Since  EOF  3  itself 
accounts  for  10  percent  of  the  total  normalized  non-seasonal  variance,  the 
biennial  spike  represents  approximately  1  percent  of  the  total  non-seasonal 
variance. 

Fig.  3  shows  the  EOF  amplitude  time  series  from  which  the  spectrum 
was  computed,  with  the  solid  line  being  a  12  month  running  mean.   Quasi- 
biennial  variations  can  easily  be  identified  throughout  almost  the  entire 
period  of  record. 

The  EOF  3  spatial  field,  also  shown  in  Fig.  3,  exhibits  two  out-of- 
phase  areas  of  maximum  variance.   It  should  be  noted  that  standard  EOF 
analyses  portray  variations  as  either  in  phase  or  180  degrees  out  of 
phase.   Therefore,  it  is  not  possible  to  represent  features  with  variable 
phase  relationships,  such  as  traveling  waves,  by  a  single  EOF.   The  next 
section  describes  the  results  of  analyses  designed  to  more  clearly  reveal 
the  phase  relationships  of  the  biennial  signal. 
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3.   Hilbert  Singular  Decomposition 

The  normalized  non-seasonal  departures  were  band-passed,  as  previously 
described,  and  the  new  data  set  subjected  to  an  analysis  which  we  refer 
to  as  Hilbert  Singular  Decomposition.   The  variance  of  the  band-passed 
data  set,  represented  in  terms  of  the  fraction  of  the  original  non- 
seasonal  variance,  is  shown  in  Fig.  4.   The  values  range  between  4  and  11 
percent,  and  average  6  percent.   As  might  be  expected,  this  value  is 
considerably  larger  than  the  1  percent  of  the  total  variance  represented 
by  the  EOF  3  quasi-biennial  spike.   A  reasonable  estimate  of  the  variance 
associated  with  quasi-biennial  variations  probably  lies  between  these  two 
values.   EOF  3  (Fig.  3)  and  Fig.  4  show  a  general  similarity,  with  maxima 
over  the  northeast  U.S.,  minima  from  the  northern  plains  to  the  south  central 
Atlantic  states,  and  a  second  maximum  extending  northwestward  from  Texas. 

The  Hilbert  Singular  Decomposition  (HSD)  analysis  is  similar  to  the 
method  used  by  Wallace  and  Dickinson  (1972)  for  analyzing  the  character- 
istics of  tropical  waves,  and  can  be  viewed  as  a  complex  form  of  EOF  analysis. 
In  order  to  perform  this  analysis  an  augmented  matrix  A  (x, t)  is  formed 
from  the  basic  data  matrix  A(x, t)  by  the  following  operation: 

A+(x,t)  =  A(x,t)  +  j  I,  WH(t-t')A(x,t)  (1) 

where  j  =   V^T 

and  W  (t-t' )  is  the  Discrete  Hilbert  Transform,  with  frequency  response 
of  unity  magnitude  and  90°  phase  shift. 

Equation  (1)  follows  directly  from  the  basic  theory  of  analytic  functions, 

Using  the  augmented  matrix  A  (x,t),  a  complex  singular  decomposition 
series  can  be  computed  in  terms  of  phase  and  amplitude  in  both  the  time 
and  space  domains  i.e. 


^V>-^MV<rH  fe^f*  el»*\ 


(2) 


The  amplitude  at  any  point  in  space  and  time  is  then  obtained  by 
multiplying  Pxp£>  while  the  phase  is  obtained  by  adding  wx  +  wt. 

The  variance  explained  by  the  first  two  HSD  components  is  listed  in 
row  4  of  Table  1.   Only  the  spatial  phase  (tox)  and  amplitude  (p  )  fields 
for  the  first  component   (Fig.  5)  are  presented.   The  amplitude  field 
shows  maxima  over  the  northeastern  and  southwestern  U.S.   The  phase  varies 
little  over  the  northeast,  but  shows  a  steady  shift  over  the  middle  third 
of  the  nation,  to  variations  over  the  southwestern  states  which  are  nearly 
180°  out  of  phase  with  those  over  the  northeast. 
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4.  Analysis  of  Selected  Climate  Division  Time  Series 

It  seems  reasonable  to  return  to  the  data  from  individual  climate 
divisions  in  order  to  determine  whether  they  exhibit  the  same  features 
found  in  the  large  scale  patterns  of  the  EOF  and  HSD  analyses.   Accordingly, 
we  examined  the  average  values  for  the  divisions  covering  the  New  England 
states  (northeast)  and  those  covering  southwestern  Wyoming,  eastern 
Utah,  western  Colorado,  and  northwestern  New  Mexico  (southwest). 

Fig.  6  shows  the  autocorrelation  functions  for  the  non-seasonal  depart- 
ures for  the  two  areas.   Fig.  7  shows  the  same  autocorrelation  functions 
computed  from  band  passed  data.   The  unf iltered  data  show  weak  indications 
of  quasi-biennial  variations,  and  from  the  band  passed  data,  the  period  is 
indicated  to  be  around  25  months.   The  stronger  quasi-biennial  signal  is 
indicated  in  the  data  from  the  northeast.   This  is  confirmed  by  the  spectra 
for  the  respective  areas  (Fig.  8).   Although  there  is  considerable  power 
around  26  months  in  the  spectrum  for  the  southwest  area,  one  would  hardly 
attach  significance  to  this  period  on  the  basis  of  data  from  that  area  alone. 

The  phase  relationship  between  the  two  areas  found  in  the  HSD 
analyses  is  confirmed  by  the  cross  correlation  analysis  (Fig.  9)  and 
cross  spectrum  analysis  (Fig.  10).   These  analyses  show  the  variations  in 
each  area  with  periods  around  26  months  to  be  almost  180  degrees  out  of 
phase. 

5.  Comparison  with  Tropical  Stratospheric  Wind  Data 

Without  prejudging  as  to  whether  quasi-biennial  variations  found  in  the 
tropical  stratosphere  and  those  found  in  surface  temperatures  are  related 
features,  we  performed  preliminary  comparisons  between  the  phase  of  the 
quasi-biennial  variations  in  surface  temperature  over  the  U.S.  and  variations 
in  the  Zonal  Wind  Component  at  50  mb  over  Balboa,  C.Z.  The  Balboa  data  were 
kindly  supplied  by  Dr.  James  Angell.   Initial  results  indicate  an  average 
phase  difference  of  approximately  6  months  between  the  EOF  time  series  and 
the  wind  variations  at  50  mb  over  Balboa.  However,  more  detailed  comparison 
failed  to  reveal  any  apparent  correspondence  between  variations  in  the  phase 
and  amplitude  of  individual  oscillations  found  in  the  two  time  series. 

6.  Conclusions 

The  analysis  described  in  2-5  were  obtained  from  data  covering  the 
period  1931-1975,  and  strictly  speaking,  apply  only  for  that  period. 
However,  the  presence  of  quasi-biennial  variations  throughout  the  entire 
period,  as  illustrated  in  Fig.  5,  suggest  this  to  be  a  relatively  stable 
climate  feature.   It  might  also  be  noted  that  we  failed  to  find  any 
convincing  evidence  of  quasi-biennial  oscillations  in  the  corresponding 
precipitation  data.   We  therefore  draw  the  following  conclusions: 

(a)  Quasi-biennial  variations  probably  account  for  1-5  percent  of 
the  total  non- seasonal  variance  of  surface  temperature  over  the 
United  States. 
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(b)  By  far  the  most  pronounced  spectral  peak  for  periods  from  2  months 
to  20  years,  in  the  EOF  amplitude  time  series,  is  associated  with 
quasi-biennial  variations.   However,  limited  analyses  suggest  that 
the  quasi-biennial  oscillations  may  be  masked  by  local  effects  in 
individual  station  data,  particularly  over  the  western  United  States. 

(c)  Both  EOF  analyses,  and  analysis  of  data  from  individual  climate 
divisions  indicates  mean  phase  differences  of  the  quasi-biennial 
oscillations  over  the  United  States  of  up  to  180  degrees. 

(d)  During  the  period  1931-1975,  there  are  indications  of  an  average 
phase  difference  of  around  6  months  between  the  Balboa,  C.Z.  50  mb 
west  wind  component  and  the  EOF  time  series  for  the  U.S.  surface 
temperature  data.   However,  little  relationship  was  noted  between 
variations  in  phase  and  amplitude  of  individual  oscillations 
exhibited  by  the  different  time  series. 
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Fig.  2.   Amplitude  time  series  spectra  for  EOF's 
1,  2  and  3. 
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Fig.  3.   EOF  3  and  corresponding  amplitude  time  series.   Scale  on 
right  of  amplitude  time  series  is  for  12  month  running 
mean  (solid  line). 
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Fig.  5.   HSD  1  phase  and  amplitude. 
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Autocorrelation  functions  for  average  of  climate  zones  over 
an  area  of  the  northeastern  (above)  and  an  area  of  the  south- 
western U.S.  (see  text). 
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Fig.    7.      Same   as   Fig.    6,    except   for  band-passed  data. 
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Fig.  8.   Spectra  for  northeastern  (above)  and  southwestern  areas, 
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Fig.  9.   Cross  correlation  functions  for  the  northeast  and  southwest  areas 
computed  from  unfiltered  (above)  and  band-pass  filtered  data. 
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Fig.    10.      Coherence   squared   and   phase  difference   for   the  northeast 
and   southwest  areas   computed   from  unfiltered,   normalized 
non-seasonal  departure  data. 
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SUMMARY  OF  NORTHERN  HEMISPHERE  700  MB  HEIGHT  DEPARTURES 

FROM  NORMAL,  1 966- 19^9 


By 


Kirby  J.  Hanson 
NOAA 


Based  on  monthly  700  mb  height  departures  from  normal  (DN)  for  the 
northern  hemisphere  from  Gilman  (l),  zonal  averages  were  calculated.* 
A  13~nx)nth  running  mean  was  applied  to  both  the  original  monthly  grid- 
point  data  and  the  monthly  zonal  values.   Only  13-month  running  mean  data 
are  presented  in  this  paper. 

Zonal  averages  of  700  mb  height  DN  for  the  period  1 965- 1 970  are 
given  in  Figure  1.   The  continuity  of  the  positive  and  negative  anomalies 
with  time  and  latitude  is  apparent. 

Analysis  of  the  regional  distribution  of  700  mb  height  DNs  over 
the  northern  hemisphere  during  the  portion  of  this  period  involving  the 
first  positive  anomaly  is  given  in  Figures  2-k.      A  composite  of  the  bi- 
monthly positions  of  the  positive  poleward  moving  anomalies  is  given  in 
Figure  5. 

The  discussion  at  the  workshop  was  centered  on  the  possible  signifi- 
cance (or  insigni f icance)of  these  map  features. 


References 

(1)  Gilman,  D. ,  1977.  (private  communication) 


-A  grid  of  ^8*4  points  north  of  1 5°N  having  a  separation  of  5°  latitude 
and  10   longitude. 
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Figure 


Legends 


1  Time  variation  of  the  zonal  700  mb  height  anomaly  for 
*»0°N.  Lat.,  1965-70. 

2  Northern  hemisphere  distribution  of  13-month  mean  700 
mb  height  anomalies  centered  on  Jan.  1 966 . 

3  Northern  hemisphere  distribution  of  13-month  mean  700 
mb  height  anomalies  centered  on  July,  1 967- 

4  Northern  hemisphere  distribution  of  13-month  mean  700  mb 
height  anomalies  centered  on  Jan.  1969. 

5  Bi-monthly  position  of  centers  of  positive  700  mb  height 
anomalies,  January,  19&6-  March, 19&9- 
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Variability  in  the  Annual  Fields  of 
Large-Scale  Atmospheric  Water  Vapor  Transport 

Richard  D.  Rosen 
Environmental  Research  §  Technology,  Inc.,  Concord,  Mass 


1.  For  the  last  several  years,  our  studies  at  ERT  have  been  aimed  at 
diagnosing  the  interannual  variability  of  the  general  circulation.  Our 
data  source  has  been  the  General  Circulation  Data  Library,  whose  collec- 
tion was  begun  at  MIT  under  V.P.  Starr  and  which  has  been  added  to  in 
important  ways  by  A.  Oort  at  NOAA/GFDL. 

2.  One  of  our  earliest  studies  along  these  lines  included  an  evaluation 
of  the  annual  mean  meridional  circulation  patterns  in  the  NH  for  each  of 
10  years  (1958-1968).   The  results  are  shown  in  Fig.  1.   Two  features  are 
immediately  apparent:   a  3-cell  pattern  in  most  years;     and  the  exis- 
tence of  two  regimes  characterized  by  the  relative  strengths  of  the  Hadley 
and  Ferrel  cells.   During  years  1-6,  the  Hadley  cell  is  relatively  weak, 
but  during  years  7-10  quite  the  opposite  holds. 

These  results  should  be  viewed  with  some  skepticism,  however.   The 
mean  circulation  is  one  of  the  most  difficult  features  of  the  general  cir- 
culation to  assess  accurately.  Moreover,  the  annual  Hadley  cell  is  the 
result  of  large  cancellations  between  summer  and  winter.   Nonetheless, 
we've  used  a  consistent  scheme  throughout  so  that  the  changes  from  one 
year  to  the  next  may  be  realistic.   If  they  are,  then  we  have  produced 
evidence  that  the  "climate  system"  takes  on  discrete  states  that  may  per- 
sist for  years  and  then  change  relatively  quickly  ("almost  intransitivity"?) 

3.  Motivated  in  part  by  this  result,  we  began  a  study  of  moisture  fluxes 
in  the  atmosphere  by  looking  at  annual  mean  conditions  for  years  drawn  from 
both  circulation  regimes.  Our  interest  in  moisture  fluxes  derives  from  the 
important  role  water  vapor  plays  in  the  energetics  of  the  atmosphere. 

4.  Somewhat  complicating  the  study  of  atmospheric  water  vapor  are  the 
problems  involved  in  accurately  sensing  it.   The  bias  discovered  in  the 
U.S.  instrument  during  the  mid  to  late  60's  is  a  well-known  example.   These 
problems  seem  not  to  impact  strongly  on  our  results,  however,  as  Fig.  2 
shows.   In  it  are  displayed  profiles  of  the  vertically  integrated  amount 

of  water  vapor  in  the  atmosphere  for  each  of  six  years:   years  1-5  from  the 
first  regime  of  meridional  circulation  patterns  of  Fig.  1,  and  year  10  from 
the  second.  Note  that  the  profiles  all  lie  fairly  close  to  each  other, 
although  the  decrease  in  moisture  content  between  1958  and  1963  first  noted 
by  Starr  and  Oort  is  visible,  particularly  in  the  region  10°-25°N. 

5.  Fig.  3  shows  Q^  =  1/g  /qu  dp  for  each  of  six  years  over  the  NH.   Notice 
the  mid-latitude  westerly  fluxes  and  the  tropical  easterly  flow.   The  latter 
achieves  a  local  maximum  in  the  Pacific,  which  is  well  captured  by  the  net- 
work of  island  stations.  Although  this  center  is  fairly  stable  in  years 
1-5,  it  has  shifted  and  weakened  considerably  in  year  10.   Indeed,  as  Fig. 

4  shows,  the  [Qx]  =  1/2tt  /O^d^  profile  for  year  10  is  noticeably  smaller 
in  the  tropics. 
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6.  The  contrasting  behavior  of  year  10  in  the  tropical  Pacific  is  also 
evident  in  the  streamfun ction  field  for  moisture  flow,  the  nature  of  which 
is  largely  determined  by  Q^.  Fig.  5  shows  the  streamfunction  field  for 
each  of  the  six  years.  Note  the  large-scale  organization  of  this  field,  the 
low  wavenumber  in  its  structure,  and  its  reduced  strength  in  the  Pacific 
during  year  10. 

It  is  tempting  to  associate  changes  in  the  tropical  circulation  like 
the  one  we've  depicted  here  to  changes  in  oceanic  forcing.   The  recent 
modeling  study  by  Charney  £j  Shukla  tends  to  support  this  view,  as  well  as 
the  belief  that  year-to-year  changes  in  mid-latitudes  result  from  the  natural 
evolution  of  a  non-linear  system. 

7.  Although  the  other  component  of  moisture  flow,  Q*  =  1/g  Jqv  dp,  is  gen- 
erally smaller  than  Q^  locally,  its  zonal  average  is_one  of  the  major  ele- 
ments of  the  global  energy  cycle.   The  profiles  of  [Q^]  in  Fig.  6  show  once 
again  a  different  character  in  year  10,  in  that  the  tropical  southward 
fluxes  extend  over  a  greater  latitude  belt  than  in  the  other  years.  This 
supports  our  earlier  finding  of  an  expanded  Hadley  cell  during  year  10.  The 
different  behavior  of  [Q^]  in  year  10  can  be  traced  to  strong  southward 
fluxes  which  extend  in  a  belt  from  the  west  coast  of  Africa  to  the  east 
coast  of  Asia  (Q,  values  in  the  tropical  Pacific  are  uncharacteristically 
small  during  year  10,  just  like  Q, ) . 

Finally,  it  is  worth  noting  that  the  range  in  [Qa]  is  of  considerable 
consequence  for  the  workings  of  the  energy  cycle.  For  example,  the  differ- 
ence in  the  moisture  flux  across  30°N  between  high  and  low  values  would,  over  a  year 
imply  a  temperature  change  in  the  atmosphere  north  of  this  latitude  of  some 
25°C.  Certainly  this  was  compensated  for  by  other  meridional  energy  ex- 
change mechanisms,  among  which  are  oceanic  heat  fluxes.  Vonder  Haar  and 
Oort  have  demonstrated  that  these  fluxes  are  not  small  -  we've  shown  here 
that  there  may  in  fact  be  enough  freedom  for  them  to  vary  considerably  from 
one  year  to  the  next.   (It  must  be  borne  in  mind,  of  course,  that  there  are 
other  parts  of  the  energy  cycle  besides  ocean  heat  fluxes  that  could  balance 
the  variability  in  [Qa].) 
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Fig.    1       Cross  sections  of  the  mass  streamfunction  for  the  mean 
meridional   circulation  for  each  of  ten  years.     Year  1  refers 
to   the  period  from  May   1958   through  April    1959,    and  so  on   for 
the  other. years.     Units   are  109  kg/s. 


24-3 


[W] ,  kg  m" 


x  10 

4.5  -i 

4- 
3.5- 

3- 
2.5- 

2- 
1.5- 

1  - 
0.5- 


10°S 


10°N 


20° 


30° 


40° 


50° 


60° 


70° 


—r 

80c 


90^ 


Fig.  2   Meridional  profiles  of  the  zonal  mean  of  W  =  1/g  Jq  dp 
for  each  of  six  individual  years. 
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Fig.    6       Meridional  profiles  of  the  total  moisture   flow  across   latitude 
walls   in   the  atmosphere  for  each  of  the  years   studied. 


24-18 


A  COMPARISON  OF  SURFACE  WIND  PARAMETERS,  CLOUD  AMOUNT, 
AND  SEA  SURFACE  TEMPERATURE  IN  THE  TROPICAL  PACIFIC 

Ellen  J.  Steiner 
CI RES,  University  of  Colorado 

Monthly  means  and  monthly  deviations  from  the  mean  of  sur- 
face wind  divergence,  vorticity,  wind  speed,  cloud  amount,  and 
sea  surface  temperature  (SST)  for  the  period  1965-1972  and  the 
region  25US-25UN,  130-280°  (190°-280c  for  the  SSTs),  were 
examined  and  compared.  The  reason  for  undertaking  such  a  study 
is  to  determine  the  interactions  between  these  variables.  The 
ultimate  purpose  is  to  add  to  the  understanding  and  aid  in  the 
forecasting  of  events  such  as  El  Nino. 

The  wind  data  used  to  compute  the  divergence  and  curl  came 
from  ship  observations  compiled  by  Wyrtki  and  Meyers  (1975). 
The  cloud  data  were  from  daily  nephanalyses  prepared  from  polar 
orbiting  satellite  data,  arid  were  averaged  and  gridded  by 
Sadler  (1969),  and  the  SSTs  were  ship  observations  provided  by 
Forrest  Miller  of  the  National  Marine  Fisheries  Service. 

Time  series  of  correlation  coefficients  between  the  fields 
were  computed,  as  well  as  maps  of  correlation  coefficients  over 
the  entire  time  record  and  for  winter  and  summer.  For  these 
calculations,  the  data  were  interpolated  onto  5  latitude  by  10 
longitude  grids.  The  high  noise  levels  inherent  in  the  data 
rendered  the  correlation  coefficients  systematically  too  low. 
However,  their  significance  was  confirmed  through  their  consis- 
tency in  time  and  their  sampling  distributions. 

Several  distinctive  features  were  manifest  in  the  mean 
fields  (Figure  1).  A  region  of  strong  convergence  and  cyclonic 
vorticity  overlapped  the  deep  near-equatorial  cloudiness  regime. 
A  southern  branch  of  convergence  also  coincided  with  a  cloudi- 
ness maximum,  and  divergence  in  the  northeastern  and  south- 
eastern portions  of  the  field  were  coincident  with  persistent 
stratocumulus.  Both  cloud  and  divergence  fields  were  weakest 
in  March  and  April.  The  SSTs  predictably  decreased  polewards, 
with  the  isotherms  shifting  polewards  in  the  summer. 

The  correlations  in  the  time  series  between  the  cloud  and 
surface  divergence  fields  both  in  the  mean  and  deviation  cases 
were  negative  (Figures  2  and  3),  demonstrating  the  relationship 
between  subsidence  and  fewer  clouds  and  convergence  and  more 
clouds.  There  was  a  striking  hemispheric  difference  in  the 
correlations  between  the  mean  cloud  and  divergence  fields,  with 
the  southern  hemisphere  correlations  far  more  irregular.  Strato- 
cumulus which  exist  in  the  face  of  divergence  systematically 
reduce  the  negative  correlations  in  the  time  series  and  empha- 
size the  importance  of  knowing  cloud  type  as  well  as  cloud 
amount  for  future  analyses. 

The  time  series  correlations  between  the  other  mean  fields 
are  not  discussed  because  meridional  gradients  (especially  in 
the  SSTs)  completely  obscured  the  subtle  relationships  sought 
for  here.  Similarly,  maps  of  the  correlations  between  the 
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mean  fields  mainly  reflected  the  seasonal  oscillation.  In  some 
cases,  the  correlations  were  found  to  reverse  sign  once  the 
seasonal  trend  was  removed.  The  seasonal  oscillation  of  the 
variables  was  in  phase  in  the  region  of  the  near  equatorial  con- 
vergence and  cloud  bands;  that  is,  warm  SSTs  occurred  together 
with  increased  cloudiness,  which  in  turn  was  concurrent  with 
cyclonic  vorticity  and  convergence.  Other  regions  were  more 
anomalous.  An  interesting  example  is  the  mean  SST  vs.  cloud 
amount  correlation  field  (Figure  4)  which  shows  two  large  nega- 
tively correlated  regions  extending  from  the  northwest  and  south- 
east. 

Seasonal  and  meridional  gradients  do  not  appear  in  correla- 
tions between  the  deviation  fields.  These  correlations  will  be 
discussed  next. 

Cloudiness  was  correlated  with  cyclonic  vorticity  (Figure  5). 
There  were  \/ery  weak  positive  correlations  between  anomalies  of 
SST  and  vorticity  when  the  SSTs  were  lagged  by  one,  two,  and  four 
months.  Also,  the  SSTs  were  very  slightly  negatively  correlated 
with  deviations  of  divergence  and  cloud  amount.  In  the  case  of 
SST  vs.  cloud  amount,  two  large  negatively  correlated  reqions 
(similar  tothose  in  the  map  of  the  mean  field  correlations). possi- 
bly dominate  the  signs  of  the  correlations  in  the  time  series 
(Figure  6).  Other  correlation  fields  showed  coherent  regional 
differences  as  well.  Anomalies  of  SST  and  wind  speed  were  slightly 
negatively  correlated,  supporting  the  tenets  that  strong  winds  lead 
to  upwelling,and  also  increase  the  depth  of  the  oceanic  mixed  layer  by 
generating  turbulence.  Anomalies  of  wind  speed  and  cloud  amount 
were  positively  correlated  (Figure  7),  perhaps  because  of  an 
association  between  strong  winds  and  sub-grid  scale  wind  speed 
variability.  Deviations  of  SST  especially,  and  of  cloud  amount 
to  a  lesser  extent,  were  found  to  be  persistent  in  time.  This 
may  account  for  the  persistence  of  the  correlations  when  the  sea 
surface  temperatures  were  lagged  against  the  other  variables. 

There  were  often  regional  differences  between  the  winter  and 
summer  correlation  fields.  The  most  striking  examples  were  the 
seasonal  SST  vs.  cloud  amount  correlations  (Figure  8).  The  win- 
ter correlation  field  shows  substantial  regions  of  positive 
correlations  in  the  east  and  northwest.  By  summer,  the  posi- 
tive correlations  in  the  east  had  contracted  and  weakened  while 
those  in  the  west  had  intensified.  Correspondingly,  the  winter 
minus  summer  correlation  difference  field  was  positive  in  the  east 
and  negative  in  the  west. 

The  relatively  small  correlations  and  in  some  cases  the  weak 
relationships  opposite  to  those  expected  (such  as  those  between 
the  divergence  and  SST,  and  vorticity  and  SST),  while  in  part  a 
function  of  the  noise  level  of  the  data,  may  also  point  towards  the 
existence  of  some  large  scale  forcing  which  overrides  local  effects. 
The  differences  between  correlations  in  the  winter  and  summer 
fields  suggest  that  this  forcing  is  a  function  of  both  season  and 
region.  Though  too  small  to  be  useful  in  forecasts, the  corre- 
lations do  attest  to  the  existence  of  empirical  relationships 
between  the  parameters.  Better  and  more  complete  data  for  recent 
years  will  be  used  to  resolve  some  of  the  questions  left  unanswered 
in  this  study.  25-2 
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Ocean  Variability  and  Climate  Change 

J.  0.  Fletcher 
NOAA/ERL,  Boulder,  Colorado 

The  ponderous  ocean,  integrator  of  atmospheric  behavior,  should  reflect 
climatic  behavior  more  coherently  than  does  the  atmosphere. 

The  100  year  marine  data  set  (1861-1960)  for  the  Pacific  and  Indian 
ocean  areas  is  being  examined  to  identify  the  nature  of  significant  climatic 
changes  and  the  coherence  of  related  sea  surface  temperature  (SST)  anomalies 
in  time  and  space.  Highlights  of  preliminary  findings  are: 

1.  The  SST  field  exhibits  two  main  modes  of  fluctuation. 

2.  One  mode  is  identified  with  short  periods  of  anomalous  warmth 
over  the  eastern  equatorial  pacific  and  occurs  at  irregular  intervals 
(3-15  years).  One  manifestation  of  this  recurring  anomaly  is  the 

"El  Nino"  phenomenon  but  coherent  behavior  is  apparent  over  the  entire 
pacific  basin  and  the  central  Indian  ocean.  Temporal  coherence  is 
strong  from  more  than  one  year  before  the  maximum  to  about  a  half 
year  after. 

3.  The  other  mode  of  SST  fluctuation  appears  to  span  more  than  a 
century  so  the  inference  that  it  is  a  recurring  mode  is  drawn  from 
proxy  data  spanning  longer  periods  (such  as  runoff  of  the  Nile  river). 

However,  its  remarkable  coherence  in  space  and  time  is  reflected  in  the 
1860-1960  data  set  as  follows:  Before  1870  SST  was  generally  warmer  than 
the  century  average  over  most  of  the  oceans,  especially  the  circumpolar  cur- 
rent, eastern  Pacific  and  North  Pacific.  In  this  generally  warm  ocean  there 
appeared  in  1870  a  cold  anomaly  around  southern  New  Zealand,  which  persisted 
and  spread  eastward  and  northward  across  the  south  Pacific  during  the  following 
45  years.  Less  distinct  but  similar  anomaly  patterns  also  developed  in  the 
north  Pacific.  This  cooling  trend  intensified  after  1900,  reaching  a  cli- 
max about  1916,  at  which  time  warm  anomalies  suddenly  appeared  at  the  same 
locations  and  subsequently  spread  in  the  same  patterns  as  did  the  previous 
cold  anomalies.  The  warm  anomalies  appeared  in  1916  (cold  in  1870)  in- 
tensified in  1940  (cold  in  1900)  and  continued  to  the  end  of  the  data  set 
(1960)  averaged  over  the  whole  ocean,  temperatures  were  warmer  before  1900 
and  after  1940,  cooler  between  1900  and  1940.  This  cool  ocean  period  co- 
incides with  the  much  discussed  northern  hemisphere  "climatic  warming".  The 
two  modes  interact.  During  periods  of  cool  ocean  at  0-30°  lat  (1900-1940)  the 
El  Nino  mode  anomaly  occurs  infrequently.  During  periods  of  warmer  ocean  it 
occurs  more  frequently.  From  Nile  runoff  data,  the  inferred  counterpart 
of  the  1900  change  of  climatic  state  is  about  1780  (120  year).  From  SST 
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data  through  1976  it  appears  that  the  cold  anomaly  originating  in  1870 
may  be  repeating  in  1975  (105  year)  but  later  data  is  needed  to  evaluate 
this. 

The  marine  data  set  also  includes  wind  direction  and  force.  The  fol- 
lowing gross  changes  in  the  wind  fields  correspond  to  the  SST  changes 
described.  Before  1870  wind  was  stronger  (than  the  mean)  nearly  every- 
where, about  20%  in  the  circumpolar  westerlies,  over  50%  in  the  western 
Pacific  monsoon  region.  The  main  features  of  the  circulation  were  farther 
north  than  normal  in  both  hemispheres.  After  1870  the  stronger  winds  in 
the  monsoon  region  but  the  westerlies  and  subtropical  highs  remained 
stronger  than  normal  until  about  1900.  After  1900  the  winds  weakened 
generally,  especially  the  southern  hemisphere  westerlies  in  the  Indian 
ocean  sector.  The  weakening  of  the  wind  field  climaxed  during  the  1918-1938 
period,  with  the  westerlies  10-15%  below  normal  in  both  hemispheres.  In 
1940-60  the  wind  field  strengthened  generally  but  the  southern  westerlies 
remained  10-15%  down. 


" 
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Sea-Air  Interaction  at  the  Pole-tide 
Frequency:  A  Preliminary  Report 


T.B.  Starr,  R.A.  Bryson  and  J.E.  Kutzbach 

Center  for  Climatic  Research,  Institute  for  Environmental  Studies 

University  of  Wisconsin,  Madison,  Wisconsin  53706 


Maximum  monthly  sea-level  pressure  anomalies  at  the  pole-tide  (Chandler) 
frequency  have  been  found  by  Bryson  and  Starr  (1977)  to  be  at  higher  mid- 
latitudes  in  the  eastern  parts  of  the  oceans,  especially  the  North  Sea- 
Scandinavia  area.  The  work  of  Kutzbach,  Chervin  and  Houghton  (1977)  suggests 
that  large  anomalies  in  those  areas  can  be  related  to  sea  surface  temperature 
anomalies  in  the  western  parts  of  the  oceans.  Assuming  that  this  mechanism 
is  operative  at  the  pole-tide  frequency,  one  might  hypothesize  that  there 
should  be  a  sea  surface  temperature  variation  at  the  same  frequency.  This 
report  is  on  the  preliminary  results  of  an  investigation  of  the  hypothesized 
relation. 

The  raw  material  for  the  study  was  a  long  set  of  monthly  mean  sea 
surface  temperature  data  by  ocean  squares. 

Fourteen  squares  from  the  Pacific  eastern  and  western  edges  were  chosen 
on  the  basis  of  completeness  and  quantity  of  data  (60  years),  and  fourteen 
more  were  chosen  reaching  across  the  North  Atlantic  from  Labrador  to  the 
North  Sea  (72  years). 

Fast  fourier  transform  spectra  were  calculated  for  the  28  time  series. 
Most  of  the  spectra  showed  more-than-expected  variance  at  both  the  a_  priori 
specified  Chandler  and  Chandler  alias  frequency  bands  (Table  1).  In  a  majority 
of  the  spectra  these  concentrations  of  variance  were  significant  beyond  the 
95%  level  compared  to  a  white  noise  hypothesis  with  "trend" frequencies  longer 
than  1  cy/15  yrs  excluded.  This  seemed  to  be  true  also  when  reddening  of  the 
spectra  was  accounted  for  by  excluding  intra-annual  variance. 

These  results  are  consistent  with  the  angular  momentum  conservation 
mechanism  suggested  by  Bryson  and  Starr  (op.  cit.).  If  that  mechanism  is 
valid  there  should  also  be  a  strong  phase  shift  across  the  oceans  at  the 
pole-tide  frequencies  for  the  mechanism  requires  lateral  shifts  of  water 
masses.  Such  a  phase  shift  was  found.  The  amplitude  of  the  thermal  wave 
at  the  pole  tide  frequencies  is  also  consistent  in  order  of  magnitude  with 
the  calculated  lateral  shifts  and  observed  temperature  gradients. 
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Interannual  Sea  Surface  Temperature  Variations 
As  Related  to  Severe  North  American  Winters 


by 

Arthur  V.  Douglas 

Climatology  Program 
Department  of  Geography 
University  of  Nebraska 
Lincoln,  Nebraska  68588 

The  unusually  cold  winters  of  1976-1977  and  1977-1978 
in  the  central  and  eastern  United  States  have  been  linked, 
inpart,  to  anomalous  sea  surface  temperatures  in  the  North 
Pacific.   The  purpose  of  this  report  is  to  demonstrate  that 
most  of  the  severe  winters  since  1904  have  been  associated 
with  similar  oceanic  states. 

In  order  to  select  winters  with  comparable  conditions 
in  the  United  States .separate  eigenvector  analyses  were  run 
on  monthly  temperature  data  for  December,  January,  and  Feb- 
ruary.  In  all  cases  the  first  eigenvector  of  monthly  temp- 
erature contrasts  conditions  in  the  eastern  United  States 
with  conditions  along  the  west  coast.   This  is  illustrated 
in  Figure  1.   An  examination  of  the  time  coefficients  for 
the  eigenvectors  indicates  that  eigenvector  1  for  December 
is  poorly  correlated  with  eigenvector  1  for  subsequent 
months.   As  a  result  of  this  lack  of  persistence  between 
December  and  subsequent  winter  months ,  it  was  decided  to 
narrow  this  study  to  the  months  of  January  and  February. 

The  past  two  winters  were  characterized  by  high  time 
coefficients  for  eigenvector  1  in  January,  and  the  winter 
of  1977-1978  was  also  characterized  by  an  equally  high 
time  coefficient  for  eigenvector  3.   The  seven  highest 
time  coefficients  for  eigenvector  1  (positive)  in  either 
January  or  February  were  in  1912,  1918,  1929,  1936,  1940, 
1977,  and  1978.   Data  presented  at  the  last  climate  diagnostics 
workshop  indicated  1905  may  have  been  very  similar  to  these 
seven  years,  and  thus  it  is  included  in  tnis  analysis.  The 
January  through  February  air  temperature  anomalies  for  these 
8  years  are  presented  in  Figures  2-9. 

Sea  surface  temperature  anomalies  for  the  Novembers 
prior  to  these  cold  winters  are  shown  in  Figures  2-9. 
Anomalously  cold  sea  surface  temperatures  are  noted  in  the 
subtropical  North  Pacific  prior  to  each  of  the  eight  cold 
winters  in  the  central  and  eastern  United  States.   These 
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cold  water  pools  suggest  enhanced  sea  surface  temperature 
gradients  over  the  subtropical  central  North  Pacific  and 
possible  increased  atmospheric  baroclinicity,  frontogenesis, 
and  cyclogenesis . 

The  relationship  between  subtropical  central  North 
Pacific  sea  surface  temperatures  and  700  mb  heights  is 
illustrated  in  Figure  10.   In  January  below  normal  temper- 
atures are  usually  associated  with  below  normal  700  mb 
heights  downstream.   In  addition,  apparently  colder  than 
normal  sea  surface  temperatures  in  the  subtropical  central 
North  Pacific  are  associated  with  above  normal  heights 
along  the  west  coast  of  North  America.   Farther  downstream, 
a  deepened  continental  trough  is  implied  during  the  cold 
water  years.   The  suggested  teleconnections  in  Figure  10 
are  very  similar  to  the  events  observed  during  the  past 
two  winters . 

A  correlation  of  November  sea  surface  temperatures 
for  the  same  location  with  subsequent  January  through 
February  air  temperatures  (Figure  10)  shows  a  striking 
resemblance  to  eigenvector  1  (Figure  1) .   It  is  suggested 
that  most  of  the  cold  winters  in  the  central  and  eastern 
United  States  were  a  partial  result  of  the  above  air-sea 
interactions  and  downstream  teleconnections. 
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LEGENDS  TO  FIGURES 


Figure  1.   The  first  and  third  eigenvectors  of  January 
temperatures  in  the  United  States.   Period 
of  analysis  1906-1978  for  32  stations. 

Figure  2.   Eastern  North  Pacific  sea  surface  temperature 
anomalies  for  November  1904  (upper)  and  United 
States  air  temperature  anomalies  for  January  through 
February  1905  (lower) .   Sea  surface  temperature 
anomalies  are  based  on  the  23  year  period 
1949-71,  and  air  temperature  anomalies  are 
based  on  the  40  year  period  1938-77. 

Figure  3.   Same  as  Figure  2  except  sea  surface  temperature 
anomalies  are  for  November  1911,  and  air 
temperature  anomalies  are  for  January  through 
February  1912. 

Figure  4.   Same  as  Figure  2  except  sea  surface  temperature 
anomalies  are  for  November  1917,  and  air 
temperature  anomalies  are  for  January  through 
February  1918. 

Figure  5.   Same  as  Figure  2  except  sea  surface  temperature 
anomalies  are  for  November  1928,   and  air 
temperature  anomalies  are  for  January  through 
February  1929. 

Figure  6.   Same  as  Figure  2  except  sea  surface  temperature 
anomalies  are  for  November  19  35,  and  air 
temperature  anomalies  are  for  January  through 
February  1936. 

Figure  7.   Same  as  Figure  2  except  sea  surface  temperature 
anomalies  are  for  November  1939,  and  air 
temperature  anomalies  are  for  January  through 
February  1940. 

Figure  8.   Same  as  Figure  2  except  sea  surface  temperature 
anomalies  are  for  November  19  76,  and  air 
temperature  anomalies  are  for  January  through 
February  1977. 

Figure  9.   Same  as  Figure  2  except  sea  surface  temperature 
anomalies  are  for  November  19  77,  and  air 
temperature  anomalies  are  for  January  through 
February  1978. 

Figure  10.  Upper  map  shows  the  correlation  between  central 

North  Pacific  sea  surface  temperatures  and  700  mb 
heights.   Lower  map  shows  the  correlation  between 
November  sea  surface  temperatures  and  January 
through  February  United  States  air  temperatures. 
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Sensitivity  Study  of  the  Annual  Water  Balance 

Peter  S.  Eagleson 

Massachusetts  Institute  of  Technology 

Cambridge,  Massachusetts   02139 

Soil  moisture  dynamics  may  be  uncoupled  from  soil  heat  flow  and  from 
atmospheric  dynamics  by  the  expedient  assumption  that  the  land  surface  tem- 
perature, Tg,  is  an  independent  climatic  variable.   The  instantaneous  water 
balance  is  then  governed  by  the  conservation  equation 

i  -  eT  -  y  =  9S/9t  (1) 

and  two  soil  moisture  flux  relations 

y  =  y(s,  t;  climate,  soil)  (2) 


and 


e  (s,  t;  climate,  soil,  vegetation),     e  <  e 
em  -  <      T  T    P       (3) 


'T 


e  otherwise 

P 


in  which  i  =  rainfall  intensity;  e  =  evapotranspiration  rate;  y  =  surface 
+  groundwater  runoff  rate;  S  =  water  storage;  t  =  time;  s  =  moisture  concen- 
tration in  the  soil  boundary  layer;  e  =  potential  evaporation  rate. 

Integrating  Equations  (2)  and  (3)  over  the  respective  storm  and  interstorm 
periods  and  using  the  probability  density  functions  for  the  independent  cli- 
matic variables  allows  estimation  of  the  expected  annual  value  of  Eq.  (1). 
This  gives,  for  a  given  climate,  soil  and  vegetation 

E[P  ]  =  E[ET  (s  )]  +  E[Rs  (s  )]  +  E[R   (s  ) ]  (A) 

A  A  8A 

in  which  P.  =  annual  precipitation;  E   =  annual  evapotranspiration;  R   = 

A  TA  SA 

annual  surface  runoff;  R   =  annual  groundwater  runoff;  and  s  =  time  average 

moisture  concentration  in  the  soil  boundary  layer. 


In  Figure  1,  the  bare  soil  evaporation  (normalized  by  the  potential 
evaporation)  is  presented  as  a  function  of  the  dimensionless  climate-soil 
parameter 

2faK(l)V(l)»a    (£+5)/2 

E  -  ri      so  (5) 

7T  m  e 
P 

in  which  3  =  reciprocal  of  mean  time  between  storms;  n  =  effective  soil 
porosity;  K(l)  =  effective  hydraulic  conductivity  of  soil  at  saturation; 
H'(l)  =  (extrapolated)  soil  moisture  potential  at  saturation;  <J)  =  dimension- 
less  diffusivity  of  soil;  m  =  -  ds  /d^(s  );  c  =  dK(s  )/ds  . 
J  0      0  o     o 

Notice  that  as  E  gets  large  the  evaporation,  E   ,  becomes  climatically 

A 
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limited  by  the  potential  evaporation,  E   ,  and  is  independent  of  the  soil 

properties.   On  the  other  hand,  when  E  gets  small,  E   falls  away  from  E 

A  PA 

due  to  limitations  of  water  supply  to  the  soil  surface  either  by  precipitation 
(small  s  as  given  by  Eq.  (4))  or  by  the  inability  of  the  soil  to  deliver 
moisture  to  the  surface.   For  small  E,  evaporation  is  very  sensitive  to  soil 
properties.   This  characteristic  dominates  water  balance  behavior  as  is  illus- 
trated in  the  following. 

For  a  given  set  of  climate,  soil  and  vegetation  parameters,  Eq.  (4)  deter- 
mines s  which  can  then  be  reintroduced  to  evaluate  each  term  of  the  water 
balance.   In  this  way,  the  water  balance  sensitivity  of  bare  soil  (i.e.,  no 
vegetation)  to  changes  in  soil  and  climate  parameters  is  explored  as  is  illus- 
trated in  Figures  2  and  3.   We  will  examine  a  sub-humid  and  an  arid  climate 
having  the  characteristics  given  in  Table  1. 

In  Figure  2,  the  normalized  water  balance  components  of  the  two  climates 
are  compared  for  a  range  of  saturated  intrinsic  soil  permeabilities,  k(l),  and 
of  soil  pore  disconnectedness,  c.   Moisture  moves  most  easily  for  large  k(l) 
and  small  c.   Notice  the  insensitivity  of  E   to  soil  properties  in  the  sub- 

A 
humid  climate  except  where  the  soil  is  impermeable.   In  the  arid  case,  E   is 
extremely  sensitive  to  the  soil  properties.  A 

To  illustrate  the  importance  of  the  potential  rate  of  evaporation,  we 
have  used  the  Santa  Paula  e  with  the  rest  of  the  Clinton  parameters  in 
Figure  3-b.   Notice  that  this  single  change  causes  a  large  qualitative  altera- 
tion in  the  Clinton  water  balance  causing  it  to  closely  resemble  the  arid 
Santa  Paula  case. 

This  illustrates  the  important  role  played  by  soil  properties  in  deter- 
mining local  recharge  of  atmospheric  moisture  and  may  help  to  explain  why  an 
increased  albedo  (such  as  accompanies  overgrazing),  through  its  reduction  of 
e  ,  reduces  the  local  recharge  of  atmospheric  moisture. 

REFERENCE 

Eagleson,  P.  S.,  "Climate,  Soil  and  Vegetation,"  Parts  1-7,  Water  Resources 
Research,  Vol.  14,  No.  5,  Oct.  1978. 

Table  1 
COMPARATIVE  CLIMATES 

Clinton,  Ma.  Characteristic  Santa  Paula,  Ca. 

94.1  cm  E[P  ]  53.4  cm 

0.15  cm/ day  e  0.27  cm/day 

365  days  Season  length,  m  212  days 

109  Number  of  storms,  m  16 

0.3  days  Storm  duration,  m  1.4  days 

0.9  cm  Storm  depth,  m^        r  3.4  cm 
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log   Ml) 


ETA  /  PA  ■  0.49 


log  k(l) 


log[k(l) 


log[k(l)] 


TRSA/PA  '0.76 


RSA/PAi0.82 


log[k(l)] 


log  [k(D] 


rRqA/PA     0.56 


rRgfl/PA     0.62 


log  [k(l)] 


a.    SUB  -HUMID   CLIMATE     (CLINTON,  MA. ) 


b.    ARID    CLIMATE    (SANTA    PAULA,    CA. ) 


FIGURE  2    WATER  BALANCE  SENSITIVITY  TO  SOIL  PROPERTIES 
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Ts0  0.78 


s0=  0.72 


og[k(l)] 


log[k(l)] 


T  Eta  /Pa-  0.45 


-log[k(l)] 


log   Ml) 


log[k(D] 


log  [k  (I )] 


rRgA'PA  °V 


a.    CLINTON    CLIMATE  b.  CLINTON    CLIMATE    WITH   SANTA    PAULA    ep 

EFFECT  ON   ANNUAL    WATER  BUDGET    DUE    TO   INCREASING    AVERAGE     RATE   OF    POTENTIAL    E  VAPOTRANSPI R  AT  ION 

FIGURE   3  WATER  BALANCE   SENSITIVITY   TO  RATE   OF  POTENTIAL   EVAPORATION 
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STUDY  OF  THE  HYDROLOGICAL  CYCLE  BASED  ON 
SATELLITE  AND  SURFACE  INFORMATION 

by 
Yu.V.  Kurilova,  P.A.  Kolosov,  E.L.  Muzylev, 
L.K.  Poplavskaya 

The  study  of  real  processes  of  moistening  anomaly 
formation  is  possible  only  on  the  basis  of  full  description 
of  the  water  exchange  processes  and  active  parameterization  of 
the  hydrological  cycle.  Such  a  parameterization  includes  the 
following  three  steps;  -  (1)  search  for  analytical  or  corre- 
lation relationships  between  internal  parameters  of  the 
hydrological  cycle  and  external  parameters  of  large-scale 
atmospheric  processes  obtained  from  surface  and  satellite 
information;   (2)  formation  of  generalized  informative  para- 
meters of  the  hydrological  cycle  having  properties  of  spa- 
tial generalization;   (3)  consideration  of  regularities  of 
real  processes  within  the  framework  of  space-time  genetic 
boundaries.  Natural  temporal  boundaries  separate  the  main 
seasons  in  the  hydrological  cycle  according  to  the  relation 
between  phase  states  of  water  and  values  of  the  radiation 
balance.  The  establishment  of  spatial  natural  boundaries 
is  made  taking  into  account  the  zones  of  different  moistening, 
and  amounts  to  isolating  the  anomalies  or  areas  of  uniform 
values  of  hydrometeorological  fields  and  accounting  for  their 
total  or  mean  characteristics  via  informative  generalized 
parameters.  The  distribution  of  zones  of  different  moiste- 
ning is  determined  by  the  relation  between  heat  and  moisture 
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which  is  governed  "by  large-scale  atmospheric  processes  and 
defined  from  both  the  thermal  characteristics,  and  subsequent 
precipitation  and  soil  moisture  storage* 

Most  suitable  information  for  obtaining  the  necessary 
generalized  characteristics  of  water  exchange  is  satellite 
information  due  to  its  properties  of  generalization  and 
global  registration.  In  the  proposed  parameterization, 
global  fields  of  cloudiness,  snow  cover,  soil  moistening 
and  anomalies  of  ocean  characteristics  are  used  as  initial 
information. 

Being  the  link  connecting  all  processes  of  the  hydro- 
logical  cycle,  the  moistening  of  the  layer  of  soil  thermal 
regulation  \A/  L^)      is  chosen  as  the  basic  parameter  for 
description  of  the  hydrological  cycle©  Having  more  inertia 
than  dry  soil,  moistened  soil  is  functioning  as  a  "memory" 
in  the  atmosphere  —  underlying-surface-of-the-continents 
system.  The  degree  of  inertia  in  heat  absorption  and  release 
is  proportional  to  soil  moistening.  Since  soil  moistening 
for  a  given  moment  is  the  result  of  the  processes  of  heat 
and  moisture  exchange  for  a  certain  preceding  period,  it 
turns  out  to  be  able  to  affect  the  direction  and  character 
of  subsequent  processes,  thus  playing  the  role  of  a  regula- 
tor for  the  processes  of  heat  and  moisture  exchange  between 
the  atmosphere  and  the  land  surface. 

Parameterization  of  the  hydrological  cycle  is  based 
on  the  water  balance  equation  for  the  layer  of  thermal 
regulation 
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A\'  Li)  ^  \VC  AVeri  +  i  F  Ox)  oJt 


B>(r)  =  x(T)-y(T)-E^; 


d) 


where  Wo  ~~  is  the  initial  soil  moistening  characterizing 
year-to-year  changes  of  soil  moisture;  \\ch   ~  is  the  water 
equivalent  of  snow  cover;  Xt*w  ,  X  ttj  f  and  (r  (ft)  ~~ 
describe  the  changes  of  precipitation,  streamflow  and  eva- 
poration respectively  for  the  period  under  consideration. 
The  thermal  regulation  layer  is  variable  with  time;  for 
its  description  we  used  data  on  moisture  storage  in  the  one- 
meter  layer,  in  accordance  with  measurements  of  the  existing 
agrometeorological  network.  Dynamic  components  of  the  water 
balance  l3  L'X)  in  equation  (1)  determine  the  processes  of 
accumulation  and  discharge  of  soil  moisture  and  form  the 
most  inert ial  water  balance  component,  i#e#  soil  moistening 
\\  (jiz)   •  Its  value  for  the  given  moment  plays  role  of  the 
initial  condition  for  the  subsequent  stages  of  the  hydrolo- 
gical  cycle. 

The  hydrological  cycle  is  different  in  various  cli- 
matic zones*  The  principal  distinctive  features  showing 
up  in  the  inertia  of  processes  are  determined  by  the  phase 
state  of  water  in  winter.  They  can  be  described  by  the 
following  characteristics:  the  relation  between  liquid 
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and  solid  precipitation  and  duration  of  the  seasonal  snow 
cover.  On  the  basis  of  these  characteristics  the  following 
four  zones  are  identified  according  to  the  annual  distribu- 
tion of  surface  albedo:  a  zone  with  settled  snow  cover; 
a  zone  where  snow  cover  is  absent;  and  two  transitional 
zones  with  the  prevalem.ee  of  solid  and  liquid  precipitation 
respectively.  Parameterization  of  the  hydrological  cycle 
according  to  equation  (1 )  was  made  for  zones  with  settled 
and  unsettled  snow  cover* 

In  accordance  with  seasonal  variations  of  the  radia 
tion  balance,  two  basic  periods  are  isolated  within  the  hydro- 
logical  cycle  -i.e.  a  warm  period  with  the  active  turbulent 
heat  and  moisture  exchange,  and  a  cold  passive  period  of  the 
supressed  exchange  processes  (Fig.  la,  lb).  The  principles 
of  estimating  the  water  balance  elements  for  these  periods 
are  different.  In  the  former,  the  dynamic  members  of  equa- 
tion (1)  are  taken  into  account;  in  the  latter,  the  resul- 
ting effect  is  estimated. 

In  the  warm  snowless  period,  according  to  the  rela- 
tion between  accumulated  heat  and  stored  moisture,  and  con- 
sequently according  to  the  character  of  the  interaction  of 
the  processes  with  the  underlying  surface,  the  following 
three  stages  are  isolated:  - 

(1 )  the  initial  stage  of  active  evaporation  in 

spring,  with  the  unlimited  soil  moisture  storage 
and  sharp  increase  of  the  radiation  balance; 
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(2)  the  summer  proper,  with  the  maximum  radiation 
balance  and  prevalence  of  transpiration; 

(3)  the  final  stage  of  attenuation  of  evaporation 
when  at  the  small  radiation  balance,  conditions 
for  autumn  soil  moistening  are  formed. 

In  the  cold  winter  period,  when  the  heat  and  moisture 
exchange  between  the  atmosphere  and  land  does  not  occur, 
due  to  heat -insulating  properties  of  snow,  and  over  ocean, 
it  is  more  active,  the  relevant  stage  of  hydrological  cycle 
may  amount  to  accumulating  oceanic  moisture  in  the  continen- 
tal snow  cover# 

A  genetically  transitional  season  preceding  early 
winter  (when  snow  is  accumulated)  and  a  pre-spring  season 
(when  snowmelt  occurs)  are  characterized  by  the  presence  of 
two  phase  states  of  water  fixed  in  satellite  images  by  inter- 
mittent white  and  black  patches.  The  physical  meaning  of 
the  influence  of  these  seasons  amounts  to  including  self- 
intensifying  mechanisms  with  a  positive  feedback  at  the 
snowline,  which  were  studied  from  satellite  data  by 
G.  Kukla.  The  main  point  of  such  mechanisms  in  snowmelt 
is  an  approximately  4-fold  increase  of  the  radiation  balance 
through  an  uneven  decrease  of  albedo  by  80  to  20  per  cent. 
The  intensity  of  inclusion  of  such  mechanisms  influencing 
the  snowmelt  rate,  may  be  defined  not  only  by  atmospheric 
conditions  but  also  by  the  water  equivalent  of  the  snow 
cover  accumulated  during  the  winter,  i.e.  the  depth  of 
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snow  cover.  Thus,  the  snow  cover  "keeps  memories"  of  the 
character  of  the  preceding  winter.  Similar,  but  weaker, 
feedback  mechanisms  can  be  observed  in  the  period  of  matura- 
tion and  harvesting,  grass  and  leaves  yellowing,  when  albedo 
increases  by  5  to  15  per  cent.  Thus,  for  regions  with  the 
settled  snow  cover,  six  seasons  of  different  duration  with 
the  prevalence  of  different  physical  processes  are  dis- 
tinguished within  the  annual  hydrological  cycle  (Fig.  2a). 
In  regions  with  the  settled  snow  cover,  where  the  entire 
cold  period  represents  an  extended  transitional  season 
(Fig.  2b) f  only  four  genetic  seasons  are  noted.  The  dura- 
tion of  these  genetic  seasons  is  subject  to  changes  from  to 
year  depending  on  the  peculiarities  of  circulation  processes. 

The  transitional  period  of  snow  cover  formation 
(season  preceding  early  winter  )  may  be  considered  as  a 
passive  period  of  the  so-called  suppressed  processes  of 
exchange  with  the  underlying  surface  since  the  evaporation 
from  snow  cover  C  (w  is  close  to  zero,  the  precipitation 

X  L.T)     on  snow  cover  is  realized  only  in  spring,  and  the 
runoff  equal  to  2  to  4  mm  (as  was  approximately  estimated 
by  the  authors  for  the  Don  catchment  area)  may  be  neglected. 
The  downward  flux  from  the  thermal  regulation  layer  may  be 
not  taken  into  account  either  since  they  participate  in  the 
active  exchange  only  in  the  succeeding  warm  period. 
Therefore,  by  the  end  of  the  period  preceding  early  winter 
the  total  value  of  soil  moistening  remains  equal  to  VVo' 

In  the  transitional  snowmelt  period  (pre-spring 
season)  evaporation  participates  unevenly,  its  value,  at 
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first,  "being  equal  to  zero  (over  the  snow  surface),  and  then 
climbing  to  the  maximum  (typical  for  the  spring  season)  over 
extremely  moistened  snow-free  areas.  It  should  be  noted 
that  in  this  period  evaporation  increases  through  both  the 
limiting  value  of  soil  moistening  in  the  spring  flood  period 
and  extention  of  snow-free  area,  and  an  increase  of  the  radi- 
ation balance  intensified  by  the  climb  of  albedo  during  the 
snowmelt  period.  The  part  of  water  which  was  not  evaporated 
distributes  among  the  surface  runoff  VCt)  and  ground-water 
storage  in  accordance  with  the  runoff  coefficient  /]  • 

Dynamic  members  of  equation  (1)  f3>  (/t)  somewhat 
alter  their  contributions  in  spring  and  summer  in  connection 
with  different  relations  between  heat  and  moisture  *•  in 
spring,  the  soil  moisture  storage  is  practically  unlimited, 
at  the  relatively  small  radiation  balance;   and  in  summer, 
when  the  radiation  balance  is  great,  evaporation  is  limited 
by  the  soil  moisture  content.  Transpiration  becomes  signi- 
ficant, its  relation  with  evaporation  being  of  the  order  of 
4  ;  1.  The  above  explains  the  reason  for  a  substantial 
excess  of  evaporation  in  summer  over  that  in  spring  (Fig#  1). 

At  the  moment,  at  which  transpiration  becomes  preva- 
lent over  evaporation  and  excesses  it  significantly  may  be 
considered  as  the  genetic  boundary  of  the  beginning  of  the 
summer  proper.  On  a  series  of  satellite  images  (especially 
multizonal  ones),  where  vegetation  and  soil  are  distinguished, 
it  is  possible  to  identify  the  contours  of  the  vegetation- 
covered  areas,  in  accordance  with  the  above  criterium. 
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Parameterization  of  components  of  the  hydrological 

equation  (1)) 

1«   Initial  moistening  \J\/0C 

It  is  reasonable  to  begin  the  study  of  the  hydrolo- 
gical  cycle  with  the  end  of  the  summer  season  since  under  the 
conditions  of  due  greatest  evaporation  rate  in  summer,  the 
relationships  with  the  antecedent  processes  fixed  by  the 
moistening  of  the  active  layer  are  broken  to  a  great  extent. 
The  minimum  soil  moistening  in  the  end  of  summer  is  assumed 
equal  to  \Woc  and  calculated  using  the  same  equation  (1). 

As  the  zero  approximation,  it  may  be  assessed  using  a  map  of 
soils  taking  into  account  their  least  water  field  capacity  : 
"Woe-  *  <*  Wi-ins 

2#  The  initial  water  equivalent  of  snow  cover  v\ic«. 

VVQi   is  determined  for  the  moment  of  snow  melting 
by  tracing  the  snowmelt  limit  on  a  series  of  subsequent 
images.  \\JCu    is  calculated  from 


"We,-  XUu^rjUf)  (2) 

where  r;   andht(  —  are  areas  where  snowmelt  occurred  during 

successive  periods;  [^   CF)-  is  a  coefficient  describing  the 
snowmelt  rate  which  may  be  considered  proportional  to  the 
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changes  of  snow-free  area  Ar  i  JC  *V  A  F  * 

Temporal  trend  of  values  A  V  obtained  from  satel- 
lites is  in  a  good  agreement  with  the  temporal  trend  of  the 
averaged  surface  temperatures  (Fig.  2).  Then 

Wch.^ZCFw-F"^  (3) 

From  surface  data  on  the  water  equivalent  of  snow  cover  and 
satellite  data  on  the  dynamics  of  areas  where  snowmelt  had 
occurred,  the  map  of  the  value  |c   was  prepared  for  the 
south  of  the  European  territory  of  the  country  (Fig.  3)» 
The  coefficient  k.  varied  within  the  range  of  the  order 
increasing  regularly  in  the  direction  from  south-west  to 
north-east,  which  conforms  with  the  boundaries  of  physiogra- 
phic zonality.  Therefore,  in  remote  determination  of  the 
water  equivalent  of  snow  cover  using  formula  (3),  zonation 
of  the  value  k   is  taken  into  account. 

3.  Parameterization  of  the  dynamic  member  X  \Jj) 
of  equation  (1)  is  aimed  at  reconstruction  of  precipitation 
fields  in  the  standard  grid  square  (2.5°  x  2.5°)  from  cloud 
fields  fixed  in  TV  images.  The  calculation  is  made  for  the 
total  daily,  10-day  and  monthly  precipitation  using  formulae 
(4)  and  (5)  respectively 
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A  r*m  ^~    rj    j      j 


W 


X-^-^Pi1!3^ 


(5) 


Here  Ni  ■  denotes  the  indices  of  groups  of  large-scale  and 
mesoscale  cloud  fields;  S ;-  ,  rj   and  i~  ;    are  the 
areas  of  cloud  cover,  probability  and  intensity  of  precipi- 
tation respectively,  for  each  group.  The  selection  of  cloud 
groups  and  calculation  of  probabilities  and  mean  intensity 
of  their  precipitation  were  made  using  prior  sampling  in 
the  analysis  of  both  the  satellite  information  on  clouds  and 
surface  rainfall  data.  The  basic  amount  of  precipitation 
is  carried  by  the  frontal  cloudiness  -  though  the  amount 
of  cyclonic  clouds  for  10  days  is  two  times  less  than  that 
of  air-mass  clouds,  the  former  carries  two-fold  more  precipi- 
tation than  the  latter.  The  averaged  and  reduced  rainfall 
rates  in  different  groups  (Table  1)  show  that  the  air-mass 
cloud  contribution  to  precipitation  is  almost  an  order  less 
than  that  of  frontal  cloudiness.  Therefore,  the  intensity 
of  the  former  may  be  given  by  a  constant.  In  case  of  large- 
scale  cloudiness,  different  phases  of  the  development  of 
the  cyclonic  cloud  system  are  identified,  and  in  compliance 
with  them  the  groups  of  large-scale  cloudiness  have  been 
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assumed  (Table  1). 

4.  Runoff  ylil 

The  methods  of  runoff  parameterization  are  different 
for  different  seasons  depending  on  characteric  properties  of 
its  formation.  For  instance,  low  flow,  summer  flow  and  win- 
ter flow  may  "be  neglected  since  they  are  defined  "by  deep 
groundwater  formed  during  the  hydrological  cycle  of  the  pre- 
ceding years.  Over  the  entire  continents  (including  those 
with  sparse  data),  the  estimates  of  spring  flow,  the  contri- 
bution of  which  to  the  hydrological  cycle  is  great,  may  be 
made  by  the  following  two  satellite  methods,,  The  first 
method  is  based  on  the  relationship  between  the  flow  at  the 
outlet  and  the  portion  of  the  catchment  area  where  the  snow 
had  melted#  Such  relationships  in  relative  units  (runoff 
in  per  cent  of  the  maximum  value;   snow  coverage  in  per  cent 
of  the  entire  catchment  area)  were  derived  by  the  authors 
earlier  for  22  catchments  of  the  Don,  Volga  and  Ob.  Fig.  4- 
confirms  such  a  relationship  for  the  catchment  of  the  Ob. 
A  relative  hydrograph  may  be  transformed  into  an  absolute 
hydrograph  via  the  maximum  flow  determined  from  the  nomogram 
of  its  dependency  on  the  water  equivalent  of  snow  cover  We 

and  soil  moisture  storage  in  the  end  of  autumn  WW . 

The  other  method  of  streamflow  routing  uses  the  gene 
tic  formula  of  flow 

v  -    J  vj  (t)  FC?)\cCT)P(t  -<c)dr       (6) 
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n 


where  apart  from  the  areas  where  snow  had  melted  away  r  C*Cj 
and  the  snowmelt  rate  \(  it)  =  k  &r   f  one  more  characte- 
ristics is  used,  namely  the  runoff  coefficient  1)     ,  which 
may  also  be  determined  with  the  aid  of  satellite  data.  The 

y 

authors  made  calculations  of  the  runoff  coefficient  Kj  =.  ■— 

*   vvch. 

as  a  relation  between  the  depth  of  snowmelt  runoff  and  water 
equivalent  of  snow  cover  estimated  by  the  above  diagnostic 
method  V\zcv».~^  Lkk~J    •  A  map  of  the  runoff  coefficient 
was  made  for  a  part  of  the  European  territory  of  the  country 
(Fig.  5c)  on  the  basis  of  preliminarily  prepared  maps  of 
the  depth  of  snowmelt  runoff  for  april  and  may  1972  (Fig#  5"b)» 
and  the  water  equivalent  of  snow  cover  on  5  March  1972 
(Fig.  5a) o  The  calculated  map  of  the  runoff  coefficient 
was  then  compared  with  the  map  of  the  depth  of  autumn  runoff 
for  September,  October  and  November  1971  (Fig.  5d),  which 
was  assumed  as  an  idex  of  the  degree  of  autumn  moistening 
since  thellatter  usually  defines  the  runoff  coefficient  to 
a  great  extent.  The  coefficient  of  correlation  of  the  value 
'f)      with  the  value  of  the  depth  of  autumn  runoff,  calcula- 
ted from  the  values  of  ))       in  the  standard  grid  point  taking 
into  account  two  soil  types  (deep  chernozem  and  leached 
podzolic  chernozem)  turned  out  to  be  equal  to  0#9»  Without 
taking  into  account  the  soil  types,  the  correlation  coeffi- 
cient was  substantially  lower. 
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5«  Estimation  of  evaporation  from  large  territories 
for  short  intervals ♦  Analysis  of  soil  moisture 

movement 

Evaporation  for  short  intervals  is  -usually  estimated 
for  a  point,  as  the  residual  member  of  the  water  "balance 
equation.  Maps  of  climatic  values  of  evaporation  for  large 
territories  were  prepared  using  formal  interpolation  of 
point  data*  To  estimate  evaporation  for  short  intervals 
from  large  territories,  the  authors  performed  modification 
of  Oldecop  analytical  relationship  derived  for  conditions 
of  climatic  averaging 

where  t  "~\  (  K  cl  IX  X  °)  ^s  determined  from  Penman- 

Budagovsky  formula  through  the  radiational  balance  R,  »  - 
moisture  deficit  a      >  air  temperature  t°  and  wind  velocity 
li  ;  X   is  precipitation;  AW  is  the  change  of  stored 
soil  moisture  in  the  onenneter  layer. 

Analysis  of  soil  moisture  movement  at  different 
levels  up  to  1  m,  together  with  determination  of  evaporation 
from  formula  (7)i  has  shown  that  the  argument  of  equation 
(7)  may  "be  interpreted  as  a  parameter  as  a  parameter  classi- 
fying the  relation  "between  heat  and  moisture  and  determining 
"both  the  trend  of  the  water  exchange  processes  in  the 
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evaporation  system  and  the  types  of  evaporation  process. 
Algorithms  of  the  computation  of  evaporation  from  formula 
(7)  is  presented  in  Table  2. 

At  ft  <'  0.5,  the  conditions  of  insufficient  moiste- 
ning are  realized.  Under  these  conditions  the  evaporation 
process  is  described  "by  a  linear  relationships  (Fig.  6)# 

At  0.5  ^  ft  sT  ftfcp,/.™.   i  the  conditions  of  optimal 

moistening  are  realized;  the  nonlinear  portion  of  Oldecop 
curve  complies  with  these  conditions.  And  when  the  value 
approximates  ifV  fc  ^H?~m.  ,  the  deviations  from  linearity 

increase  sharply  (Fig.  6). 

At   KL  2    Kjcp-wi  the  process  tends  to  saturation,  a 
attaining  the  maximum  possible  magnitude  E  =  E  which  results 
from  the  small  radiation  balance.  Thus,  at  'ft  >   n  £■,*»,'<,  , 
the  conditions  of  excessive  moistening,  and  evaporation  is 
governed  only  by  the  complex  of  atmospheric  conditions. 

The  values  of  parameter  ft,/  also  fix  different  move- 
ments of  soil  moisture  (Table  3)»  Since  at  the  insufficient 
moistening,  the  decrease  of  the  value  r'V  from  0.5  to  0 
(and  even  to  negative  values)  is  in  line  with  the  decrease 
of  the  downward  moisture  flux  from  the  one-meter  layer 
practically  to  zero  and  even  to  the  occurrence  of  upward 
moisture  fluxes  from  the  lower  layers,  at  the  increase  of 
the  value  ]\>    from  0.5  up  to  fVoam,  >  the  increase  of 

downward  moisture  flux  is  observed  as  compared  with  the  case 
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of  ft,  ^  0.5  which  is  mainly  defined,  by  an  increase  of  preci- 
pitation (Table  3)» 

At  the  excessive  moistening  ft,  -;.35'  rV^pi..^.  ,  the 
saturation  of  the  one-meter  layer  and  downward  moisture 
flux  to  the  lower  layers  takes  place.   Only  the  late  autumn 
period  with  small  radiation  balance  and  abundant  antecedent 
moistening,  responsible  for  the  saturation  of  the  one-meter 
layer  through  the  soil  lens,  presents  an  exception  (Table  3)» 

Computation  of  evaporation  from  formula  (7)  is  liable 
to  have  results  similar  to  those  obtained  using  the  heat 
balance  method.  And  the  relative  error  for  the  conditions 
of  optimal  and  excessive  moistening  is  not  more  than  5  per 
cent.  The  difference  amounts  to  25  per  cent  only  under  the 
conditions  of  insufficient  moistening  for  which  the  heat 
balance  method  proves  inaccurate.  Maps  of  the  parameter  KU 
(Fig.  7)  provide  the  patial  zoning  according  to  the  types 

of  evaporation  processes  and  soil  moisture  flux. 

Q  f ' it  \ 

Another  parameter  -  Budyko  radiation  index  — -  ~'j(X) 

La 

also  takes  into  account  (though  more  roughly)  the  relation 
between  heat  and  water  resources.  This  proves  the  simila- 
rity  of  maps  of  parameters  YY,    and  —    (Figs.  7  and  8). 

Good  agreement  between  the  isolines  of  parameters  ru  and 

R 

7—;     with  the  precipitation  isolines  shows  that  basic 

deviations  from  the  zonality  in  the  character  of  moistening 

are  defined  by  large-scale  cyclonic  cloud  systems.  There- 

fore  parameter   - —    ,  the  numerator  and  denominator  of 

Ux 
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which  may  be  expressed  by  the  cloud  field,  can  be  determined 
only  using  satellite  TV  images.  Then  taking  in  account  the 
similarity  in  the  fields  of  parameters  ft   and  ~[~^   ,  i.e. 


»K(&) 


assuming  that  \\-  4  y^z)   come  to  the  principle  of  satellite 

X  —  A  Vv 
assessment  of  the  fields  of  the  parameter  fL-  • — p — ~~  * 

b  o 

and  consequently  to  the  background  assessment  of  soil  mois- 
ture fields  V\/  ;  and  using  formula  (7)  we  come  to  the  back- 
ground assessment  of  evaporation. 

Conclusion 

The  proposed  parameterization  (1)  describing  the 
dynamics  of  the  annual  continental  hydrological  cycle  was 
studied  as  applied  to  the  zones  of  snowmelt  feeding.  The 
zone  under  consideration  is  in  moderate  latitudes.  It  is 
a  zone  of  interaction  between  the  equatorial  and  polar 
regions,  and  is  characterized  by  a  multiphase  hydrological 
systems  Therefore  the  hydrological  cycle  is  formed  here 
by  the  natural  moisture  integrators  (  W,v  and  VW^  )  during 
sharply  distinguished  genetic  periods.  These  integrators 
define  the  prehistory  of  the  cycle  and  become  the  initial 
conditions  for  the  hydrological  cycle  of  the  succeeding 
year,  thus  relating  the  cycles  of  adjacent  years.  Dynamic 
parameters  of  equation  (1)  characterize  the  governing  in- 
fluence of  the  atmosphere  upon  and  the  role  of  soil  moisture 
in  the  formation  of  the  cycle  and  should  be  computed  inde- 
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pendently.  Physical  regularities  of  the  moisture  exchange 
to  be  found  out  in  a  similar  parameterization  will  prove 
essential  for  the  study  of  year-to-year  variations  of  the 
water  exchange. 
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Table  1 

Jrobabilit 

;ies  and  intens 

iity  of  precipitation  for 

9  groups 

of  clouds  obtained  from  a 

two-month  sample 

Groups  of 
clouds 

Probability 
of  precipi- 
tation in 
the  portions 
of  a  unit 

Intensity 
of  preci- 
pitation 
i  in  mmd 

Intensity  of 
precipitation 
in  mmd  per  1 
cloud  ball 

Cyclonic 

formations 

1.  Wave 

1.0 

12.3 

1.0 

2.  Young  cyclone 

1.0 

16.7 

1.5  -  2.1 

3.  Front:  - 

(a)  continuous 
cloudiness 

1.0 

7.8 

0.8  -  1.3 

(b)  broken 

cloudiness 

0.9 

3.7 

0.2  -  1.0 

(c)  the  ends  of 
the  front  are 
broken 

OoO 

0.0 

mm 

4.  Point  of 
occlusion 

1.0 

10.3 

0.9  -  1.2 

Air-mass  clouds 

5.  Cumulonimbus 
clouds 

0.8 

0.22 

0.16 

6©  Stratiform 
clouds 

0,6 

0.1 

0.01 

7.  Cumuliform 

0.4 

0.1 

0.01 

8«  Strato-cumuliform 

clouds             - 

- 

~ 

9.  Cirrus 

0.1 

0.01 

- 
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Table     2 


Algorithm  of  computation  of  evaporation 


■h 


Value 
n 


Value  R 
Kcal 


cm  month 


Portion  of 

Oldecop 

curve 


Formula  of 
calculation 
of  evapora- 
tion 


Formula  of 
calculation  of 
upward  and  down- 
ward moisture 
fluxes 


LOihiO^ 


any 


2.    (a) 


(b) 


any 


i^<ft^ 


5.  n  £n 


t<j>um 


§•  r\  >  n  up 


5.  n  <o 


R>0.5 


R<0.5 


6.    n  <  o 


R<0.5 


R>0.5 


linear 


l-trr*    ~&W 


X-AvV 


v  0   (within  the 
range  of  the 
measurement  error) 


weak  non- 
linearity 


strong  non- 
linearity 


strong  non- 
linearity 


plateau 


plateau 


linear 


T>  ^  x-aW-E 


E^E0ttx^ 


D  =l  X  -  Aw  -  b 


t=  Eo 


D  -  X  -  AW  -  E0 


E»  E 


IX  s-*2>~A\AMEo-* 


_       _     j  i    X+  AW 

t-t0"t«  T~ 


*U*-E>«AW+E-X 
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Table  3 


Examples  of  types  of  evaporation  and  soil  moisture 


R 


movement 

1 

2(a) 

2(b) 

3 

4 

5 

6 

ft> 

0.26 

0.80 

1.34 

2.22 

2.63 

-1.22 

-0.12 

Kcal 

3«6 

1.1 

1.1 

3.3 

0.2 

0 

1.3 

2 

cm  month. 

X   mm 

12.1 

37.3 

37.0 

176.6 

42.3 

39.5 

5.5 

A  \A/  mm 

-9 

17 

6 

86 

8 

51 

10 

^-^mm 

21.1 

20.3 

31.0 

90.6 

34.3 

-11.5 

-4.5 

lu  _   mm 

79.7 

25.5 

23.2 

40.7 

13.0 

9.4 

38.0 

•—   mm 

20.6 

16.9 

20.2 

39.7 

13.0 

9.4 

14.8 

L*-   mm 

0 

0 

0 

0 

0 

20.9 

19.3 

From  mm 

0.5 

3.4 

10.8 

50.9 

21.3 

0 

0 

Askaniya- 

Krasno- 

Cherno- 

Lubny 

Kherson 

Poltava 

Kherson 

Nova 

armeisk 

byl 

Septem- 

November 

November 

October 

Septem- 

October 

Octo- 

ber 

Winter 

Winter 

Winter 

ber 

Winter 

ber 

Winter 

wheat 

wheat 

wheat 

Fall  plow- 

wheat 

Winter 

wheat 

after 

after 

after 

ed  land 

after 

wheat 

after 

corn 

barley 

corn 

after 

winter 

after 

pea 

winter 

wheat 

flax 

! 

V 

/heat 
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List  of  Figures 
of  the  paper  "Study  of  the  Hydrological  Cycle  Based  on 
Satellite  and  Surface       Information"  by  Yu.V.  Kurilova, 
P.A.  Kolosov,   E.L.  Muzylev,   L.K.   Poplavskaya 

Fig.   1#     Scheme  of  the  hydrological  cycle  for  the  zones  of 

settled  (a)  and  unsettled  (b)   snow  cover.     The  dimen- 
sions of  arrows  corresponding  to  the  components  of 
the  water  balance  are  proportional  to  their  values. 

Fig©   2.     Changes  of  snowmelt  rate   (a)  and  mean  temperature   (b) 
on  the  territory  of  the  Don  and  the  Volga  in  the 

pre-spring  period.     March  1972. 

1 .  Snow  patches  tlne^  sag*/-  izjia  aXsUL  .  2.Ccai-uitou!>  §nov/ 
cov/et^  -Lhe*v  skvow  bo/tches.  3  Con.tln/.ouii  .snow  cove^j 
then,  SnoW-fuzz   a^ea.  .  4.  J"he  7)on .    5.  P^o/CaS/Ow 

Fig.  3«  Distribution  of  the  coefficient    (formula  (3)) 

based  on  the  series  of  the  satellite  snow  cover  maps 
(March  1972)  for  catchment  areas  20000-60000  km2. 

Fig.  4.  Relationship  between  the  portion  of  the  snow-free 
basin  area  F  (in  per  cent)  and  runoff  (in  per  cent 
of  Gl   max)  at  the  outlet  for  catchment  areas  of  the 
Ob  river. 

(1)  the  Ishim;T  =  19  days; 

(2)  the  Tobol,  fa  15  days 

(3)  the  Iset,  'XT  a  2  days 

(4)  the  Tura,  "tT  *  9  days 
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Pigo  5»  Maps  illustrated  the  calculation  of  the  coefficient 
of  the  snowmelt  depth: 

(a)  the  water  equivalent  of  snow  cover  vych  on  5 
March  1972;   (b)  snowmelt  depth  for  April,  May  1972; 
(c)  coefficient  of  the  runoff  depth I J  ;   (d)  depth 
of  the  autumn  runoff  for  September,  October,  November 
1971. 

Fig.  6.  Set  of  Oldecop  curves  for  different  months  and  the 
different  conditions  of  moistening. 

Fig.  7«  Maps  of  isolines  of  monthly  total  precipitation 

and  values  of  the  parameter  tl   • 

o 
Fig.  8.  Maps  of  isolines  of  the  Budyko  radiational  index  j^- 
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The  National  Weather  Service 
Extended  Streamflow  Prediction  Techniques: 
Description  and  Applications  During  1977 

Dr.  John  C.  Schaake,  Jr.,  NOAA,  MS,  W22 
Silver  Spring,  Maryland  20910 

Streamflow  forecasts  are  needed  for  a  variety  of  water  resource  management 
purposes.  They  are  needed  for  future  periods  ranging  from  a  few  days  to 
several  months.  As  the  duration  of  the  forecast  period  increases  so  also 
does  the  uncertainty  in  the  forecast-.  Recognizing  this  uncertainty,  the 
National  Weather  Service  has  developed  an  extended  streamflow  prediction 
model  (ESP)  to  produce  probabilistic  estimates  of  future  streamflow. 

The  ESP  forecast  model  includes  within  it  a  conceptual  hydrologic  model 
representing  the  movement  of  water  through  the  hydrologic  cycle.  Inputs  to 
the  hydrologic  model  are  climatic  variables:  precipitation,  temperature, 
and  potential  evapotranspiration.  The  hydrologic  model  responds  to  these 
climatic  inputs  by  accumulating  and  ablating  snow  cover  during  the  winter 
months  by  simulating  the  movement  of  water  through  soil  and  ground  water 
systems  and  by  representing  the  flow  of  water  in  the  surface,  stream,  and 
river  systems.  This  hydrologic  model  must  be  calibrated  for  a  specific 
location  on  the  basis  of  historical  climatic  and  streamflow  data.  Model 
calibration  is  needed  because  model  parameters  vary  from  place  to  place  and 
can  only  be  estimated  on  the  basis  of  historical  data. 

To  use  the  conceptual  model  in  streamflow  forecasting,  the  initial  condition 
of  the  hydrologic  system  must  be  known  and  forecasts  of  the  climatic  vari- 
ables are  needed  as  well.  Initial  conditions  can  be  known  by  operating  the 
model  continuously  in  time  with  observed  data,  using  streamflow  data  to 
update  model  state  variables  so  that  the  model  is  in  tune  with  the  most 
current  data. 

Since  ESP  uses  a  simulation  approach,  specific  time  series  of  future  input 
climatic  variables  are  needed.  Since  the  future  input  series  are  unknown, 
it  becomes  necessary  to  work  with  a  sample  of  equally  likely  input  series. 
The  ESP  model  must  be  operated  starting  from  the  same  current  initial  condi- 
tion for  each  of  the  input  series  in  the  sample  and  produces  a  streamflow 
output  for  each  input  series.  These  streamflow  output  series  can  then  be 
used  to  derive  the  forecast  probability  distribution  for  the  streamflow 
variables  of  interest. 

If  there  were  no  skill  in  forecasts  of  precipitation,  temperature,  or  evapo- 
transpiration, an  appropriate  sample  of  climatic  input  time  series  would  be 
offered  by  the  historical  record.  In  other  words,  historical  series  of  pre- 
cipitation, temperature,  and  evapotranspiration  starting  from  the  current  day 
in  each  of  the  past  years  could  be  used  as  input  to  the  hydrologic  model. 
For  each  of  these  years  the  model  would  simulate  the  streamflow  that  would 
have  occurred  in  those  years  if  the  current  hydrologic  conditions  were  to 
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have  existed  at  that  time.  In  the  limit  as  the  forecast  period  becomes  long, 
the  simulated  streamflow  will  approach  the  observed  streamflow  because  the 
effect  of  the  initial  conditions  dissipates  with  time.  It  follows  then  that 
very  long  term  streamflow  forecasts  from  the  ESP  model  will  be  the  same  as 
given  by  the  historical  streamflow  records  (within  the  limits  of  model  error). 

This  approach  was  taken  to  forecast  the  lowest  flows  on  record  in  the  state 
of  California  during  1977.  As  an  example,  the  American  River  natural  inflow 
to  Folsom  Dam  near  Sacramento,  California,  for  the  month  of  May  1977,  was 
predicted  to  be  between  45»000  and  90,000  acre-feet  with  the  most  probable 
value  of  60,000  acre-feet.  The  historical  median  flow  for  May  is  540,000 
acre-feet.  This  forecast  was  made  on  April  15,  1977.  The  main  factor  domi- 
nating this  particular  forecast  was  the  initial  condition  of  snow-water 
equivalent  in  the  Sierra  Nevada  Mountains  and  depleted  soil  moisture  as 
computed  by  the  hydrologic  model. 

Following  the  western  drought  of  1977  and  the  application  of  ESP  in  California, 

another  drought  occurred  in  the  east  • but  on  a  much  more  limited  spacial 

scale.  This  occurred  in  northern  Virginia  and  the  water  supply  for  parts  of 
the  metropolitan  Washington,  D.C.,  area  were  severely  threatened.  This  was 
the  driest  year  in  the  last  27  years.  The  ESP  model  was  used  to  assess  the 
situation  and  to  analyze  alternative  water  conservation  measures  that  would 
bring  the  risk  within  acceptable  levels.  The  ESP  application  began  late  in 
September  of  1977. 

One  of  the  significant  problems  in  the  northern  Virginia  case  was  the  tremen- 
dous uncertainty  in  the  future  quantities  of  streamflow  to  be  expected.  The 
Occoquan  Reservoir  serving  northern  Virginia  was  almost  empty.  Inflow  to  the 
reservoir  from  natural  ground  water  seepage  was  less  than  the  rate  water  was 
being  taken  from  the  reservoir.  Before  runoff  could  occur  from  the  land 
surface  it  was  necessary  for  the  soil  moisture  and  ground  water  systems  to 
be  replenished  by  rainfall .  At  the  same  time,  natural  evapotranspiration  was 
continuing  to  deplete  the  water  remaining  in  the  soil  moisture  and  ground 
water  systems.  The  soil  moisture  and  ground  water  deficiencies  totaled  about 
5.5  inches.  This  was  of  considerable  concern  because  the  mean  October  rain- 
fall of  2.83  inches  is  rnuch  less  than  this  moisture  deficit.  Preliminary 
studies  had  shown  that  the  historical  streamflow  during  the  driest  years  in 
the  past  when  combined  with  the  depleted  contents  of  the  Occoquan  Reservoir 
would  not  suffice  to  meet  the  needs  for  water  at  the  current  demand  rate. 

Applying  the  ESP  technique  to  the  northern  Virginia  case,  it  was  important 
to  consider  if  the  future  precipitation  series  were  adequately  represented 
by  the  historical  series  or  if  there  were  some  information  available  in  this 
particular  year  that  would  indicate  the  forthcoming  October  would  be  wetter 
or  drier  than  normal.  Therefore,  we  looked  very  carefully  at  the  forecast 
information  available  from  the  National  Weather  Service  long  range  prediction 
unit  for  October.  Don  Gilman  did  not  believe  that  a  very  accurate  quantitative 
precipitation  forecast  could  be  made  for  the  month  of  October,  but  he  did 
believe  that  the  best  possible  way  to  use  the  available  data  was  through  an 
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anti-analog  approach.  According  to  this  approach  it  would  be  possible  to 
judge  that  certain  of  the  historical  years  were  not  representative  of  the 
current  year  because  the  forecast  of  upper  air  circulation  pattern  for 
October  was  vastly  different,  in  a  sense  opposite,  from  the  historical  upper 
air  circulation  patterns  for  October  in  a  number  of  years.  A  total  of  5 
years  were  excluded  from  the  historical  record  and  the  remaining  record  was 
judged  to  be  representative  of  the  type  of  precipitation  pattern  that  might 
be  expected  in  the  current  year.  The  mean  of  this  conditional  distribution 
was  only  about  a  tenth  of  an  inch  above  the  mean  for  the  climatic  distribution 
for  October  and  the  standard  deviation  of  the  condition  of  distribution  was 
about  7$  greater  than  the  climatic  distribution. 

The  ESP  technique  was  applied  in  northern  Virginia,  using  the  complete  his- 
torical precipitation  series.  It  was  found  there  is  about  a  10$  chance 
that  reservoir  levels  in  the  Occoquan  Reservoir  would  fall  below  the 
threshold  that  would  require  the  water  authorities  to  declare  a  so-called 
stage  three  emergency  which  would  prohibit  all  water  use  not  essential  for 
maintaining  the  health  and  safety  of  the  community.  Conversely  there  was  a 
90$  chance  of  avoiding  a  stage  three  emergency  and  there  was  a  considerable 
chance  that  the  reservoir  would  refill  before  the  end  of  the  1977  calendar 
year.  But  there  was  a  10$  risk  of  reaching  a  stage  three  emergency.  This 
was  judged  by  the  community  to  be  unacceptable.  Therefore,  additional 
studies  were  made  with  the  ESP  technique  to  determine  what  conservation 
measures  were  required  to  reduce  the  risk  to  a  more  acceptable  level.  Although 
political  representatives  of  the  community  were  unable  to  define  an  unaccept- 
able risk  level,  it  was  judged  that  a  20$  reduction  in  withdrawal  from  the 
reservoir  would  reduce  the  risk  level  to  about  2  or  3  percent.  Plans  were 
made,  therefore,  to  reduce  withdrawals  until  the  reservoir  had  recovered  to 
a  normal  level.  Heavy  rains  began  to  occur  in  October  and  lasted  through 
November.  By  mid-November  the  reservoir  had  refilled  and  the  emergency  was 
over.  Analysis  of  the  precipitation  for  October  and  November  showed  that 
there  was  only  a  five  percent  chance  of  having  such  a  wet  fall.  Northern 
Virginia  was  indeed  very  lucky. 
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ISSUE  DATE       April,  15,    1977 

EXTENDED  STREAMFLOW  PREDICTIONS 


The  extremely  dry  winter  throughout  California  is  expected  to  result  in 
the  lowest  flows  of  record  through  most  of  the  state.  Actual  flows  will 
be  primarily  dependent  upon  the  variation  from  normal  temperatures  and 
precipitation. 
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FORECASTS  


BASIN 

FORECAST 
PERIOD 

MOST 
PROBABLE 

REASOh 
MAXIMUM 

IABLE 
MINIMUM 

HISTORIC 
MEDIAN 

1  CLEAR  LAKE  INFLOW 

1  Includes  evaporation 

i  affects 

1 

i 

j 

April 

-  -6 

0 

-7      \ 

May 

-18 

-8 

-21 

June 

-25 

-23 

-28    j 

Julv 

-32 

-25 

-34 

\  / 

Auk 

-31 

-25 

-34 

A 

Sept 

-21 

-14 

-24 

/    \ 

Oct 

-12 

1 

-15 

/       \ 

Nov 

0 

18 

-8 

v           \ 

!  TRINITY  RIVER 

Inflow  to 
J  TRINITY  DAM 

I 

April 

35 

50 

28 

210 

Mav 

20 

35 

14 

230 

June 

5 

20 

4 

100 

JulV 

2 

10 

1 

28 

Aug. 

2 

8 

1 

10 

f 

SeDt 

2 

8 

1 

8 

Oct 

5 

15 

3 

13 

Nov 

6 

Lx     30 

4 

f    29 

1 

April 

200 

270 

180 

615 

!  SACRAMENTO  RIVER 

Mav 

175 

220 

160 

462 

j   Natural  Inflow  to 

June 

160 

200 

145 

1   292 

j  SHASTA  DAM 

July 

155 

195 

140 

|   224 

■ 

Aue 

150 

180 

140 

206 

Sept 

145 

175 

135 

200 

Oct 

155 

210 

135 

!   227 

1 

Nov 

160 

1   280 

140 

274     | 

! 
1 

April 

95 

125 

75 

I    720 

j  FEATHER  RIVER 

May 

80 

125 

70 

1   580     1 

]   Natural  Inflow  to 

June 

50 

100 

40 

I   290 

i  -OROVILLE  DAM 

Julv 

45 

90 

35 

135 

i 

Aug 

45 

90 

35 

98 

1 

Sept 

45 

95 

35 

i    86 

i 

t 

Oct 

55 

100 

35 

t    89 

! 

Nov 

65 

i   150       35 

1   126 
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FOREC 

ij 1j  .  .  .  . 

BASIN 

FORECAST 
PERIOD 

MOST 
PR03ABLE 

REASO 
MAXIMUM 

NA3LE 
MINIMUM 

HISTORIC 

MEDIAN 

FEATHER  RIVER 
North  Fork 

Inflow  to 
LAKE  ALMANOR 

April 

25 

40 

20 

100 

May 

25 

40 

20 

105 

June 

20 

35 

15 

63 

Jn  1  v 

20 

35 

15 

LZ 

Aug 

20 

35 

15 

42     j 

Seot 

20 

40 

10 

38     ! 

Oct 

30 

50 

15 

39 

Nov 

35 

60 

20   J 

43 

YUBA  RIVER 
Natural  Inflow 
above 
SMARTVILLE 

April 

45 

120 

30 

4T0" 

Mav 

45 

80 

20 

440 

June 

12 

35 

7 

200 

July 

4 

15 

2 

45 

Aug 

2 

10 

1 

23 

Sept 

3 

20 

3 

20 

Oct 

4 

30 

3 

27 

Nov 

10 

50 

5 

48 

AMERICAN  RIVER 
Natural  Inflow 

to 
FOLSOM  DAM 

April 

70 

100 

55 

480 

May 

60 

90 

45 

5^0 

June 

16 

35 

10 

280 

July 

2 

15 

1 

59 

Aug 

1 

10 

1 

17 

Sept 

1 

20 

1 

12 

Oct 

4 

30 

3 

19 

Nov 

8 

50 

6 

44    i 

AMERICAN  RIVER 
Middle  Fork 

Natural  Inflow 

near 
AUBURN 

April 

34 

50 

27 

205 

Mav 

23 

34 

17 

226 

June 

4 

10 

3 

105 

July 

1 

7 

0 

24 

Aug 

0 

5 

0 

7 

Sept 

0 

7 

0 

5 

Oct 

2 

15 

1 

6 

Nov 

4 

25 

2 

15     I 

COSUMNES  RIVER 

above 
MICHIGAN  BAR 

April 

4 

9 

.    3 

60     ] 

.  Mav 

2 

6 

1 

39     i 

June 

0 

2 

0 

15 

July 

0 

1 

0 

2.6 

Aue 

0 

0 

0 

0.7 

Sept 

0 

0 

o 

0.4 

Oct 

1 

2 

o 

1.1- 

Nov 

3 

5 

0 

3.4 
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Figure  10. — Position  Analysis:   Simulated  storages  for  the  1930-1931 
drought,  at  production  rates  of  32,  40,  and  48  Mgal/d. 
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Figure  13. — Cumulative  distribution  function  for  S _.  $ 
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Analysis  based  on  U9  years  of  record. 
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EXCEEDANCE  PROBABILITY  FOR  MINIMUM  STORAGE,  IN  PERCENT 

Figure  lU. —Cumulative  distribution  function  for  S .    ,  Uo  Mgal/d 
production  rate,  from  Position  Analysis  based 
on  22  years  with  dry  late-summer  periods. 
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!   FIGURE  2 


TABLE  1 

PROBABILITIES  OF  FALLING  BELOW  STORAGE  LEVELS 
AT  40  MGD  (CURRENT  PRODUCTION  RATE)  AT  OCCOQUAN* 


Storage  Level 

93  feet 

88  feet 

83  feet 

Volume  of  Storaqe  Remaining 

1.7  B.G. 

1.1  B.G. 

0.7  B.G 

Probability 

13% 

10% 

8% 

*  Purchase  of  1.5  B.G.  from  Lake  Manassas  is  assumed. 
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TABLE  2 

PROBABILITY  OF  FALLING  BELOW  STORAGE  LEVELS 
AT  32  MGD  PRODUCTION  AT  OCCOQUAN* 


Storage  Level 

93  feet 

88  feet 

83  feet 

Volume  of  Storage  Remaining 

1.7  B.G. 

1.1  B.G. 

0.7  B.G 

Probability 

5% 

3% 

Z% 

*  Purchase  of  1.5  B.G.  from  Lake  Manassas  is  assumed. 
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APPLICATION  OF  THE 
UNIVERSITY  OF  WISCONSIN 
MODEL  TO  CLIMATE  PREDICTION 

Gerald  J.  Dittbemer 

USAF  Environmental  Technical  Applications  Center 

Scott  AFB,  IL  62225 


Studies  of  past  climates  indicate  that  there  is  no  single  cause  of 
climatic  variation.  Evolution  of  the  atmosphere,  variations  of  the 
earth's  orbital  parameters,  solar  variations,  continental  drift,  vol- 
canic activity  variations,  and  still  other  causes  are  all  thought  to 
be  important  at  different  time  scales.  For  scales  such  as  the  coming 
season  and  the  coming  year  many  of  these  potential  causes  of  variation 
can  be  considered  negligible. 

The  climate  prediction  model  was  constructed  based  on  two  sub- 
models. The  low  frequency  sub-model  (Bryson  and  Dittbemer,  1976  and 
Dittbemer,  1978)  deals  with  trends  and  differences  in  decadal  or 
pentadal  means  assuming  until  otherwise  demonstrated  that  these  vari- 
ations are  determined  by  extrinsic  physical  factors  such  as  changes  in 
volcanic  activity,  anthropogenic  aerosols  and  carbon  dioxide  and  as- 
suming that  there  is  no  inherent  random  variation.  Unexplained  vari- 
ance of  the  observed  climate  from  the  model  was  regarded  as  a  short- 
coming of  the  model  or  of  the  observations  and  not  as  a  fundamentally 
unexplainable  random  variation  in  climate. 

In  the  high  frequency  sub-model  (Bryson  and  Starr,  1977)  it  was 
assumed  that  part  of  the  interannual  variation  of  monthly  mean  pat- 
terns of  pressure,  temperature  and  precipitation  is  associated  with 
the  complex  set  of  atmospheric  tidal  oscillations  related  to  the 
Chandler  Motion  of  the  earth's  pole.  Model  construction  uses  regres- 
sion to  link  these  input  factors  to  monthly  mean  temperature  and 
monthly  total  precipitation  at  21  stations  in  the  United  States. 

The  importance  of  estimating  effectiveness  in  climate  prediction 
methods  cannot  be  overemphasized.  Accordingly,  a  technique  was 
developed,  based  primarily  on  that  described  in  Barnett  and  Preisen- 
dorfer  (1977)  to  estimate  the  effectiveness  of  the  University  of  Wis- 
consin Model.  For  each  station,  regression  coefficients  were  computed 
on  a  dependent  data  set  (1882  -  1950)  and  used  to  calculate  values  for 
1951  through  1955.  The  procedure  was  then  repeated  using  data  from 
1883  -  1951  to  calcuate  values  for  1952  through  1956.  Repeating  this 
procedure  through  1974  produced  23  one  year  ahead  calculations,  22 
2-year  ahead  calculations,  etc.  Statistical  significance  was  esti- 
mated for  a  three  category  (H,M,L)  system  (each  category  having  a 
probability  of  one  third)  using  the  binomial  distribution.  Values  of 
significance  were  determined  as  in  the  following  example.  For  the 
case  where  there  is  a  total  of  11  correct  categories  out  of  24  at- 
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tempts  we  find  that  the  probability  of  11  or  more  correct  out  of  24  is 
0.14  (for  a  probability  of  occurrence  by  chance  of  1/3).  Subtracting 
this  from  unity  gives  a  value  of  0.86.  This  can  be  considered  a  mea- 
sure of  forecast  significance,  skill  or  confidence. 

Analysis  of  the  confidence  values  for  21  stations  in  the  United 
States  indicates  the  following: 

1 .  There  are  indicatims  that  the  confidence  or  skill  of  this 
model  may  not  be  primarily  a  function  of  forecast  lead  time. 

2.  This  method  may  not  be  applicable  to  all  stations  in  all 
months;  however  it  is  possible  that  for  some  stations  in  some  months 
the  method  may  demonstrate  significant  success. 

3.  The  geographic  distribution  of  calculated  confidence 
values  suggests  that  the  method  may  be  applicable  to  stations  in  the 
Dakotas  for  January  mean  temperature,  the  eastern  United  States  for 
April  mean  temperatures,  the  southern  and  eastern  United  States  for 
January  total  precipitation  and  Texas  for  April  total  monthly  pre- 
cipitation. 

Intriguing  relationships  between  the  mapped  confidence  patterns 
and  potentially  related  factors  (e.g.,  airmass  dominance,  importance 
of  local  effects,  possibilities  for  improving  the  forecast  methods, 
etc.)  provide  fruitful  ground  for  continued  study 
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PREDICTABILITY  OF  AIR  TEMPERATURE  OVER 
NORTH  AMERICA  -  A  PROGRESS  REPORT 


T.  P.  Barnett 


Introduction 


It  has  often  been  stated  that  properties  of  the  ocean  and  atmospheric  field  in 
the  Pacific  Ocean  can  be  used  to  predict  atmospheric  variables  over  North  America. 
This  is  a  progress  report  on  a  quantitative  attempt  to  verify  such  hypotheses. 

The  method  of  linear  prediction  has  been  used  to  test  a  variety  of  hypotheses. 
The  details  of  the  method,  i.e.,  mathematics  significance  test,  etc.,  are  given  in 
Barnett  and  Hasselmann  (in  press)  and  will  not  be  repeated  here.  Suffice  it  to  say, 
the  technique  attempts  to  take  a  series  of  specified  predictors  and  see  how  well  they 
can  be  used  to  forecast  a  given  variable  at  different  lead  times.  The  crucial  new 
element  in  the  approach  has  been  to  balance  predictive  skill  with  statistical 
significance.  Thus  it  is  possible  to  make  models  with  higher  skill  than  those  shown 
here.  However,  they  will  surely  have  less  significance.  In  the  following  work  we 
attempt  to  construct  models  that  are  significantly  different  from  the  null -hypothesis 
type  model  at  the  90%  confidence  level. 

The  predictors  used  in  this  study  are  shown  in  Figures  1  and  2.  Areal  averages 

for  the  various  regions  have  been  constructed  from  the  historical  weather  map  series 

■ 

in  the  case  of  sea  level  pressure  and  from  a  recent  updated  version  of  the  Marine  Deck 
(TDF11)  in  the  case  of  sea-surface  temperature.  Thus,  for  each  of  the  index  regions 
shown  in  Figures  1  and  2  we  have  developed  a  time  series  of  either  pressure  or  sea- 
surface  temperature  anomalies  by  season  for  the  period  1899  through  1972.  The  regions 
over  which  the  average  were  computed,  by  the  way,  were  defined  semi-objectively  by 
methods  of  variance  analysis  and  also  on  physical  grounds.  The  time  series  from  the 
index  regions  have  been  compared  where  possible  with  independent  data  sets.  The 
cross  comparisons  have  shown  the  long  time  histories  to  be  surprisingly  good  in  their 

33-1 


quality  and  continuity. 

Specifications 

Before  trying  to  predict  anything,  it  is  first  important  to  ask,  given  a  set 
of  information,  can  one  use  it  to  simultaneously  specify  the  properties  of  the  variable 
that  eventually  we  want  to  predict.  With  this  in  mind,  the  sea  level  pressure  variables 
in  winter  have  been  used  to  specify  the  corresponding  winter  air  temperature  anomalies 
at  42  stations  over  the  North  American  continent.  The  results  of  this  specification 
are  shown  in  Figure  3.  The  numbers  show  the  percent  of  variance  accounted  for  in  the 
specification  (62%  of  the  variance)  leave  a  large  measure  of  variation  unaccounted  for. 
The  principal  feature  of  the  diagram  is  the  region  where  skills  of  less  than  30%  were 
found  in  the  Great  Plains  area.  Based  on  this  we  might  expect  our  ability  to  predict 
in  these  areas  will  be  relatively  small. 

We  have  next  taken  the  sea-surface  temperature  anomaly  data  (Figure  2)  for 
winter  and  attempted  to  specify  the  winter  air  temperature  over  North  America.  The 
results  are  shown  in  Figure  4.  It  is  seen  that  the  large  region  of  low  specification 
ability  is  enhanced  and  deepened  by  using  SST  variables.  Furthermore,  the  percent  of 
variance  accounted  for  throughout  the  continental  United  States  is  substantially 
lower  using  the  sea-surface  temperature.  Thus  it  would  appear  the  sea-surface 
temperature  is  not  nearly  as  skillful  at  specifying  air  temperature  anomalies  as  is 
sea  level  pressure. 

One  Season  Predictability 

At  this  early  stage  in  the  work  we  have  considered  only  the  situation  of 
predicting  winter  air  temperatures  from  fall  sea-surface  temperature  and  sea  level 
pressure  anomalies.  Furthermore,  we  have  only  used  the  limited  set  of  predictors 
shown  in  Figures  1  and  2.  Thus  the  results  of  the  following  predictability  experiment 
should  be  considered  as  very   preliminary  in  nature. 
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Prior  to  computing  the  predictabilities,  however,  it  is  important  to  estimate 
the  amount  of  skill  that  a  prediction  based  solely  on  persistence  would  have.  The 
results  are  shown  in  Figure  5  where  it  is  seen  that  the  amount  of  variance  accounted 
for  by  persistence  is  generally  less  than  20%  and  then  only  significant  at  the  90% 
level  in  the  extreme  northwest  part  of  the  continental  U.S.  The  apparent  skills  of 
12  and  20  in  the  central  and  eastern  U.S.  respectively  are  almost  certainly  chance 
values  (indeed  at  the  90%  level  we  should  expect  at  least  4  stations  to  be  significant 
by  chance). 

We  next  attempted  to  predict  the  air  temperature  over  North  America  during 
winter  from  the  fall  sea  level  pressure  anomalies  of  the  Pacific.  The  results  are 
shown  on  Figure  6.  It  is  seen  that  there  is  extremely  low  predictability  over  most 
of  the  continental  U.S.  A  region  of  increasing  predictability  exists  as  one  moves 
from  the  Dakotas  into  southern  Canada.  Also,  as  indicated  in  the  upper  left  corner 
of  Figure  6,  substantial  predictability  exists  at  Nome  and  Juneau,  Alaska.  Thus 
while  the  results  may  be  encouraging  for  those  interested  in  Canadian  winter  forecasts 
they  cannot  be  considered  encouraging  for  those  interested  in  predicting  air 
temperatures  over  North  America  during  the  winter.  Principal  predictors  responsible 
for  the  skill  shown  in  Figure  6  come  from  the  sub  tropical  ridge,  specifically  regions 
6  and  8;  although  region  7  does  contribute  to  the  predictability  also.  Thus  it 
appears  a  major  part  of  the  skill  shown  in  Figure  6  is  associated  with  tropical 
events  during  the  preceding  fall  and  summer  (the  time  structure  of  the  predictor  field 
is  also  available  through  these  studies  [see  cf  Banrett  and  Hasselmann]  but  will  only 
be  mentioned  briefly  here). 

We  next  used  information  on  fall  sea-surface  temperature  anomalies  to  try  to 
predict  winter  air  temperature  anomalies  over  North  America.  The  results  are  shown 
on  Figure  /.  Again  we  observe  that  over  most  of  the  continental  United  States  there 
is  little  or  no  predictability.  Low  levels  of  predictability  exist  in  the  same  region 
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as  they  did  for  sea  level  pressure,  but  again  are  generally  smaller  in  magnitude 
than  available  from  the  pressure  data.  An  interesting  feature  in  the  variance  field 
shown  in  Figure  7  is  the  apparent  existence  of  significant  predictability  along  the 
west  coast  of  the  United  States.  In  fact,  study  of  predictors  for  this  region 
showed  that  the  skill  was  due  to  fall  water  temperatures  immediately  off  the  coast. 
Since  these  water  temperatures  persist  for  a  period  of  seasons,  we  thus  are  seeing  the 
effect  of  air  blowing  over  warm  or  cold  water  and  affecting  the  stations  immediately 
at  the  coast.  Note  that  the  effect  does  not  penetrate  yery   far  inland  since  both 
Sacramento  and  Spokane  have  no  significant  levels  of  predictability. 

The  principal  predictors  associated  with  the  sea-surface  temperature  field  come 
almost  exclusively  from  regions  4  and  5,  and  to  some  extent  region  6.  The  former 
regions  are  associated  with  tropical  sea-surface  temperature  anomalies  exclusively, 
and  particularly  those  generally  associated  with  El  Nino.  The  importance  of  region  6 
as  a  predictor  and  seems  to  imply  land/sea  contrast  is  of  principal  importance  in 
determining  the  resulting  predictability.  It  was  clear  that  the  often  mentioned 
region  of  SST  variability  in  the  central  Pacific,  i.e.,  SST  1  on  Figure  2,  was  not 
important  as  a  predictor  in  the  present  study. 

Conclusions 

Preliminary  results  suggest  there  is  very  little  skill  in  using  Pacific 

predictors  to  forecast  from  fall  to  winter  air  temperatures  over  continental 

United  States.  Significant  and  substantial  skill,  however,  does  exist  for  forecasting 

this  same  variable  in  the  region  of  southern  Canada.  What  skill  that  was  found  came 

almost  exclusively  from  the  tropical  regions  of  the  Pacific  Ocean,  thus  suggesting 

that  there  must  be  direct  teleconnective  effects  between  these  areas  and  the  high 

latitudes  over  North  America.  The  mid-latitude  regions  of  the  central  Pacific  Ocean 

appear  to  play  a  minor  role  in  the  demonstrated  abilities  to  predict  from  fall  to 

winter. 
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Caveats 

The  above  results  must  be  considered  as  tentative  for  a  variety  of  reasons 
which  are  listed  below. 

1.  There  may  be  other  predictors  that  would  have  given  better  forecast  skill. 
However,  the  predictor  regions  used  are  the  ones  most  often  described  in  the  scientific 
literature  as  having  importance  and  predictive  ability  over  North  America. 

2.  It  may  be  argued  that  the  linear  model  used  here  will  fail  in  strongly 
nonlinear  situations.  However,  work  by  a  number  of  authors  has  shown  that  linear 
prediction  models  work  well  even  for  quasi -nonlinear  systems.  Furthermore,  most 
workers  seem  to  feel  that  quasi -linear  models  of  climate  are  more  than  adequate  to 
represent  our  current  data  set.  Thus  the  linear  model  may  not  be  as  much  of  a 
limitation  as  one  might  imagine. 

3.  The  skill  scoring  method  used  to  estimate  predictability  was  the  amount  of 
variance  accounted  for  by  the  forecast  model.  This  is  a  tough  skill  score  and  one  that 
is  unforgiving  of  large  errors  in  model  hindcast.  A  more  lenient  skill  score  could 
perhaps  have  produced  a  result  that  would  have  appeared  more  pleasing  to  the  eye. 

4.  The  significance  test  used  to  evaluate  the  models  is  a  new  approach  to 
the  significance  problem  (Banrett  and  Hasselmann).  There  are  those  that  might  argue 
that  a  significance  test  over  a  hindcast  is  not  a  true  measure  of  model  validity. 

In  fact  some  might  argue  the  only  way  significance  can  be  demonstrated  is  making 
a  model  of  half  of  a  record  and  then  using  it  to  predict  the  independent  part  of  the 
record.  This  becomes  a  matter  of  philosophy  as  to  which  way  to  proceed.  At  any  rate, 
the  same  study  as  described  above  has  been  carried  out  over  the  30-year  data  set, 
1947  through  1976.  The  results  are  substantially  the  same. 

5.  The  results  indicated  above  are  for  only  one  season  and  one  lag  forecasts. 
Perhaps  we  will  be  more  successful  in  other  seasons  or  at  other  lags. 

6.  Other  predictands,  e.g.,  those  from  the  Atlantic,  may  increase  our  abilities 
to  predict  air  temperatures  and  other  variables  over  North  America.  Certainly  during 
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other  seasons  there  is  no  reason  to  believe  the  Pacific  is  the  sole  source  of  all  weather 
over  North  America. 

7.  The  predictor  data  has  been  subjected  to  a  variety  of  error  checks.  However 
there  are  some  checks  which  it  still  must  undergo  and  so  we  must  assume  that  its  quality 
is  only  reasonably  well  established.  Thus  unknown  errors  in  the  predictor  data  may  be 
giving  us  lower  levels  of  predictability  than  actually  exist. 
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A  Sea  Level  Pressure  Analog  Model  for  Predicting  Temperature 
Over  the  Northern  Hemisphere 

J.  G.  Schofield 
USAFETAC/DNC  Scott  AFB,  II.  62225 

T.  P.  Barnett 
Scripps  Institution  of  Oceanography,  La  Jolla,  Ca.  92093 

Barnett  and  Preisendorfer  introduced  a  new  scheme  for  using  analog  methods 
in  predicting  climatic  fluctuations  at  the  October  1977  Climate  Diagnostics 
Workshop.  Their  scheme  involves  the  use  of  empirical  orthogonal  functions  to 
represent  data  fields  and  defines  a  "climate  state  vector"  that  is  described 
by  the  eigenvectors  of  the  Northern  Hemisphere  climate  field.  Climatic  analogs 
are  chosen  based  on  the  behavior  of  this  "climate  state  vector." 

It  was  shown  that  the  climate  state  vector  idea  resulted  in  statistically 
significant  forecasting  capabilities.  However,  the  previous  work  did  not  make 
it  clear  which  of  the  fields  in  the  climate  state  vector  was  responsible  for 
the  predictability  nor  indicate  what  regions  of  the  globe  were  the  most  sensi- 
tive predictors.  In  the  present  work  we  have  tried  to  home  in  on  the  fields 
and  geographic  regions  that  might  be  responsible  for  predictability.  As  our 
first  start  we  have 'considered  only  one  field  from  which  to  make  the  analog 
forecast.  That  field  was  the  hemispheric  sea  level  pressure  field  extending 
from  1899  to  present.  We  have  developed  a  sequence  of  models  using  first  the 
entire  hemispheric  field,  then  regional  fields  and  finally  weighted  regional 
fields  in  order  to  determine  the  key  predictor  regions.  Clearly  this  same  type 
of  analysis  will  have  to  be  done  with  other  fields  used  to  define  the  climate 
state  vector,  e.g.,  sea-surface  temperature. 

A  model  based  on  the  hemispheric  pressure  field  was  used  to  predict 
seasonal  surface  temperature  anomalies  for  42  North  American  stations  over  the 
75  year  period,  1902-1976.  A  similar  analysis  has  been  made  for  81  stations 
in  Eurasia.  A  regional  model  that  uses  spatially  limited  (150°E  to  20°W) 
pressure  field  information  was  also  developed  for  predicting  temperature 
anomalies  over  the  North  American  continent  only.  Finally,  an  attempt  was  made 
to  find  an  "optimum"  analog  selection  scheme. 

A  skill  score  was  computed  for  the  models  by  counting  the  number  of  correct 
predictions  at  each  station.  The  temperature  anomalies  for  each  station 
[F(x,t)]  were  grouped  into  terciles,  or  groups  of  upper,  middle,  and  lower 
thirds,  denoted  by  Above  normal,  Normal,  and  Below  normal.  The  two  fields 
F(x,t)  and  F(x,t')  match  at  station  x  if  their  numerical  values  fall  into  a 
common  tercile.  The  total  number  of  matches  for  a  fixed  forecast  time  t  over 
each  station  (x)  gives  a  measure  of  the  skill  of  the  model  at  forecast  time  t. 
The  global  skill  is  the  average  skill  score  over  x  for  all  forecast  times  t. 
The  total  number  of  matches  for  fixed  x  and  all  forecast  times  t  gives  the 
temporal,  or  local  skill  at  each  station  x.  Probability  limits  for  the  local 
skill  scores  are  computed  using  the  binomial  distribution 
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P(r)  =  TTTsWr  ar(l-a)m-r 

where  r  is  the  number  of  matches  and  m  is  the  number  of  trials  (75  in  our 
example).  A  score  of  29  matches  over  75  years  (38.7%)  is  significant  at  the 
92%  probability  level. 

Figure  1  shows  the  local  skill  scores  for  the  hemispheric  model.  Shaded 
areas  represent  significant  local  skill.  A  comparison  with  the  North  American 
model  (Figure  2)  shows  an  increase  in  the  areas  of  significant  skill.  For 
both  models,  the  global  skill  score  is  relatively  low.  The  fact  that  signifi- 
cant local  skill  can  be  obtained  in  certain  areas,  such  as  the  East  Coast,  but 
not  in  others  (the  North  Central  U.S.  for  example)  suggests  that  global  skill 
scores  may  not  be  entirely  representative  of  the  model's  capabilities. 

An  attempt  was  made  to  find  an  "optimum"  method  for  selecting  regional 
analogs.  This  scheme  is  based  on  the  correlation  between  the  sea  level  pressure 
EOF  amplitudes  (modes  1-6)  and  the  North  American  temperature  field  EOF  ampli- 
tudes (modes  1-6).  Figure  3  shows  the  relationship  between  the  winter  pressure 
EOF  amplitudes  and  the  summer  temperature  EOF  amplitudes.  In  this  example, 
we  selected  modes  1,  2,  and  4  for  our  "optimum"  analog  scheme  and  weighted  the 
amplitudes  according  to  their  relative  contributions  in  the  pressure-temperature 
correlation.  Mode  1,  for  example,  is  weighted  more  heavily  than  modes  2  and  4. 

The  results  of  the  "optimum"  selection  model  are  depicted  in  Figure  4. 
Areas  of  significant  local  skill  appear  similar  to  the  results  obtained  by  the 
North  American  model  (Figure  2),  however  the  number  of  stations  with  significant 
local  skill  is  slightly  lower.  The  global  skill  score  for  the  "optimum"  analog 
scheme  is  considerably  lower  than  the  score  for  the  North  American  model. 

In  summary,  it  was  found  that  the  regional  analog  model  for  North  America 
showed  a  significant  increase  in  skill  over  the  Northern  Hemisphere  model.  The 
attempt  to  find  an  "optimum"  analog  scheme  did  not  meet  with  success,  however 
research  in  this  area  is  continuing. 
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Figure  3 
"OPTIMUM"  Analog  Scheme 

Based  on  correlation  of  pressure  EOF  amplitudes  with  temperature 
EOF  amplitudes: 


Winter 

Pressure 

EOF 

Amp. 

(modes  1-6) 


Summer  Temperature 

EOF  Amp.   (modes  1-6)  -*■ 

(.23)  (.25)  -.04 

(.29)   .15   -.14 

-.11  (-.28) 

.07   -.06  (-.29) 

-.02 

-.01    .18    .04 

-.11   (-.32)   .18 

.17    .08 

-.02 

Use  pressure  EOF  modes  1,  2,  4  in  analog  model. 

The  optimum  metric  is 

d2(t,f)  =  W]  (A^O-A^f))2  +w2  (A2(t)-A2(f))2  + 

w3  (A4(t)-A4(t'))2 

min  [d]  at  t,t*    so  t*+6t  is  forecast 

for  t+5t 


w1  =  0.6979 
w2  =  0.2337 
w3  =  0.0683 


34-5 


1 

* 

•»    *    * 


) 


*  i 

*|  i        * 

«i     j  •    * 
*i  *  *>. «      * 

*,  i  *    *•  *  + 

«1  !  I       *  * 

i  *  •*'  i         + 

!*• 

!  *  I  * 

i  j      * 

!'*  !        •: 

!         i  *■  I  * 

!  1*1* 


i  * 

■■  * 
I  * 


♦  .  « 


O! 

I 

I 


-* 
m 


o 


Solar  Variations  and  Weather/Climate 

A.  B.  Pittock 

Laboratory  of  Tree-Ring  Research,  University  of  Arizona,  Tucson  85721 
(On  leave  from  CSIRO  Atmospheric  Physics,  Victoria  3195,  Australia) 

In  a  recent  more  extensive  review  paper  (Pittock  1978)  I  concluded 
that  despite  a  massive  literature  on  the  subject,  there  is  at  present 
little  or  no  convincing  evidence  of  statistically  significant  or  prac- 
tically useful  correlations  between  sunspot  cycles  (11-  or  22-yr. 
periods)  and  the  weather  or  climate.  Further,  that  the  evidence  to 
hand  suggests  that  if  in  the  future  more  data  and  better  analyses 
enable  the  detection  of  statistically  significant  relationships,  these 
will  account  for  so  little  of  the  total  variance  in  the  climatic  record 
as  to  be  of  little  practical  value. 

It  is  conceivable  that  a  more  complex  hypothesis  could  be  developed 
capable  of  accounting  more  subtly  for  a  greater  proportion  of  the  vari- 
ance, perhaps  one  which  allows  for  climatic  discontinuities  or  flips 
from  one  climatic  state  to  another  (suggested  by  data  such  as  in  figs. 
1  and  2)  such  that  solar  influences  have  different  effects  in  different 
epochs.   Such  a  hypothesis  would  not  be  established  as  a  viable  theory, 
however,  until  it  was  used  to  make  detailed  predictions  which  were  ver- 
ified on  independent  data. 

On  shorter  time  scales,  characterized  by  the  period  of  solar  ro- 
tation, the  passage  of  solar  magnetic  sector  boundaries,  solar  flares, 
and  related  events,  there  is  again  a  voluminous  and  confusing  litera- 
ture. Many  claimed  relationships  are  again  based  on  inadequate  data 
and/or  analysis  but  there  appear  to  be  some  statistically  significant 
results  (Wilcox  et  al,  1976;  Hines  and  Halevy,  1977).  Other  related 
studies,  however,  have  drawn  negative  or  inconclusive  results,  so  that 
in  the  absence  of  detailed  and  testable  physical  hypotheses  as  to 
causal  mechanisms  the  few  statistically  convincing  results  are  not  yet 
so  persuasive  as  to  enable  us  to  dogmatically  assert  that  solar-weather 
relationships  on  this  time  scale  have  been  "proven"  to  exist. 

Hope  of  progress  is  greater  on  the  shorter  time  scale  because  it 
is  much  easier  to  gather  independent  data  sets  with  which  to  establish 
more  credible  statistics  and  to  test  mechanisms  and  predictions.   So 
far  statistics  on  the  skill  of  conventionally-based  numerical  weather 
forecasts  following  key  days  (Larsen  and  Kelley,  1977)  show  only  a  mar- 
ginally significant  drop  in  skill  on  days  which  might  be  related  to 
solar- induced  effects.  More  extensive  testing  might  prove  profitable 
here.   Short-term  transient  effects  seem  to  me  intuitively  more  plaus- 
ible than  significant  effects  on  the  11-  and  22-yr.  time  scales. 

On  the  longer  time  scales  of  centuries,  where  the  possibility  of 
physically  significant  changes  in  the  solar  "constant"  exist,  there  is 
some  evidence  suggestive  of  a  relationship  between  solar  variations 
and  climate  (Eddy,  1977) .  Long  proxy  records  of  solar  activity  via 
carbon  14  in  tree  rings,  and  of  climatic  variables,  exist  or  can  be 
developed  for  some  localities.  These  may  lead  to  the  establishment 
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of  statistically  significant  relationships,  but  there  is  a  danger  of 
circularity  in  relating  solar  variations  via  C-14  variations  to  climate, 
since  climatic  change  might  affect  C-14  concentrations. 

The  literature  on  solar  cycle-weather/ climate  relationships  dis- 
plays a  great  lack  of  rigor  and  particularly  of  the  "common-sense" 
application  of  relatively  simple  statistical  considerations.   The 
effects  of  auto-correlations,  smoothing,  and  the  a  posteriori  selection 
of  data,  on  the  statistical  significance  of  results  have  very  often 
been  overlooked.   Elaboration  of  hypotheses  to  fit  the  data,  without 
testing  on  independent  data,  is  not  uncommon.   Contradictory  and  some- 
times incredible  results  have  too  often  been  uncritically  accepted,  and 
there  has  been  a  widespread  lack  of  appreciation  of  the  non-stationary 
behavior  of  climatic  data  series. 

Many  of  these  faults  are  evidently  the  result  of  over-enthusiasm 
and  self-delusion,  and  the  lack  of  inter-disciplinary  communication 
and  criticism.   There  is  evidence  (Gani,  1975)  that  these  faults  are 
widespread  in  the  climatological  literature  where  a  frustratingly  inad- 
equate data  base  often  raises  the  temptation  to  overstate  the  statisti- 
cal or  practical  significance  of  a  result.   This  may  be  particularly 
the  case  with  unrepresentative  or  highly  auto-correlated  oceanic  data 
and  with  short  runs  of  upper  atmospheric  data.   Comments  and  some  gra- 
tuitous advice  will  be  found  in  my  published  review. 
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Captions 

Figure  1.   Water  level  in  Lake  Victoria  1899-1973  as  measured  at  the 
Jinja  gauge  (height  in  meters).   Note  the  apparent  nonstationarity  of 
the  data  set,  and  see  the  discussion  by  Rodhe  and  Virji  (1976). 


Figure  2.   Variations  of  the  correlation  coefficient,  calculated  over 
moving  15-year  intervals,  between  the  mean  sea  level  pressure  at  the 
indicated  stations  in  August  and  the  mean  annual  sunspot  relative 
number.   In  each  case  the  15-year  correlation  coefficient  is  for  the 
15  years  following  the  indicated  date  (after  Bell,  1977) .   Note  that 
the  correlation  coefficients  vary  from  +0.5  to  -0.5  over  different 
15-year  intervals. 
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Additional  comment 

It  may  be  worthwhile  in  the  context  of  short-term  climatic  studies 
(which  have  a  potential  for  forecasting  which  could  be  of  great  econom- 
ic and  social  value)  to  remind  ourselves  of  the  traditional  scientific 
attitude  that  a  hypothesis  is  merely  a  hypothesis  until  it  is  proven 
beyond  reasonable  doubt.   In  contrast  to  this,  there  are  undoubtedly 
situations  of  a  practical  nature  involving  decision-making  under  uncer- 
tainty in  which  one  should  act  on  a  hypothesis  as  if  it  were  true 
unless  it  is  proven  beyond  reasonable  doubt  that  it  is  false.   Such 
decisions  are  properly  made  by  estimating  the  consequences  (or  value) 
of  a  decision  either  way  and  weighting  these  by  the  statistical  prob- 
abilities of  what  statisticians  call  a  type  I  or  type  II  error  (see 
eg.,  Pittock,  1972).   I  suspect  that  some  solar-cycle  enthusiasts, 
convinced  of  the  potential  value  of  a  solar  activity-climate  link  in 
forecasting,  have  confused  type  I  and  type  II  errors  in  their  eager- 
ness to  see  the  practical  application  of  the  solar  activity  hypothesis. 
As  climatologists  in  general  address  themselves  increasingly  to  the 
problem  of  climatic  forecasting,  our  enthusiasm  may  well  tempt  us  into 
a  similar  confusion.   For  the  sake  of  our  science  the  distinction  must 
be  consciously  and  rigorously  maintained. 


Additional  reference: 

Pittock,  A.  B.,  Evaluating  the  risk  to  society  from  the  SST:  some 
thoughts  occasioned  by  the  AAS  report,  Search  3,  285-289,  1972. 
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Status  Report  on  the  Geophysical  Data  Processing  of  the  Seasat  Microwave  Sensors 

John  A.  Ernst  and  John  W.  Sherman,  III  -  NOAA/NESS 

The  Seasat  spacecraft  failed  on  the  night  of  October 
9,  1978  after  approximately  100  days  in  orbit.   During 
this  period  the  four  microwave  sensors,  consisting  of 
an  imaging  radar,  a  scatterometer ,  a  radar  altimeter 
and  a  multichannel  radiometer,  collected  data  which 
will  satisfy  approximately  80%  of  the  NOAA  program. 
objectives . 

This  status  report  represents  only  the  initial  findings 
on  the  sensor's  ability  to  provide  geophysical  information 
on  the  ocean  surface  conditions.   These  data  are  essen- 
tially unvalidated  at  this  juncture,  although  there  are 
comparisons  to  surface  data  in  several  cases. 

Hurricane  Fico,  which  occurred  in  the  mid-Pacific  during 
July  1978,  provided  wind  and  wave  extremes  for  early 
assessment  of  the  scatterometer  and  radiometer,  although 
the  coverage  was  not  simultaneous.   The  early  results 
are  consistent  with  the  expectations  of  these  instruments, 
with  the  initial  scatterometer  results  more  quantitative 
than  the  radiometer  results.   For  the  scatterometer  pass 
over  Fico  on  July  20,  1978,  the  winds  ranged  from  several 
meters  per  sec  to  as  great  as  about  50  m  sec~l.  The 
Seasat  data  is  shown  in  the  sequenced  figures  attached. 
Subject  to  the  atmospheric  attenuation  effects  near  the 
eye  of  Fico,  the  winds  agree  relatively  well  with  surface 
measurements.   The  importance  of  making  the  atmospheric 
corrections  are  obvious . 

The  radiometer  data  has  been  processed  for  Hurricane 
Fico  for  July  12,  1978  in  terms  of  antenna  temperatures, 
not  actual  brightness  temperatures  or  some  other  geo- 
physical unit  (wind  speed,  surface  temperature,  precip- 
itation, etc.)   Even  with  the  modest  resolution  of  this 
instrument,  details  of  the  eye  can  be  readily  seen  and  can 
be  detected  down  to  resolutions  on  the  order  of  50-60  km. 

Data  samples  using  the  imaging  radar  information  from 
Fico  were  processed  from  about  10  km  x  10  km  areas  to 
establish  the  radar's  ability  to  image  waves.   Such 
images  are  possible  when  there  is  a  radial  wave  component 
of  velocity  with  respect  to  the  spaceborne  radar.   Under 
certain  conditions,  wave  imaging  is  difficult  when  only 
an  along-track  wave  component  of  velocity  exists.   The 
selected  sample  illustrated  shows  a  successful  case  and 
the  effects  of  special  processing  techniques. 

Also  shown  for  the  radar  is  a  sample  showing  internal 
waves  off  the  Baja  California  Coast,  the  effects  of  a 
storm  squall  moving  off  the  Florida  Peninsula  and  sea 

36-1 


ice  characteristics  west  of  Banks  Island. 

Sea  ice  from  Antarctica  has  also  been  imaged  by  the 
radiometer  and  a  selected  sample,  again  using  antenna 
temperatures  only,  is  illustrated. 

Lastly,  the  radar  altimeter  is  designed  to  measure  pre- 
cisely the  displacement  from  the  satellite  to  the  ocean 
surface.   GEOS-3  earlier  had  demonstrated  the  geodetic 
utility  of  such  measurements  with  an  rms  precision  of 
about  50  to  100  cm  for  its  radar  altimeter.   The  Seasat 
altimeter  is  expected  to  perform  at  a  precision  approaching 
10  cm, and  based  on  early  noise  analyses  and  the  classic 
history  of  the  Puerto  Rica  Trench  results  from  SKYLAB, 
and  GEOS-3,  the  Seasat  altimeter  appears  intuitively  to 
be  near  its  goal.   The  set  of  three  paired  graphs  compare 
three  different  algorithms  to  obtain  significant  wave 
height  (SWH)  from  GEOS-3  altimetric  data. 
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